        gh1HpHm=-(cb**2*L3*v1*Rmat(1,1)) - L1*sb**2*v1*Rmat(1,1) + 
     &  cb*L4*sb*v2*Rmat(1,1) + cb*RL5*sb*v2*Rmat(1,1) + 
     &  cb*L4*sb*v1*Rmat(1,2) + cb*RL5*sb*v1*Rmat(1,2) - 
     &  cb**2*L2*v2*Rmat(1,2) - L3*sb**2*v2*Rmat(1,2) - 
     &  cb**2*IL5*sb*v1*Rmat(1,3) - 
     &  cb*IL5*sb**2*v2*Rmat(1,3)
