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

HarmonicTrap1D[];
sizes={1024};
timeStep=0.001;
k=-0.4;
contactInteractionFactor=4 Pi k;
max=10ho;
min=-10ho;
initialize[];
gridReport

calcAllTable
p1=plotProjections

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

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

reportMoments

Turn off harmonic potential

newPotential[free];
plot:=AppendTo[lplots,{time,plotProjections}];
reset[];
myPrint=Null&;
AbsoluteTiming[evolve["rte",40,200,200]]

plotM["sigmaX"]

Show[lplots[[All,2]],PlotRange->All]