• Remix
  • Share
  • New Entry

on 27 Oct 2021
  • 2
  • 8
  • 0
  • 0
  • 274
N=4;
L=3;
%M=1;
r=meshgrid(0:.3:50);
t=r'/9*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(1,cos(t));
P=squeeze(G(2,:,:)).*F.*exp(-r/2).*r.^L;
contourf(r.*sin(t),r.*cos(t),abs(P),1e3,'edgec','n');
caxis([0 4]);
axis equal off;
axis(30*[-1 1 -1 1]);
colormap jet;
Remix Tree
Load full remix tree