
S=RK2(f,300,0.05)
scf(0)
clf()
plot2d([0:0.05:300],S(1,:))

