        gh1h1h2h2=-3*L1*Rmat(1,1)**2*Rmat(2,1)**2 - 
     &  L3*Rmat(1,2)**2*Rmat(2,1)**2 - 
     &  L4*Rmat(1,2)**2*Rmat(2,1)**2 - 
     &  RL5*Rmat(1,2)**2*Rmat(2,1)**2 + 
     &  2*cb*IL5*Rmat(1,2)*Rmat(1,3)*Rmat(2,1)**2 - 
     &  cb**2*L3*Rmat(1,3)**2*Rmat(2,1)**2 - 
     &  cb**2*L4*Rmat(1,3)**2*Rmat(2,1)**2 + 
     &  cb**2*RL5*Rmat(1,3)**2*Rmat(2,1)**2 - 
     &  L1*sb**2*Rmat(1,3)**2*Rmat(2,1)**2 - 
     &  4*L3*Rmat(1,1)*Rmat(1,2)*Rmat(2,1)*Rmat(2,2) - 
     &  4*L4*Rmat(1,1)*Rmat(1,2)*Rmat(2,1)*Rmat(2,2) - 
     &  4*RL5*Rmat(1,1)*Rmat(1,2)*Rmat(2,1)*Rmat(2,2) + 
     &  4*cb*IL5*Rmat(1,1)*Rmat(1,3)*Rmat(2,1)*Rmat(2,2) + 
     &  4*IL5*sb*Rmat(1,2)*Rmat(1,3)*Rmat(2,1)*Rmat(2,2) + 
     &  4*cb*RL5*sb*Rmat(1,3)**2*Rmat(2,1)*Rmat(2,2) - 
     &  L3*Rmat(1,1)**2*Rmat(2,2)**2 - 
     &  L4*Rmat(1,1)**2*Rmat(2,2)**2 - 
     &  RL5*Rmat(1,1)**2*Rmat(2,2)**2 - 
     &  3*L2*Rmat(1,2)**2*Rmat(2,2)**2 + 
     &  2*IL5*sb*Rmat(1,1)*Rmat(1,3)*Rmat(2,2)**2 - 
     &  cb**2*L2*Rmat(1,3)**2*Rmat(2,2)**2 - 
     &  L3*sb**2*Rmat(1,3)**2*Rmat(2,2)**2 - 
     &  L4*sb**2*Rmat(1,3)**2*Rmat(2,2)**2 + 
     &  RL5*sb**2*Rmat(1,3)**2*Rmat(2,2)**2 + 
     &  4*cb*IL5*Rmat(1,1)*Rmat(1,2)*Rmat(2,1)*Rmat(2,3) + 
     &  2*IL5*sb*Rmat(1,2)**2*Rmat(2,1)*Rmat(2,3) - 
     &  4*cb**2*L3*Rmat(1,1)*Rmat(1,3)*Rmat(2,1)*
     &   Rmat(2,3) - 4*cb**2*L4*Rmat(1,1)*Rmat(1,3)*
     &   Rmat(2,1)*Rmat(2,3) + 
     &  4*cb**2*RL5*Rmat(1,1)*Rmat(1,3)*Rmat(2,1)*
     &   Rmat(2,3) - 4*L1*sb**2*Rmat(1,1)*Rmat(1,3)*
     &   Rmat(2,1)*Rmat(2,3) + 
     &  8*cb*RL5*sb*Rmat(1,2)*Rmat(1,3)*Rmat(2,1)*
     &   Rmat(2,3) - 6*cb**2*IL5*sb*Rmat(1,3)**2*Rmat(2,1)*
     &   Rmat(2,3) + 2*cb*IL5*Rmat(1,1)**2*Rmat(2,2)*
     &   Rmat(2,3) + 4*IL5*sb*Rmat(1,1)*Rmat(1,2)*Rmat(2,2)*
     &   Rmat(2,3) + 8*cb*RL5*sb*Rmat(1,1)*Rmat(1,3)*
     &   Rmat(2,2)*Rmat(2,3) - 
     &  4*cb**2*L2*Rmat(1,2)*Rmat(1,3)*Rmat(2,2)*
     &   Rmat(2,3) - 4*L3*sb**2*Rmat(1,2)*Rmat(1,3)*
     &   Rmat(2,2)*Rmat(2,3) - 
     &  4*L4*sb**2*Rmat(1,2)*Rmat(1,3)*Rmat(2,2)*
     &   Rmat(2,3) + 4*RL5*sb**2*Rmat(1,2)*Rmat(1,3)*
     &   Rmat(2,2)*Rmat(2,3) - 
     &  6*cb*IL5*sb**2*Rmat(1,3)**2*Rmat(2,2)*Rmat(2,3) - 
     &  cb**2*L3*Rmat(1,1)**2*Rmat(2,3)**2 - 
     &  cb**2*L4*Rmat(1,1)**2*Rmat(2,3)**2 + 
     &  cb**2*RL5*Rmat(1,1)**2*Rmat(2,3)**2 - 
     &  L1*sb**2*Rmat(1,1)**2*Rmat(2,3)**2 + 
     &  4*cb*RL5*sb*Rmat(1,1)*Rmat(1,2)*Rmat(2,3)**2 - 
     &  cb**2*L2*Rmat(1,2)**2*Rmat(2,3)**2 - 
     &  L3*sb**2*Rmat(1,2)**2*Rmat(2,3)**2 - 
     &  L4*sb**2*Rmat(1,2)**2*Rmat(2,3)**2 + 
     &  RL5*sb**2*Rmat(1,2)**2*Rmat(2,3)**2 - 
     &  6*cb**2*IL5*sb*Rmat(1,1)*Rmat(1,3)*Rmat(2,3)**2 - 
     &  6*cb*IL5*sb**2*Rmat(1,2)*Rmat(1,3)*Rmat(2,3)**2 - 
     &  3*cb**4*L2*Rmat(1,3)**2*Rmat(2,3)**2 - 
     &  6*cb**2*L3*sb**2*Rmat(1,3)**2*Rmat(2,3)**2 - 
     &  6*cb**2*L4*sb**2*Rmat(1,3)**2*Rmat(2,3)**2 - 
     &  6*cb**2*RL5*sb**2*Rmat(1,3)**2*Rmat(2,3)**2 - 
     &  3*L1*sb**4*Rmat(1,3)**2*Rmat(2,3)**2
