gh3h3h3=-3*(L4*R31*R32^2*v1 - 2*cb*IL5*R31*R32*R33*v1 + cb^2*L4*R31*R33^2*v1 + 
  R31*R32^2*RL5*v1 - cb^2*R31*R33^2*RL5*v1 - IL5*R32^2*R33*sb*v1 + 
  cb^2*IL5*R33^3*sb*v1 - 2*cb*R32*R33^2*RL5*sb*v1 + 
  L1*R31*(R31^2 + R33^2*sb^2)*v1 + L4*R31^2*R32*v2 + L2*R32^3*v2 - 
  cb*IL5*R31^2*R33*v2 + cb^2*L2*R32*R33^2*v2 + R31^2*R32*RL5*v2 - 
  2*IL5*R31*R32*R33*sb*v2 - 2*cb*R31*R33^2*RL5*sb*v2 + L4*R32*R33^2*sb^2*v2 + 
  cb*IL5*R33^3*sb^2*v2 - R32*R33^2*RL5*sb^2*v2 + 
  L3*(R31*(R32^2 + cb^2*R33^2)*v1 + R31^2*R32*v2 + R32*R33^2*sb^2*v2))
