【求助】微分方程与普通方程组成的方程组求解

romanticer

以下方程系数可能有些复杂，请大牛们看完后给出代码，数字可换用简单的。谢谢
dSdt=-(7.37*10^7/(1+2.48*10^-3*A)+1.46*10^7/(1+2.48*10^-3*A)/H+3570/(1+2.48*10^-3*A)/H^2)*S*O-1.09*10^15/(1+13*H)*S*B

dOdt=-2.33*10^-5*(7.37*10^7/(1+2.48*10^-3*A)+1.46*10^7/(1+2.48*10^-3*A)/H+3570/(1+2.48*10^-3*A)/H^2)*S*O

dBdt=-1.08*10^15/(1+13*H)*S*B

A=1.24*(1+0.0129/H+7.76*10^-10/(H^2))
H 满足 H+2.17*10^10/(1+62(1+1.709*10^9*H))*H=10^-14/H+1.46*10^7/(1+2.48*10^-3*A)*S/H+7140/(1+2.48*10^-3*A)*S/H^2+2*(0.2-10^6*S)/(H^2/10.2+1+H/0.0102)+(0.2-10^6*S)/(H/1000+1+0.0102/H)

t=0时, H=10^-6.17, S=2*10^-7, O=4.7*10^-11, B=7.1*10^-5
求t从0到3600， H，S，O，B的数值解
求代码以供学习
总的来说，这是三个常微分方程和一个普通方程组成的方程组，知道初值，请问如何求数值解；或者说是变系数常微分方程组，其中系数被另一个复杂的分式方程限制。

另弱问这种情况ode45好使吗？还是我太弱。。。

WF1987

搜微分代数方程组   可解