code: def P_c(x, L=L, E=E, G=G):
      return (
          (4.013 * E / (6 * L**2))
          * x[2]
          * x[3] ** 3
          * (1 - 0.25 * x[2] * (E / G) ** 0.5 / L)
      )
  
file_path: mystic/examples2/beam.py
imports: from mystic.penalty import quadratic_inequality
label_desc: Use image node_id 5 for calling a global function with example usage: P_c(x, L, E, G) and returns: unknown