Use differential equations in models (#580)

Use differential equations in models (#580)

Solve dydt=ky\frac{dy}{dt} = ky to get y=y0ekty = y_0 e^{kt}.