gh2h2h2=-3*(L4*R21*R22^2*v1 - 2*cb*IL5*R21*R22*R23*v1 + cb^2*L4*R21*R23^2*v1 + 
  R21*R22^2*RL5*v1 - cb^2*R21*R23^2*RL5*v1 - IL5*R22^2*R23*sb*v1 + 
  cb^2*IL5*R23^3*sb*v1 - 2*cb*R22*R23^2*RL5*sb*v1 + 
  L1*R21*(R21^2 + R23^2*sb^2)*v1 + L4*R21^2*R22*v2 + L2*R22^3*v2 - 
  cb*IL5*R21^2*R23*v2 + cb^2*L2*R22*R23^2*v2 + R21^2*R22*RL5*v2 - 
  2*IL5*R21*R22*R23*sb*v2 - 2*cb*R21*R23^2*RL5*sb*v2 + L4*R22*R23^2*sb^2*v2 + 
  cb*IL5*R23^3*sb^2*v2 - R22*R23^2*RL5*sb^2*v2 + 
  L3*(R21*(R22^2 + cb^2*R23^2)*v1 + R21^2*R22*v2 + R22*R23^2*sb^2*v2))
