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