gh2h2h2h3=-3*(L4*R21*R22^2*R31 - 2*cb*IL5*R21*R22*R23*R31 + 
  cb^2*L4*R21*R23^2*R31 + 
  L4*R21^2*R22*R32 + L2*R22^3*R32 - cb*IL5*R21^2*R23*R32 + 
  cb^2*L2*R22*R23^2*R32 - cb*IL5*R21^2*R22*R33 + cb^2*L4*R21^2*R23*R33 + 
  cb^2*L2*R22^2*R23*R33 + cb^4*L2*R23^3*R33 + R21*R22^2*R31*RL5 - 
  cb^2*R21*R23^2*R31*RL5 + R21^2*R22*R32*RL5 - cb^2*R21^2*R23*R33*RL5 - 
  IL5*R22^2*R23*R31*sb + cb^2*IL5*R23^3*R31*sb - 2*IL5*R21*R22*R23*R32*sb - 
  IL5*R21*R22^2*R33*sb + 3*cb^2*IL5*R21*R23^2*R33*sb - 
  2*cb*R22*R23^2*R31*RL5*sb - 2*cb*R21*R23^2*R32*RL5*sb - 
  4*cb*R21*R22*R23*R33*RL5*sb + L4*R22*R23^2*R32*sb^2 + 
  cb*IL5*R23^3*R32*sb^2 + L4*R22^2*R23*R33*sb^2 + 
  3*cb*IL5*R22*R23^2*R33*sb^2 + 2*cb^2*L4*R23^3*R33*sb^2 - 
  R22*R23^2*R32*RL5*sb^2 - R22^2*R23*R33*RL5*sb^2 + 
  2*cb^2*R23^3*R33*RL5*sb^2 + L1*(R21^2 + R23^2*sb^2)*
   (R21*R31 + R23*R33*sb^2) + L3*(R21*(R22^2 + cb^2*R23^2)*R31 + 
    R21^2*(R22*R32 + cb^2*R23*R33) + R23*(R22*R23*R32 + R22^2*R33 + 
      2*cb^2*R23^2*R33)*sb^2))
