% normh.m % %% return the inverse of the square root of the %% normalization integral int(R(z)^2 z^2 d z) %% radhyd0(n,l,'z') is the unnormalized radial function %% function str2 = normh(n,l) str1=radhyd0(n,l,'zz'); str1=symmul(str1,str1); str1=symmul(str1,'zz^2'); str1=simplify(str1); str1=simplify(expand(str1)); str1=int(str1,'zz'); str2=subs(str1,0,'zz'); str2=['-(' str2 ')']; str2=simplify(str2); str2=sympow(str2,-1/2); str2=simplify(str2);