亲们想问一下这个MATLAB脚本
代码怎么错了,一直不太懂:
clc;
clear;
q=1.6*10^(-19); ap=1.582*10^6;
e0=8.85*10^(-12); er=11.9;
A=1.8*10^35; an=7.03*10^5;
bp=2.036*10^6; bn=1.231*10^6;
u=1360;
syms NB t x Wc
k0=-q*NB/(er*e0);
E0=-(k0*Wc);
E=k0*x+E0;
anzhu=7.03*10^5*exp(-bn./(k0*x+E0));
apzhu=bn*exp(-bp./(k0*x+E0));
f=exp(t)/t;
Ei1=int(f,t,-inf,-(bn/E0));
Ei2=int(f,t,-inf,-(bp/E0));
F0=(an*bn*Ei1-ap*bp*Ei2+(an*exp(-bn/E0)-ap*exp(-bp/E0))*E0)/k0;
Ei3=int(f,t,-inf,-bn/(k0*x+E0));
Ei4=int(f,t,-inf,-bp/(k0*x+E0));
MOM=int(anzhu*exp(-ap/k0*(bp*Ei3+E0*exp(Ei4))-F0),x,0,Wc);
eqn=[MOM==1];
var=[Wc];
[Wc]=solve(eqn,var);
Ronsp=Wc./(q*u*NB);
BVpp=int(k0*x+E0,x,0,Wc);
Wc
NB