examp3_33.m

来自「matlab的大量实例代码」· M 代码 · 共 11 行

tic,[t,y]=ode45(@c3xstf1,[0,1],[1;0;-1]); toc

x1=exp(-2*t); x2=exp(-40*t).*cos(40*t); x3=exp(-40*t).*sin(40*t);
y1=[0.5*x1+0.5*x2+0.5*x3, 0.5*x1-0.5*x2-0.5*x3, -x2+x3];
plot(t,y,t,y1), [max(abs(y-y1)),max(diff(t))]

tic,[t2,y2]=rk_4('c3exstf1',[0,1],[1;0;-1]); toc

plot(t,y1,t2,y2)