% radhyd0.m % %% return the unnormalized radial function as a string % procedure rad(n,l,r); sub(x=2*r/n,lagguere(n+l,2*l+1,x)); function str2 = radhyd0(n,l,r) str1=[ '2*' r '/' num2str(n) ]; % The string for subs; str1=simplify(str1); % rad=lagguere(n+l,2*l+1,'x'); rad=subs(rad, str1, 'x'); % % procedure radhyd0(n,l,r); rad(n,l,r)*r^l*e^(-r/n); % str2=[ r '^' num2str(l) '*exp(-' r '/' num2str(n) ')' ]; str2=simplify(str2); str2=[ '(' rad ')*' str2]; str2=simplify(str2);