exc7eqa.m

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

   function dx=exc7eqa(t,x)
   r=sqrt(x(1)^2+x(3)^2);
   dx=[x(2); -x(1)/r^3; x(4); -x(3)/r^3];