untitled2.m

来自「简单的解常微分方程的matlab程序！压缩包」· M 代码 · 共 6 行

function y=workci(T);
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,y] = ode45(@rigid,[0 12],[0 1 1],options);
[m,n]=size(y);
cn=y(m,:)
plot(T,y(:,1),'-',T,y(:,2),'-.',T,y(:,3),'.')