Well there is the program but first some comments.
This Algol 68 program verifies numerically a few of the derived equations but also illustrates the failure of one simple attempt at an n-dimensional extension.
Here is a guide to the program for those who don’t know Algol 68.
The Cast
ma is an array of coordinates for 20 unit mass points.
omev is the axial vector describing the angular velocity.
ccl is the angular momentum of the collective mass points.
m2 is the second moments of the set of masses.
ti is the conventional tensor of inertial computed from the second moments.
cxl is the angular momentum computed from the tensor of inertial and angular velocity as L = Iω.
Good News
These two computations of angular momentum L agree numerically.
Bad News
The obvious ways to compute L from the conventional tensor of inertia while avoiding axial vectors, all fail.
Axial vectors do not generalize beyond 3D except as anti-symmetric tensors.
Note the rank 4 tensor it2 at end that is computed from 2nd moments which gets the right answer.