平板光波导芯层中所有导模的传播常数求解的MATLAB程序

clear
close all
n1 = 1.56;              %芯层折射率
n2 = 1.2;               %衬底折射率
n3 = 1;                 %覆盖层折射率
d = 3e-6;               %芯层厚度
lambda = 1.55e-6;       %光波波长
k = 2*pi/lambda;

V12 = sqrt（n1^2-n2^2）*k*d;
V13 = sqrt（n1^2-n3^2）*k*d;

F = @（x）（x*（sqrt（