def inertia_matrix(dyadic, rf):
    return Matrix([[dot(dot(dyadic, i), j) for j in rf] for i in rf])
Return the inertia matrix of a given dyadic for a specified
dataset/ETHPy150Open pydy/pydy/examples/Kane1985/Chapter3/Ex6.14.py/inertia_matrix