小推廣的Runge-Kutta方法正確嗎？

同學們：
有一個數學問題請教：“Runge–Kutta methods”
可以參考wikipedia：
http://en.wikipedia.org/wiki/Runge%E2%80%93Kutta_methods
其中有一句話“Now pick a step-size h>0 and define”
事實上，星係的曲線是可以繞圈的，所以上麵的h>0 有時必須h小於0.
我的程序是這樣處理有時大於零，有時小於零（下麵是程序）。我的處理正確嗎？

k1=tan(alp); dec=sin(alp);
if( cos(alp)>=0 ){goe1=fabs(goe1);}else{goe1=-fabs(goe1);}
xh=X+0.5*goe1; yh=Y+0.5*goe1*k1;
sl2=x1v(xh,yh,dec);
k2=tan(sl2); dec=sin(sl2);
if( cos(sl2)>=0 ){goe2=fabs(goe2);}else{goe2=-fabs(goe2);}
xh=X+0.5*goe2; yh=Y+0.5*goe2*k2;
sl2=x1v(xh,yh,dec);
k3=tan(sl2); dec=sin(sl2);
if( cos(sl2)>=0 ){goe3=fabs(goe3);}else{goe3=-fabs(goe3);}
xh=X+goe3; yh=Y+goe3*k3;
sl2=x1v(xh,yh,dec);
k4=tan(sl2); dec=sin(sl2);
if( cos(sl2)>=0 ){goe4=fabs(goe4);}else{goe4=-fabs(goe4);}
Xp=X+(1./6.)*(goe1+2*goe2+2*goe3+goe4);
Yp=Y+(1./6.)*(goe1*k1+2*goe2*k2+2*goe3*k3+goe4*k4);
alp=x1v(Xp,Yp,dec);

其中的goe1等等就是上述的h
請幫忙證明是否正確。謝謝！！