fvd.m

利用龙格库塔法（等步长

function rhs=fvd(t,u)

rhs=[(1*10^(-14)*(exp((u(1)-5/sqrt(2))/0.0259)-1)+0.01*u(1))/(sqrt(-1)*2*pi*50*0.001)                       
    ];