NotebookEvaluate[NotebookDirectory[]<>"/split_step.nb"];

HarmonicTrap3D[0.2,1,1];
sizes={64,16,16};
timeStep=0.001;
k=-0.3;
contactInteractionFactor=4 Pi k;
max=6ho;
min=-max;
myPrint=Null&;
plotProjections:=plotProjectionsOpts[Filling->0];
initialize[];
gridReport

calcAllTable
p1=plotProjections;

AbsoluteTiming[evolve["ite",10,10]]

calcAllTable
p2=plotProjections;
Show[p1,p2,PlotRange->{{-6,6},All}]

psi[0.2,0.3,0.1]

Plot[Abs[psi[x,0,0]],{x,-12,12}]

psi0=psi;
x0=1;
psi1[{x_,y_,z_}]=Piecewise[{{psi0[x-x0,y,z],Abs[x-x0]<12}},0];
Plot[Abs[psi1[{x,0,0}]],{x,min[[1]],max[[1]]},PlotRange->All,Filling->0]
wf=Map[psi1,mesh[],{dim}];

meanX

reset;
k=-0.5;
contactInteractionFactor=4 Pi k;
AbsoluteTiming[evolve["rte",50,500,500]]

plotM["meanX"]

m1=listM["meanX"];
m2=listM["sigmaX"];
c2=m2^2-m1^2;
c2[[All,1]]=m1[[All,1]];
c2[[All,2]]=Sqrt[c2[[All,2]]];
ListPlot[c2,Joined->True,PlotRange->All]

lplots[[1;;;;50,2]]