function main()
clc
close all
global c u %定义全局变量
c=3;
u=7;
[T,Y]=ode45(@fun,[0 100],[3 2 12]); %根据fun定义的微分方程组，使用ode45函数求解微分方程组
%绘制x随时间t变化曲线
figure
plot(T,Y(:,1)) %绘制曲线
xlabel('t') %添加横轴标注t
ylabel('x') %添加纵轴标注x
title('x的值') %添加标题
%绘制y随时间t变化曲线
figure
plot(T,Y(:,2)) %绘制曲线
xlabel('t') %添加横轴标注t
ylabel('y') %添加纵轴标注y
title('y的值') %添加标题
%绘制z随时间t变化曲线
figure
plot(T,Y(:,3)) %绘制曲线
xlabel('t') %添加横轴标注t
ylabel('z') %添加纵轴标注z
title('z的值') %添加标题
figure
plot3(Y(:,1),Y(:,2),Y(:,3)) %plot3绘制三维曲线
xlabel('x')
ylabel('y')
zlabel('z')
title('x y z 三维曲线')
msgbox('MATLAB编程答疑，请加QQ: 1530497909','MATLAB答疑','help')
web http://url.cn/TKcdXk -browser
end
function dy=fun(t,y)
%fun定义微分方程组
global c u %全局变量，使用之前需要先声明
dy=zeros(3,1);
dy(1)=y(2)-y(1);
dy(2)=y(1)*y(3)+c;
dy(3)=u-exp(y(1)*y(2));
end
