def get_inertia_tensor(elements, coordinates):
pow2 = coordinates**<NUM_LIT:2><EOL>molecular_weight = np.array(<EOL>[[atomic_mass[e.upper()]] for e in elements])<EOL>diag_1 = np.sum(molecular_weight * (pow2[:, <NUM_LIT:1>] + pow2[:, <NUM_LIT:2>]))<EOL>diag_2 = np.sum(molecular_weight * (pow2[:, <NUM_LIT:0>] + pow2[:, <NUM_LIT:2>]))<EOL>diag_3 = np.sum(molecular_weight * (pow2[:, <NUM_LIT:0>] + pow2[:, <NUM_LIT:1>]))<EOL>mxy = np.sum(-molecular_weight * coordinates[:, <NUM_LIT:0>] * coordinates[:, <NUM_LIT:1>])<EOL>mxz = np.sum(-molecular_weight * coordinates[:, <NUM_LIT:0>] * coordinates[:, <NUM_LIT:2>])<EOL>myz = np.sum(-molecular_weight * coordinates[:, <NUM_LIT:1>] * coordinates[:, <NUM_LIT:2>])<EOL>inertia_tensor = np.array([[diag_1, mxy, mxz], [mxy, diag_2, myz],<EOL>[mxz, myz, diag_3]]) / coordinates.shape[<NUM_LIT:0>]<EOL>return (inertia_tensor)<EOL>
Return the tensor of inertia a molecule.

Parameters
----------
elements : numpy.ndarray
    The array containing the molecule's elemental data.

coordinates :  numpy.ndarray
    The array containing the Cartesian coordinates of the molecule.

Returns
-------
numpy.ndarray
    The tensor of inertia of a molecule.