c7impode.m

高等引用数学的matlab求解的配套程序

   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)];