我有这个 Mathematica 问题,它在一个盒子中以数值方式解决了与时间相关的薛定谔方程:
sig = 1;
bsize = 655;
s = NDSolve[{D[u[t, x], t] == I/2*D[u[t, x], x, x],
u[0, x] == (1/(2*Pi*sig^2))^(1/4)*Exp[-x^2/(2*sig)^2],
u[t, -bsize] == u[t, bsize]},
u, {t, 0, 10}, {x, -bsize, bsize}];
wx[t_, x_] := u[t, x] /. s
px[t_, x_] := Abs[wx[t, x]]^2
Manipulate[Plot[px[t, x], {x, -10, 10}, PlotRange -> {{-10, 10}, {0, 1}}], {t,0, 10}]
的解是:
但是,如果我将盒子大小bsize增加到bsize=656我得到这个不连贯的解决方案:
问题:这里发生了什么?

