Solution (source code)

= Solution

Minimal coupling gives $H=(p-qA)^2/(2m)$ and $i\hbar\partial_t\psi=H\psi$.