验算《光波导技术物理基础》p7的例子
由于计算TE模，使用如下函数文件
%fc.m计算平板波导的模式%
function y=fc(x)
n1=2.2899;n2=1.5040;n3=1.0000;bc=1.06;h=1.5;
m=0; %可使用不同值代入
k0=2*pi/bc;
y=k0*h*sqrt(n1^2-x^2)-(m*pi+(atan(sqrt((x^2-n2^2)/(n1^2-x^2)))) +(atan(sqrt((x^2-n3^2)/(n1^2-x^2)))));
x=[y];
屏幕输入x=[2];
fsolve(‘fc’,x) %二乘拟合
选取初始值N＝2，结果如下：
m 计算得到的N值 书本上的N值
0 2.2681 2.264
1 2.2016 2.201
2 2.0873 2.086
3 1.9192 1.916
4 1.6889 1.685
5 1.4115（小于n2）
