schro_p.m

clear;
load x.dat;  load psip.dat;  load p_plot.dat
temp = size(p_plot);
plot(x,psip);
title('Hit any key to continue');
xlabel('Position');  ylabel('real(psi) and imag(psi)');
pause;
plot(x,p_plot(:,1:3:temp(2)-1),x,p_plot(:,temp(2)));
xlabel('position'); ylabel('Prob. density');
title('Probability density at various times');
pause;
mesh(p_plot,[90 30]);
xlabel('position'); ylabel('time');