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