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（