10楼: Originally posted by xds42 at 2011-11-09 16:57:05:
哪位高手帮忙，用matlab算算吧，金币再另送。

function AA

z=0:0.01:0.2;
[z,x]=ode15s(@odefun,z,[0 463]);
plotyy(z,x(:,1),z,x(:,2))
legend('z-Xa','z-T')

function dXT=odefun(z,x)
%x(1)=Xa
%x(2)=T
dXa=1/8.95012e-6*5.94e10*exp(-4632.97/x(2))*(1-x(1))/(1+x(1));
dT=1/0.00133894*(-34680*(-463+x(2))+(5.2272e10-66*exp(-4632.97/x(2))*(1-x(1)))/(1+x(1)));
dXT=[dXa;dT];