
T1=RK2(f,%pi/6,0.1)
T2=RK2(f,%pi/4,0.1)
T3=RK2(f,%pi/3,0.1)

scf(0)
clf()
plot2d(T1(1,:),T1(2,:))
plot2d(T2(1,:),T2(2,:))
plot2d(T3(1,:),T3(2,:))

c=cible(3000)

T4=RK2(f,c,0.1)
scf(1)
clf()
plot2d(T4(1,:),T4(2,:))
