c7impode.m

来自「薛定宇高等应用数学问题的MATLAB求解例子源码」· M 代码 · 共 6 行

   function dy=c6impode(t,x)
   dx=inline(['[p(1)*sin(x(4))+p(2)^2+2*x(1)*x(3)-x(1)*p(1)*x(4);',...
              'x(1)*p(1)*p(2)+cos(p(2))-3*x(3)*x(2)]'],'p','x');
   ff=optimset; dx1=fsolve(dx,x([1,3]),ff,x);
   dy=[x(2); dx1(1); x(4); dx1(2)];