怎么用编程画氢原子的电子轨道

Question

这些软件都没有学过，希望各位能帮帮忙，不是计算机专业的但是需要用到这个，再次麻烦各位了

Answer 1

不要求归一化的话可以这样画：

subplot(2,2,1)
l=1;m=0;
[Xm,Ym,Zm]=getLegendreSurf(l,m);
surf(Xm,Ym,Zm)

subplot(2,2,2)
l=2;m=0;
[Xm,Ym,Zm]=getLegendreSurf(l,m);
surf(Xm,Ym,Zm)

subplot(2,2,3)
l=2;m=2;
[Xm,Ym,Zm]=getLegendreSurf(l,m);
surf(Xm,Ym,Zm)

subplot(2,2,4)
l=3;m=2;
[Xm,Ym,Zm]=getLegendreSurf(l,m);
surf(Xm,Ym,Zm)



function [Xm,Ym,Zm]=getLegendreSurf(l,m)
dx=pi/30;
col=0:dx:pi;
az=0:dx:2*pi;
[phi,theta]=meshgrid(az,col);

Plm=legendre(l,cos(theta));

if l~=0
    Plm=reshape(Plm(m+1,:,:),size(phi));
end

a=(2*l+1)*factorial(l-m);
b=4*pi*factorial(l+m);
C=sqrt(a/b);
Ylm=C.*Plm.*exp(1i*m*phi);

[Xm,Ym,Zm]=sph2cart(phi,pi/2-theta,abs(real(Ylm)));
end

matlab代码