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