/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created by wxMaxima version 0.7.4 ] */

/* [wxMaxima: input   start ] */
load(orthopoly);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
declare([p,q],integer);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
sLag(p,q,u):=diff(laguerre(p,u),u,q);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
N(n,l,Zs):=-sqrt((2*Zs/n)^3*(n-l-1)!/(2*n*(n+l)!));
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
R(n,l,Zs,u):=N(n,l,Zs)*(u)^l*%e^(-u/2)*sLag(n+l,2*l+1,u);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Y(l,m,theta,phi):=spherical_harmonic(l,m,theta,phi);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
psi(n,l,m,Zs,u,theta,phi):=R(n,l,Zs,u)*Y(l,m,theta,phi);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Psi(n,l,m,Zs,r,theta,phi):=subst(2*Zs*r/n,u,psi(n,l,m,Zs,u,theta,phi));
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Zs1:Z/a0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Psi(1,0,0,Zs1,r,theta,phi);
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$