func6p3.asv

来自「数值分析最常用的四十种算法」· ASV 代码 · 共 3 行

options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode23(@sub6p3,[0 12],[0 0 0],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')