clear;
n1=1.615;
d1=0.24;
n2=1.336;
d2=0.29;
no=1.40;
ne=1.45;
d=3;
lambda=1555;
syms theta1 theta2
%
eta1=n1*cos(theta1);
>> eta2=n2*cos(theta2);
>> delta1=2*pi*(1/lambda)*n1*d1*cos(theta1);
>> delta2=2*pi*(1/lambda)*n2*d2*cos(theta2);
>>  m1=[cos(delta1) i*(1/eta1)*sin(delta1)
i*(1/eta1)*sin(delta1) cos(delta1)];
>> m2=[cos(delta2) i*(1/eta2)*sin(delta2)
i*(1/eta2)*sin(delta2) cos(delta2)];
>> %周期单元的
>> N=5;
>> M1=(m2*m1)^N;
>> M2=(m1*m2)^N;
 n=ne-no;
>> B=2*pi*n/lambda;
>> arf=90*pi/180;%
>>  A=sqrt((arf)^2+(B/2)^2);
>> syms epsilon1 epsilon2;
>> gama=epsilon1/epsilon2;
>> k11=11.9;
k33=13.6;
>> k=k33-k11;
epsilon0=1;
>> theta=0:pi;
>> subs(theta,'C',1);
 V=4:12;


diff(theta,'z')=-(epsilon0*(epsilon2+epsilon1*(sin(theta))^2)*V/d*sqrt(epsilon0*epsilon2*k11))*sqrt((C+1/((1+gama*(sin(theta))^2))/(1+k*(sin(theta))^2)));   
[color=FF0000]??? Error using ==> mpower
Matrix must be square.[/color]
J=[cos(arf) -sin(arf);sin(arf) cos(arf)]*[cosA-i*(1/(2A))*B*sinA (1/A)*(arf)*sinA;-(1/A)*(arf)*sinA cosA+i*(1/((2*A))*B*sinA];
??? J=[cos(arf) -sin(arf);sin(arf) cos(arf)]*[cosA-i*(1/(2A))*B*sinA (1/A)*(arf)*sinA;-(1/A)*(arf)*sinA cosA+i*(1/((2*A))*B*sinA];
                                                          |
Error: Missing MATLAB operator.