clear?all

tic
eps0=8.85e-12;
mu0=4*pi*1e-7;
ll=sqrt（eps0/mu0）;
c=3e+8;
lamd0=800e-9;
i=sqrt（-1）;
w0=2*pi*c/lamd0;

ct0=0*pi/180;

c0=cos（ct0）;
NN=5;

na=1.7;
a=142e-9;
nb=3;%（3:5.5）硅光電池
%nb=1.66;
b=200e-9;
%4*lamb
n0=2.86;
%?k0=0.15;
k0=0;
nd=n0+k0*i;
d=500e-9;
%l=（a+b+d）*NN*2;
l=（a+b）*2*NN+d;

scta=sin（ct0）/na;
sctb=scta*na/nb;
sctd=scta*nb/nd;

ca=sqrt（1-scta*scta）;
cb=sqrt（1-sctb*sctb）;
cd=sqrt（1-sctd*sctd）;

deltaA=（ca*na*a/c）;
deltaB=（cb*nb*b/c）;
deltaD=（cd*nd*d/c）;

????etaA=（ll*na*ca*ca）;
????etaB=（ll*nb*cb*cb）;
????etaD=（ll*nd*cd*cd）;
%?????etaA=（ll*na/ca/ca）;
%?????etaB=（ll*nb/cb/cb）;
%?????etaD=（ll*nd/cd/cd）;


n=1;nw（1）=nb;m=-1;k=0;
dx=0.0001;

shui=[];shuo=[];heni=[];heno=[];dw=[];ww0=[];
p=0;q=0;
for?x=8/11:dx:8/3
????Ma=[cos（deltaA*x*w0）?-i*sin（deltaA*x*w0）/etaA;-i*etaA*sin（deltaA*x*w0）?cos（deltaA*x*w0）];??%%?x=w/w0
????Mb=[cos（deltaB*x*w0）?-i*sin（deltaB*x*w0）/etaB;-i*etaB*sin（deltaB*x*w0）?cos（deltaB*x*w0）];
????Md=[cos（deltaD*x*w0）?-i*sin（deltaD*x*w0）/etaD;-i*etaD*sin（deltaD*x*w0）?cos（deltaD*x*w0）];
????
M=（Mb*Ma）^NN*Md*（Mb*Ma）^NN;

????t=2*ll*c0/（M（11）*ll*c0??+??M（12）*ll*ll*c0*c0??+??M（21）+??ll*M（22））;
????r=（M（11）*ll*c0??+??M（12）*ll*ll*c0??-??M（21）??-??ll*M（22））/（ll*M（11）*c0+??ll*ll*M（12）*c0+M（21）+ll*M（22））;

????%t=2*ll/（M（11）*ll+M（12）*ll*ll+M（21）+ll*M（22））;
????%r=（M（11）*ll+M（12）*ll*ll-M（21）-