GrossPitaevskiiEquation Fortran program that solves the three-dimensional Gross - Pitaevskii equation for a system of N bosons confined in an harmonic trap. $ E = mc^2 $