code: def test_inertia_of_point_mass():
      r, s, t, m = symbols("r s t m")
      N = ReferenceFrame("N")
  
      px = r * N.x
      I = inertia_of_point_mass(m, px, N)
      assert I == m * r**2 * (N.y | N.y) + m * r**2 * (N.z | N.z)
  
      py = s * N.y
      I = inertia_of_point_mass(m, py, N)
      assert I == m * s**2 * (N.x | N.x) + m * s**2 * (N.z | N.z)
  
      pz = t * N.z
      I = inertia_of_point_mass(m, pz, N)
      assert I == m * t**2 * (N.x | N.x) + m * t**2 * (N.y | N.y)
  
      p = px + py + pz
      I = inertia_of_point_mass(m, p, N)
      assert I == (
          m * (s**2 + t**2) * (N.x | N.x)
          - m * r * s * (N.x | N.y)
          - m * r * t * (N.x | N.z)
          - m * r * s * (N.y | N.x)
          + m * (r**2 + t**2) * (N.y | N.y)
          - m * s * t * (N.y | N.z)
          - m * r * t * (N.z | N.x)
          - m * s * t * (N.z | N.y)
          + m * (r**2 + s**2) * (N.z | N.z)
      )
  
imports: from sympy import symbols
  from sympy.testing.pytest import raises
  from sympy.physics.mechanics import inertia, inertia_of_point_mass, Inertia, ReferenceFrame, Point
label_desc: Use image node_id 2 for calling a global function with example usage: test_inertia_of_point_mass() without return types