[MATLAB数学相关] MATLAB求解微分方程组

[lxh.3888]

哪位高手能帮忙解一下,用matlab龙格库塔求解下列微分方程组,万分感谢!
未知量：Tu、Tb、Mb
其他的参数先假设已知，在后面可以采用子程序调用；
dTu/dt=1/((m-Mb)*Cpu)*(Vu*dP+dQu);
dTb/dt=1/(Rb*Mb)*(V*dP-(Rb*Tb-Ru*Tu)*dMb-Ru/Cpu*(Vu*dP+dQu));
-dP*(Vu*(Cvu/Cpu-(Ru/Rb)*(Cvb/Cpu)))=((330.22-Rb*Tb-1716+Ru*Tu)-Cvb*(Tu-(Ru/Rb)*Tb)*dMb/dt+(Cvu/Cpu-(Cvb/Cpu)*(Ru/Rb))*dQu-(dQu+dQb);

谢谢，急用！