Solution
= Solution
Eliminating $P_2=1-P_1$ gives $\dot P_1=\beta-(\alpha+\beta)P_1$. Therefore
$$
P_1(t)=\frac\beta{\alpha+\beta}
+\left(p-\frac\beta{\alpha+\beta}\right)e^{-(\alpha+\beta)t}.
$$
With $r=\alpha p-\beta(1-p)$,
$$
\boxed{P_1=\frac{\beta+re^{-(\alpha+\beta)t}}{\alpha+\beta},
\qquad P_2=\frac{\alpha-re^{-(\alpha+\beta)t}}{\alpha+\beta}.}
$$