Question

matlab编程求解一阶微分方程组？我编了程序但是出错了，在线等高手解答。

题目不难。因为我是新手，所以不懂怎么调试，请大家不吝赐教。
function zODE
clear all;clc
X0=[0,898];
[z,X] = ode45(@Equations,0:0.1:3.5,X0)
plot(z,XA,'k-',z,T,'k-')
legend('XA','Location','NorthWest','T','Location','NorthEest')
xlabel('z'),ylabel('XA','T')
%---------------------------------
function dXdz=Equations(t,X)
XA=X(1);T=X(2);
dXAdz=4.104e6/(21+XA)*exp(-10983/T)*(1-XA-3.03e-7*XA^2/(21+XA)*exp(14520/T));
dTdz=0.07605*(1100-T)-137*(4.104e6/(21+XA)*exp(-10983/T)*(1-XA-3.03e-7*XA^2/(21+XA)*exp(14520/T)));
dXdz=[dXAdz;dTdz];
请大家帮忙指正，在线给我留言就更好了！感谢…

Answer 1

你的程序没有很大的错误，只是你的画图的语句错了
clear all;clc
X0=[0,898];
[z,X] = ode45(@Equations,0:0.1:3.5,X0)
plot(z,X(:,1),'k-',z,X(:,2),'k-')
这样写就行了

追问

为什么改了还是错呢？虚心求解…

Answer 2

一行一行的代码执行，如果出错，会自动提示的，这个比C语言好调试多了。既然不会，那就更要学会调试了，这是最基本的知识。