hydrowavefnplots1.m

物理公式代码库 常见的物理公式可以从库中直接调用 方便

% Plot the first few radial wavefunctions for the hydrogen atom with l=1.

x=[0:.05:10];
y = zeros(5,201);

for n=1:5
    y(n,:) = -x .* exp(-x) .* polyval(AssociatedLaguerrePoly(n,1),x);
end

figure;
plot(x,y(1,:),x,y(2,:),x,y(3,:),x,y(4,:),x,y(5,:)),
    xlabel('\rho'),...
    ylabel('R_{n,1}(\rho)'),title('Hydrogen Radial Wavefunctions; n\leq 5, l=1')
h = legend('n=1', 'n=2', 'n=3', 'n=4', 'n=5', 4);
grid on;