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