Trojan wavepacket viewer (click the green text to open Java application)

Copy and paste the red code below into the upper Java applet supercalculator window; NN=1-Infinity is the quantum number; eps=0-0.11556 is the electric field:

NN=60;
eps=0.04442;
gg=(27+2*eps**3+3*3**0.5*(27+4*eps**3)**0.5)**(1/3);
q=1/(eps/3 + 2**(1/3)*eps**2/gg/3 + gg/3/2**(1/3))**3;
q=(Signum(q-8/9)+1)/2*q-0.9999999*(Signum(q-8/9)-1)/2;
s=(Abs(1-q)*(1+2*q))**0.5;
A=(Abs(1-q)*(8+4*q-9*q**2-8*s))**0.5/3/q*NN;
B=((1+2*q)*(8+4*q-9*q**2-8*s))**0.5/3/q*NN;
Plot3D(Exp(-B*((x**2+y**2)**0.5-q**(-1/3))**2-2*A*(1-x/(x**2+y**2)**0.5)),{x,-1.5,1.5},{y,-1.5,1.5})

(Expressions are according to circular symmetry augmented formula improved to Phys. Rev. Lett. 73, 1777 (1994))