MATLAB: Solve system of diferential equations. how matlab can use 4) and 5) to solve the system

1) d(na(t))/dt=-ko*exp(-ER/T(t))*na
2) d(nw(t))/dt=-ko*exp(-ER/T(t))*nw
3) d(nb(t))/dt=-ko*exp(-ER/T(t))*nb
4) d(h(t))/dt=15*10^3*(6*T(t)-360)
5) h(t)=nb(t)*T(t)*250-na(t)*T(t)*125-nw(t)*T(t)*75+95945*na(t)+242875*nw(t)-430450*nb(t)

Best Answer

Solve 5) for T and insert the expression in 1) - 4).
Then use ODE15S to solve 1) - 4).
Best wishes
Torsten.