gh2h2h3=-(L4*R22^2*R31*v1) + 2*cb*IL5*R22*R23*R31*v1 - cb^2*L4*R23^2*R31*v1 - 
 2*L4*R21*R22*R32*v1 + 2*cb*IL5*R21*R23*R32*v1 + 2*cb*IL5*R21*R22*R33*v1 - 
 2*cb^2*L4*R21*R23*R33*v1 - R22^2*R31*RL5*v1 + cb^2*R23^2*R31*RL5*v1 - 
 2*R21*R22*R32*RL5*v1 + 2*cb^2*R21*R23*R33*RL5*v1 + 2*IL5*R22*R23*R32*sb*v1 + 
 IL5*R22^2*R33*sb*v1 - 3*cb^2*IL5*R23^2*R33*sb*v1 + 
 2*cb*R23^2*R32*RL5*sb*v1 + 4*cb*R22*R23*R33*RL5*sb*v1 - 
 L1*(3*R21^2*R31 + R23^2*R31*sb^2 + 2*R21*R23*R33*sb^2)*v1 - 
 2*L4*R21*R22*R31*v2 + 2*cb*IL5*R21*R23*R31*v2 - L4*R21^2*R32*v2 - 
 3*L2*R22^2*R32*v2 - cb^2*L2*R23^2*R32*v2 + cb*IL5*R21^2*R33*v2 - 
 2*cb^2*L2*R22*R23*R33*v2 - 2*R21*R22*R31*RL5*v2 - R21^2*R32*RL5*v2 + 
 2*IL5*R22*R23*R31*sb*v2 + 2*IL5*R21*R23*R32*sb*v2 + 
 2*IL5*R21*R22*R33*sb*v2 + 2*cb*R23^2*R31*RL5*sb*v2 + 
 4*cb*R21*R23*R33*RL5*sb*v2 - L4*R23^2*R32*sb^2*v2 - 
 2*L4*R22*R23*R33*sb^2*v2 - 3*cb*IL5*R23^2*R33*sb^2*v2 + 
 R23^2*R32*RL5*sb^2*v2 + 2*R22*R23*R33*RL5*sb^2*v2 - 
 L3*(R22^2*R31*v1 + cb^2*R23*(R23*R31 + 2*R21*R33)*v1 + 
   R32*(R21^2 + R23^2*sb^2)*v2 + 2*R22*(R21*R32*v1 + R21*R31*v2 + 
     R23*R33*sb^2*v2))
