examp7_21.m

来自「matlab例题全部代码」· M 代码 · 共 2 行

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