求助

[caihaoren]

求解方程组：
(m1c1+m2c2+m2c2x1)*dTch/dt=W1 ；(里面是Tch关于时间t的求导）
Tout=Tch+(Tin -Tch)*e^(-aA/mwcw)
W1=cwmw(Tin-Tout)
已知量：m1=273.2；c1=0.469；m2=12.9；c2=2；Tin=80；cw=4.2；mw=0.83；
a=120,A=1.86
x1=0.000004Tch4+0.0014Tch3-0.1573Tch2+7.9031Tch-145.7
求Tch与时间t的关系（用matlab求解出图像）
哪位牛人能帮我把主要程序写出来就ok了，谢谢大家！