关于常微分方程组的求解，各位帮个忙吧
 

y''+0.0175*y'+39.47842*y=-0.1439
z''+2.1549*z'+(y''-9.8)*sin(z+2.274*pi)./0.06=[0.1392*cos(z+1.5*pi)-0.3809*sin(z+1.5*pi) ]./0.0033
用四阶龙格库塔法求解。其中y 和z是关于时间t的函数。在t=0时，且y’=y=z’=z=0。

我自己用的是ode45的命令先是求解了第一条方程式，但是发现这样的话无法满足第二条，主要是第二条公式中含有y'' 提示说是将二元微分方程组拆分成四个单独的一阶一维微分方程，然后再用龙格库塔同步求解，可是我在运行中总是出现error。

希望各位高手能够帮我解一下。很急


我的QQ 381559634