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);