Ornstein-Uhlenbeck moment equations (source code)

= Ornstein-Uhlenbeck moment equations
{c}
{title2=$M'=0,\quad m'=-m,\quad E'=dM-2E$}

<Integration by parts> against $1$, $v_i$ and $|v|^2/2$ gives these equations for <mass>, unnormalized <momentum> and <kinetic energy>. With $M_0\ne0$, the mean is $u(t)=e^{-t}u_0$, and $E(t)=dM_0/2+(E_0-dM_0/2)e^{-2t}$. Normalizing by <mass> is invalid if that <mass> is zero.