text: import sympy.physics.mechanics as me
  import sympy as sm
  import math as m
  import numpy as np
  
  f = sm.S(3)
  g = sm.S(9.81)
  a, b = sm.symbols('a b', real=True)
  s, s1 = sm.symbols('s s1', real=True)
  s2, s3 = sm.symbols('s2 s3', real=True, nonnegative=True)
  s4 = sm.symbols('s4', real=True, nonpositive=True)
  k1, k2, k3, k4, l1, l2, l3, p11, p12, p13, p21, p22, p23 = sm.symbols('k1 k2 k3 k4 l1 l2 l3 p11 p12 p13 p21 p22 p23', real=True)
  c11, c12, c13, c21, c22, c23 = sm.symbols('c11 c12 c13 c21 c22 c23', real=True)
  e1 = a*f+s2-g
  e2 = f**2+k3*k2*g
  