fvdp.m

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

function rhs=fvdp(t,u)
global epsilon
rhs=[u(2);
    -epsilon*(u(1)*u(1)-1)*u(2)-u(1)];