m =

m

m

Beware: m^{jk} is not the “tensor of inertia”!
See below!
Sometimes m^{jk} is called the covariance matrix.

These moments are a slight generalization over the moments of Wikipedia on 2011 Jan 8.
Notational disaster has struck.
Heretofore context has made it clear whether a superscript was a power or a coordinate index.
The superscripts above are coordinate indexes, as in tensor notation, and my equations work in any number of dimensions but for moments of degree only up to two.
My notation can handle higher moments but each requires a new equation.
The equation for the third moments m^{ijk} is obvious.
This article uses almost exactly the same expression but where the exponents are powers.
Their equation is limited to 2D but describes higher moments.

See note on summation convention.
The coördinates of the center of mass (CM) are at m^{j}/m.

These 10 numbers characterize the rigid body inertially.
With no external forces, these 10 numbers together with initial position (3 numbers), velocity (3 numbers) and either angular velocity or angular momentum (3 numbers each), suffice to predict the future position, orientation, velocity and angular velocities.
The velocity of the CM and angular momentum are constant, of course.
The orientation can be given by an orthogonal matrix but really needs only 3 independent numbers.
The development below assumes that the point of the body at the origin is constrained not to move.
This matches the tumbling problem above if m^{j} = 0, a situation that may be brought about by translating the coördinate system.

In 3 dimensions there are 3 independent values in an antisymmetric tensor and there are 3 independent values in a vector.
Axial vectors are thus able to code the information that requires an antisymmetric tensor in higher dimensions.
The tensor of inertia, I_{ij} = δ_{ij}m^{kk} − m^{ij}, derived from the 2nd moments, fits well with these axial vectors for angular velocity, ω, and momentum, L, to make the convenient and simple formula “L = Iω”.
The tensor of inertia looses this magic beyond 3D however but the 2nd moments remain in control.

This code corroborates that L = Iω is correct in 3D when I is viewed as a symmetric matrix and L and ω are both (axial) vectors. Here is a generalized formula in n dimensions.

This program corroborates the conventional 3D relation L = Iω.

This differs from this only in style and notation.

The equations so far look very simple but the tensor of inertia must presumably be expressed in the laboratory frame and thus varies with time. If A

The time rate of change of A^{i}_{j} depends on the angular velocity ω_{ik} thus:
d/dt A^{i}_{j} = A^{i}_{k}ω_{jk}

A calculation cycle can proceed thus: We know L and I which are constant.

- We know the initial A.
- We calculate J from A and I
- We solve L = Jω for ω
- We compute a new A using ω (and δt)