gh2h2h2h2=-3*(2*L4*R21^2*R22^2 + L2*R22^4 - 4*cb*IL5*R21^2*R22*R23 + 
  2*cb^2*L4*R21^2*R23^2 + 2*cb^2*L2*R22^2*R23^2 + cb^4*L2*R23^4 + 
  2*R21^2*R22^2*RL5 - 2*cb^2*R21^2*R23^2*RL5 - 4*IL5*R21*R22^2*R23*sb + 
  4*cb^2*IL5*R21*R23^3*sb - 8*cb*R21*R22*R23^2*RL5*sb + 
  2*L4*R22^2*R23^2*sb^2 + 4*cb*IL5*R22*R23^3*sb^2 + 2*cb^2*L4*R23^4*sb^2 - 
  2*R22^2*R23^2*RL5*sb^2 + 2*cb^2*R23^4*RL5*sb^2 + 
  2*L3*(R22^2 + cb^2*R23^2)*(R21^2 + R23^2*sb^2) + L1*(R21^2 + R23^2*sb^2)^2)
