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

HarmonicTrap3D[];
sizes={32,32,32};
max=4ho;
min=-max;
timeStep=0.0001;
k=-0.57;
contactInteractionFactor=4 Pi k;
initialize[];
gridReport

calcAllTable
plotProjections

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

calcAllTable
plotProjections

reportEnergies

reportMoments