Adjusted for extremely small numbers 10^-101 breaking max. length for output
Plots wave functions with correct axes labels
Took care of symmetry, changed tanTheta routine to check if cot2Theta is neg/positive and implemented rewrite
Changed initial step i=0 -> i=1, and changed h=rhomax/N-1 -> h=rhomax/N accordingly