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