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