Numerical solution of the 1D Gross-Pitaevsky equation with Crank-Nicolson alogrithm. Unitarity respected thanks to Cayley's form.