code: def test_remove_load():
      P1 = Point("P1")
      P2 = Point("P2")
      B = Body("B")
      f1 = B.x
      f2 = B.y
      B.apply_force(f1, P1)
      B.apply_force(f2, P2)
      assert B.loads == [(P1, f1), (P2, f2)]
      B.remove_load(P2)
      assert B.loads == [(P1, f1)]
      B.apply_torque(f1.cross(f2))
      assert B.loads == [(P1, f1), (B.frame, f1.cross(f2))]
      B.remove_load()
      assert B.loads == [(P1, f1)]
  
imports: from sympy import Symbol, symbols, sin, cos, Matrix, zeros, simplify
  from sympy.physics.vector import Point, ReferenceFrame, dynamicsymbols, Dyadic
  from sympy.physics.mechanics import inertia, Body
  from sympy.testing.pytest import raises
label_desc: Use image node_id 15 for calling a global function with example usage: test_remove_load() without return types