求大神帮帮,我的matlab程序为什么运行不了?syms x1 x y z k;N0=2000;lamda=sqrt(8);alpha0=0.42;alpha2=(2*(1-alpha0))/pi;alpha3=(9*(1-alpha0))/5;bata1=6*(11-2*log(2))/7/(5*pi-18*(1-alpha0).^2);bata2=alpha2*alpha3*bata1;alpha4=2*(alpha0-1/12
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/17 08:18:36
求大神帮帮,我的matlab程序为什么运行不了?syms x1 x y z k;N0=2000;lamda=sqrt(8);alpha0=0.42;alpha2=(2*(1-alpha0))/pi;alpha3=(9*(1-alpha0))/5;bata1=6*(11-2*log(2))/7/(5*pi-18*(1-alpha0).^2);bata2=alpha2*alpha3*bata1;alpha4=2*(alpha0-1/12
求大神帮帮,我的matlab程序为什么运行不了?
syms x1 x y z k;
N0=2000;
lamda=sqrt(8);
alpha0=0.42;
alpha2=(2*(1-alpha0))/pi;
alpha3=(9*(1-alpha0))/5;
bata1=6*(11-2*log(2))/7/(5*pi-18*(1-alpha0).^2);
bata2=alpha2*alpha3*bata1;
alpha4=2*(alpha0-1/12)/pi;
alpha5=75/32*sqrt(2.5)*pi^2.5*alpha4;
bata3=alpha4*alpha5/(1-2*alpha0/pi);
ebcsx=alpha0-alpha2*atan(alpha3*x1*(bata1-x1)/(bata2-x1));%x1^1.5*erfc(sqrt(x-y)))*(2/z.^2+1/k.^2)*(2/3*y*((y-2)*0.5*betainc(1-y/x,0.5,y-1)-(y+1)*0.5*betainc(1-y/x,0.5,y))+((x/3-y/3-0.5)*sqrt(1-y/x)*exp(y-x)-0.25*sqrt(pi/x*erfc(sqrt(x-y)))));
E2=1/3/(betainc(1-y/x,0.5,y+1)+(y/x).^y*exp(x-y)*sqrt(pi)/2/x.^0.5*erfc(sqrt(x-y)))*(2y.^2+landa*k.^2)*(0.5*betainc(1-y/x,0.5,y+1)+(y/x).^y*exp(x-y)/4/x.^2*((3+2*x-2*y)*sqrt(1-y/x)*exp(y-x)+0.5*sqrt(pi/x)*erfc(sqrt(x-y))));
E3=2*N0.^r*3/5*2.^(-1.5*r)*(3*lamda*N0).^(1/3-r/2)*ebcsx*(3*pi).^r/(2*pi*z.^2*k*beta(1-x/y,1.5,y+1)+(y/x).^y*exp(-x+y)+(y/x).^y*exp(x-y)*sqrt(pi)*erfc(sqrt(x-y)/2/x.^1.5))*(0.5*beta(1-y/x,1.5,y*(r+1)+1)+(y/x).^y*exp(x-y)*sqrt(1-y/x)/2/x/(r+1)*exp((r+1)*(y-x))+(y/x).^y*exp(x-y)/4/x/(r+1)*sqrt(pi/x/(r+1))*erfc(sqrt((r+1)*(x-y))));
Etotal=E1+E2+E3;
D1=diff(Etotal,x);
D2=diff(Etotal,y);
D3=diff(Etotal,z);
D4=diff(Etotal,k);
求大神帮帮,我的matlab程序为什么运行不了?syms x1 x y z k;N0=2000;lamda=sqrt(8);alpha0=0.42;alpha2=(2*(1-alpha0))/pi;alpha3=(9*(1-alpha0))/5;bata1=6*(11-2*log(2))/7/(5*pi-18*(1-alpha0).^2);bata2=alpha2*alpha3*bata1;alpha4=2*(alpha0-1/12
难道以前是能够运行的,现在不能了 ,还是一直都不能,调试嘛