MATLAB Mini Hack

Back to Gallery
  • New Entry

Atomic Orbital of Hydrogen (4,3,0)

Daniel Pereira on 27 Oct 2021
  • 6
  • 16
  • 3
  • 0
  • 280
N=4;
L=3;
M=0;
r=meshgrid(0:.5:50);
t=r'/25*pi;
a=N-L-1;
% Short version of laguerre polynomial
F=0;
for i=0:a
F=F+(-1)^i*nchoosek(N+L,a-i)*(2*r/N).^i/prod(1:i);
end
G=legendre(L,cos(t));
P=squeeze(G(1+M,:,:)).*F.*exp(-r/2).*r.^L;
contourf(r.*sin(t),r.*cos(t),abs(P),1e3,'edgec','n');
caxis([0 5]);
axis equal off;
axis(35*[-1 1 -1 1]);
colormap hot;
Remix Tree
Load full remix tree