>>119
ttp://www-solid.eps.s.u-tokyo.ac.jp/~ataru/edu/ensyu080609.pdf

ここらへん？

それにしても酷過ぎて追うのやめちまったよ……

k2=h*f1(t+h/2.0,x+k1/2.0,vx+l1/2.0,vy+l1/2.0,y+k1/2.0);

なんでlを２回使ってるの？

double f1(double t,double x,double y,double vx,double vy){return(vx);}

なんで(t,x,y,vx,vy)の順番に並んでいないの？

double g1(double t,double x,double y,double vx,double vy){return(-x/(r*r*r));}

rはいつ決まるの？pow( (x*x+y*y), 1.5 ) じゃないのかなー。