【SOS】在MATLAB中用4阶RK方法求解此一阶常微分方程组
 

我是新手谁能帮我看看这个微分方程组怎么解？谢谢了

dy/dt=-acosx-bcos[kx0-(k-1)x]
dx/dt={asinx+vsin[kx0-(k-1)x]}/y

其中x0是x下脚标0，就是说x的初始值。。。
a b都是已知的。。。

x0=m y0=n 就是随便已知一个数字就行，哪位大侠能帮忙实现啊？

仿真步长为0.1 时间上下限为0 100