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