examp7_21.m

来自「《高等应用数学问题MATLAB求解》的PDF版及习题解答」· M 代码 · 共 2 行

options=odeset; options.Mass=@c7eqdae2; options.RelTol=1e-6;
      [t,x]=ode45(@c7fimp2,[0,10],[0;0],options); plot(t,x)