quadeg6.m

来自「数值类综合算法 常用数值计算工具包（龙贝格算法、改进欧拉法、龙格库塔方法、」· M 代码 · 共 23 行

%微积分例6：刚性微分方程组（ode15s)
%需用模型函数quadeg6fun.m
clear;close;
[t,y]=ode45('quadeg6fun',[0,10],[2,1]');
plot(t,y);
text(1,1.1,'y1');
text(1,0.1,'y2');

pause

[t,y]=ode45('quadeg6fun',0,400,[2,1]');
tstep=length(t), minh=min(diff(t)), maxh=max(diff(t))

pause

[t,y]=ode15s('quadeg6fun',[0,400],[2,1]');
plot(t,y);
text(100,0.5,'y1');
text(1,0.1,'y2');
length(t), min(diff(t)),max(diff(t))