code: def test_remove_load():
      E = Symbol("E")
      I = Symbol("I")
      b = Beam(4, E, I)
  
      try:
          b.remove_load(2, 1, -1)
      # As no load is applied on beam, ValueError should be returned.
      except ValueError:
          assert True
      else:
          assert False
  
      b.apply_load(-3, 0, -2)
      b.apply_load(4, 2, -1)
      b.apply_load(-2, 2, 2, end=3)
      b.remove_load(-2, 2, 2, end=3)
      assert b.load == -3 * SingularityFunction(
          x, 0, -2
      ) + 4 * SingularityFunction(x, 2, -1)
      assert b.applied_loads == [(-3, 0, -2, None), (4, 2, -1, None)]
  
      try:
          b.remove_load(1, 2, -1)
      # As load of this magnitude was never applied at
      # this position, method should return a ValueError.
      except ValueError:
          assert True
      else:
          assert False
  
      b.remove_load(-3, 0, -2)
      b.remove_load(4, 2, -1)
      assert b.load == 0
      assert b.applied_loads == []
  
imports: from sympy.core.function import expand
  from sympy.core.numbers import Rational, pi
  from sympy.core.singleton import S
  from sympy.core.symbol import Symbol, symbols
  from sympy.sets.sets import Interval
  from sympy.simplify.simplify import simplify
  from sympy.physics.continuum_mechanics.beam import Beam
  from sympy.functions import SingularityFunction, Piecewise, meijerg, Abs, log
  from sympy.testing.pytest import raises
  from sympy.physics.units import meter, newton, kilo, giga, milli
  from sympy.physics.continuum_mechanics.beam import Beam3D
  from sympy.geometry import Circle, Polygon, Point2D, Triangle
  from sympy.core.sympify import sympify
label_desc: Use image node_id 8 for calling a global function with example usage: test_remove_load() without return types