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