Solving indexed ODE in Mathematica


I would like to solve an ODE with indexed variables symbolically. Is there a way to do that? $ $ \frac{dp_{\beta\alpha}}{dt}=m\left(p_{1\alpha}p_{\beta2}+p_{2\alpha}p_{\beta1}\right)$ $ m is called mass and is constant. t is time.