code: def test_apply_loads_on_multi_degree_freedom_holonomic_system():
      """Example based on: https://pydy.readthedocs.io/en/latest/examples/multidof-holonomic.html"""
      W = Body("W")  # Wall
      B = Body("B")  # Block
      P = Body("P")  # Pendulum
      b = Body("b")  # bob
      q1, q2 = dynamicsymbols("q1 q2")  # generalized coordinates
      k, c, g, kT = symbols("k c g kT")  # constants
      F, T = dynamicsymbols("F T")  # Specified forces
  
      # Applying forces
      B.apply_force(F * W.x)
      W.apply_force(k * q1 * W.x, reaction_body=B)  # Spring force
      W.apply_force(c * q1.diff() * W.x, reaction_body=B)  # dampner
      P.apply_force(P.mass * g * W.y)
      b.apply_force(b.mass * g * W.y)
  
      # Applying torques
      P.apply_torque(kT * q2 * W.z, reaction_body=b)
      P.apply_torque(T * W.z)
  
      assert B.loads == [
          (B.masscenter, (F - k * q1 - c * q1.diff()) * W.x)
      ]
      assert P.loads == [
          (P.masscenter, P.mass * g * W.y),
          (P.frame, (T + kT * q2) * W.z),
      ]
      assert b.loads == [
          (b.masscenter, b.mass * g * W.y),
          (b.frame, -kT * q2 * W.z),
      ]
      assert W.loads == [(W.masscenter, (c * q1.diff() + k * q1) * W.x)]
  
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 16 for calling a global function with example usage: test_apply_loads_on_multi_degree_freedom_holonomic_system() without return types