syms P(t); P(t) = dsolve(diff(P,t) == 1.5-1.5*P, P(0)==2); P(t) = simplify(P)