We can use deSolve package in R for ordinary differential equations (ODE), however, I can't find a way to solve two nested ODE equation, suppose b(t)' = beta - k*b(t); a(t)' = alpha -b(t)*gamma; where ' means differentiation. How can we solve a and b then? as a' is a function of b, we have to solve a and b simulataneously. package in R for ordinary di