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