fvdp1.m

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

function rhs=fvdp(t,u)

rhs=[u(3);
    u(3)*u(4)
    ];