matlab中求大神帮助，不能运行出结果，总出现下图，没有图像，求解，谢谢

Question

n=1.5;T=((1*1.55)/(2*n))*10^(-6);L=20*10^(-3);n1=5*10^(-5);wave0=1.5498*10^(-6);wave=wave0:0.1:1.5502*10^(-6);w=(pi*n1)/wave;beita3=2*n*pi/wave-pi/T;S=abs(sqrt(w^2-(beita3)^2));r=w^2*(sinh(S*L))^2/((beita3)^2*(sinh(S*L))^2+S^2*(cosh(S*L))^2);axis([1.5498*10^(-6) 1.5502*10^(-6) 0 1]);plot(wave,r,'b-');

Answer 1

n = 1.5;
T = ((1*1.55)/(2*n))*10^(-6);
L = 20*10^(-3);
n1 = 5*10^(-5);
wave0 = 1.5498*10^(-6);
wave = wave0:0.001*10^(-8):1.5502*10^(-6);  % 间隔太大，取小点
w = (pi*n1)./wave;
beita3 = 2*n*pi./wave-pi/T;                 % 所有向量的乘除要用点乘（.*）、点除（./）
S = abs(sqrt(w.^2-(beita3).^2));            % 所有向量的开方（.^） 
r = w.^2.*(sinh(S.*L)).^2./((beita3).^2.*(sinh(S.*L)).^2+S.^2.*(cosh(S.*L)).^2);
plot(wave,r,'b-');
axis([1.5498*10^(-6) 1.5502*10^(-6) 0 1]);