龙格库塔法的算法疑问，求大神解答！

调用的时候这样，[t,x,y,vx,vy]=runge_kuta(@fun_lx,@fun_ly,[0 2],0.01,0,240.8,10688,0);
fun_lx,fun_ly另外写，图像结果似乎没问题，但当我把时间取0到3时，结果就变得奇怪了，似乎哪里不对

是不是我的龙格库塔算法里的x 套错了？求指导

matlab自带龙格库塔函数，ode23和ode45都是