fvdp2.m

来自「利用龙格库塔法（等步长」· M 代码 · 共 4 行

function rhs=fvdp2(t,u)
R=2
rhs=[1-(R^2*(u(1)^2-u(2)^2))/(u(1)^2+u(2)^2)^2;
    (-2*u(1)*u(2)*R^2)/(u(1)^2+u(2)^2)^2];