``pseudo.m''

load psis_real
psis=psis_real;

%# Plot all-electron atom wave functions
figure(0);
axis([0 5]);

r=psis(:,1);
plot(r,psis(:,2),';1s;',r,psis(:,3),';2s;',r,psis(:,4),';3s;', ...
     r,psis(:,5),';4s;',r,psis(:,6),';2p;',r,psis(:,7),';3p;', ...
     r,psis(:,8),';4p;',r,psis(:,9),';3d;');

%# Compare pseudo and all-electron wave functions
figure(1);
load psis_pseudo
plot(r,psis(:,5),'r-;4s;',r,psis(:,8),'g-;4p;', ...
     r,psis_pseudo(:,2),'rx;1s(ps);',r,psis_pseudo(:,3),'gx;2p(ps);');