[fanxing39]

clear
[u,s]=meshgrid(0:.1:5, -10:.4:10);
z=1;
sech=2*u./(1+u.^2);tanh=(u.^2-1)./(u.^2+1);
R1=0.5*pi^0.25*sech.^0.5;
R2=exp(-s.^2./(2*(1+u.^2)));
R3=exp(z^2/2*tanh+z./(sqrt(2)*u).*sech);
R4=1+tanh-(z-s./(sqrt(2)*u)).^2.*sech.^2;
R=R1.*R2.*R3.*R4;
surf(R)

这样写容易调试