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