matlab中的mpower,求助，Matlab程序运行Error using ==> mpower

%%%%%%%已知波长(m) , π       电场(V/m)%%%%%%%%%%%%%%%%%%%%%%%%%%%%

B=632.8*10^-9;  W=3.14;  E0=0.9*10^5;

%%%%%%%%%%变量 %%%%%

T=[0:1:300];

%%%%%%%%%%%%LiNbO3结构参数%%%%%%%%%%%

r=0.025;%%%%%PPLN长度(m)%%%%%

L=0.00075;%%%%%正极长度(m)%%%%

D=0.75;%%%%占空比系数%%%%

M=L/D;%%%%极化周期%%%%

m=2;%%%%%准相位匹配阶数%%%%%

kg=2*W*m/M;%%%周期极化晶体的倒格矢%%%

A10=0 ;A20=1;

Gm=(1-cos(2*W*m*D)+i*sin(2*W*m*D))/(i*W*m);    G0=2*D-1;

k0=2*W/B; %%%%%%%入射光波矢%%%%%

%%%%%%%%%%%%%%%%%%%%%%o光和e光折射率%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

n1=sqrt(4.9130+(1.173*10^5+1.65*10^-2*T^2)/(B^2-(2.12*10^2+2.7*10^-5*T^2)^2)-2.78*10^-8*B^2);

n2=sqrt(4.5567+2.605*10^-7*T^2+(0.970*10^5+2.70*10^-2*T^2)/(B^2-(201+5.4*10^-5*T^2)^2)-2.24*10^-8*B^2);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%