gh1h1h1h1=-3*(2*L4*R11^2*R12^2 + L2*R12^4 - 4*cb*IL5*R11^2*R12*R13 + 
  2*cb^2*L4*R11^2*R13^2 + 2*cb^2*L2*R12^2*R13^2 + cb^4*L2*R13^4 + 
  2*R11^2*R12^2*RL5 - 2*cb^2*R11^2*R13^2*RL5 - 4*IL5*R11*R12^2*R13*sb + 
  4*cb^2*IL5*R11*R13^3*sb - 8*cb*R11*R12*R13^2*RL5*sb + 
  2*L4*R12^2*R13^2*sb^2 + 4*cb*IL5*R12*R13^3*sb^2 + 2*cb^2*L4*R13^4*sb^2 - 
  2*R12^2*R13^2*RL5*sb^2 + 2*cb^2*R13^4*RL5*sb^2 + 
  2*L3*(R12^2 + cb^2*R13^2)*(R11^2 + R13^2*sb^2) + L1*(R11^2 + R13^2*sb^2)^2)
