Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- [Psi] = NDSolve[{I D[u[t, x], t] == -D[u[t, x], {x, 2}] +
- I Sin[x] u[t, x], u[0, x] == Exp[-x^2] Exp[I x], u[t, -Pi] == u[t, Pi]}, u, {t, 0, 5}, {x, -Pi, Pi}]
- DensityPlot[Evaluate[Abs[u[x, t]] /. [Psi]], {x, -Pi, Pi}, {t, 0, 5}, PlotPoints -> 200, Mesh -> False]
Add Comment
Please, Sign In to add comment