gh3h3h3h3=-3*(2*L4*R31^2*R32^2 + L2*R32^4 - 4*cb*IL5*R31^2*R32*R33 + 
  2*cb^2*L4*R31^2*R33^2 + 2*cb^2*L2*R32^2*R33^2 + cb^4*L2*R33^4 + 
  2*R31^2*R32^2*RL5 - 2*cb^2*R31^2*R33^2*RL5 - 4*IL5*R31*R32^2*R33*sb + 
  4*cb^2*IL5*R31*R33^3*sb - 8*cb*R31*R32*R33^2*RL5*sb + 
  2*L4*R32^2*R33^2*sb^2 + 4*cb*IL5*R32*R33^3*sb^2 + 2*cb^2*L4*R33^4*sb^2 - 
  2*R32^2*R33^2*RL5*sb^2 + 2*cb^2*R33^4*RL5*sb^2 + 
  2*L3*(R32^2 + cb^2*R33^2)*(R31^2 + R33^2*sb^2) + L1*(R31^2 + R33^2*sb^2)^2)
