求助数值解复数空间超越方程
学生向各位求教复数空间超越方程的数值解
即 已知各参数求满足方程fl=fr的beta_1.
e1=1.53;e2=1;e3=-5.73+0.22*i;d=0.2;wavelength=0.441;k0=2*pi/wavelength
k1=(k0^2*e1-beta_1.^2).^0.5;
alpha2=(beta_1.^2-k0^2*e2).^0.5;
alpha3=(beta_1.^2-k0^2*e3).^0.5;
D1D2=(alpha2./e2)./(k1./e1); 
C2C1=(alpha3./e3)./(k1./e1); 
fl=k1*d;
fr=(atan(D1D2)+atan(C2C1))./k1;