rigid_example.m

来自「自己编的四节龙格琨塔程序」· M 代码 · 共 4 行

%script.m
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options,10);
plot(T,dy(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')