ff5.m

来自「研究生期间的数值计算的书配套程序源代码 奉献了」· M 代码 · 共 2 行

    [T,Y]=ode45('rigid',[0 12],[0 1 1]);
     plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')