Newton's Law of Cooling Model

> IVP := {Diff(T(t),t)=alpha*(M-T(t)), T(t_0)=T_0};
Newton := collect(simplify(rhs(dsolve(IVP, T(t)))), exp);

IVP := {T(t_0) = T_0, Diff(T(t),t) = alpha*(M-T(t))...

Newton := M+(-M+T_0)*exp(-alpha*(t-t_0))