%所给条件中,由于i、m的值都未给定,所以在处理时,均按=1处理,你可以在后期的
%计算中,自行带入常量,重新生成波函数图像。
%题中,r01是常矢量,计算时按模长带入。r01=1.15453,并用r1代替
% ћ(约化普朗克常数)符号替代为h
h=1.05457266*1e-34;
r1=1.15453;
t=0:0.1:2*pi;
r=0:0.1:2*pi; %r的范围后期自行更改
f=(sqrt((1)/(pi*h)))./(sin(t).*(1-cot(t))).*exp(-1/(2*h)*(r-cot(t)*r1).^2)...
.*exp(sin(t)/(2*h).*(cos(t)*r1^2-2*r1*r));
plot3(r,t,f)
%由于参数i、m等值不争取,所以函数图像难免会有出入,可以自行更改参数
纯手打,望采纳!!