Solution (source code)

= Solution

With total mass $M=\sum_i m_i$, define
$$
R=\frac1M\sum_i m_ix_i,
\qquad
P=\sum_i m_i\dot x_i,
\qquad
L=\sum_i x_i\times m_i\dot x_i.
$$
The <angular momentum> about the centre of mass is
$$
L_{\rm CoM}=\sum_i(x_i-R)\times m_i(\dot x_i-\dot R).
$$
Expanding, using $\sum_i m_i(x_i-R)=0$ and $P=M\dot R$, gives
$$
\boxed{L_{\rm CoM}=L-R\times P}.
$$