runge-kutta法求解常微分方程组

x(t)'=(t-1)/0.001
I2(t)'=cos(x)
初值为：x(0)=0 I2(0)=1

这种套的常微分方程如何用runge-kutta法求解，直接套用runge-kutta法求解出来的图像是错的。各位程序界大神们帮帮忙。谢谢大家