FORM by J.Vermaseren,version 3.1(Sep  4 2005) Run at: Thu May 17 08:04:27 2007
    #+
    
    * Lagrangian, Eq.18:
    *  Quantization and renormalization of the manifest left-right symmetric 
    *   model of electroweak interactions.
    *    P. Duka, J. Gluza, M. Zralek, Annals Phys.280:336-408,2000.
    *         e-Print: hep-ph/9910279
    
    * daje macierz masowa (w polach niefizycznych) 
    #define masy 
    
    off statistics;
    s g, gp;
    s i, j, k;
    cf DLDeltaLh, DLDeltaL, DRDeltaRh, DRDeltaR, DPhih, DPhi; 
    cf ad,d;                                                                    
    cf B, WLh, WL, WRh, WR, DeltaLh, DeltaL, DeltaRh, DeltaR, Phih, Phi;
    
    ************************************************************
    
    g exp = sum_(i, 1, 2, sum_(j, 1, 2, DLDeltaLh(i, j)*DLDeltaL(j, i)
                                        +DRDeltaRh(i, j)*DRDeltaR(j, i)
                                        +DPhih(i, j)*DPhi(j, i)));
    
    id DLDeltaLh(i?, j?) = d(DeltaLh(i, j))
                             -i_*g*sum_(k, 1, 2, WLh(i, k)*DeltaLh(k, j)
                             -DeltaLh(i, k)*WLh(k, j))
                             +i_*gp*B*DeltaLh(i, j);
    id DLDeltaL(i?, j?) = d(DeltaL(i, j))
                           -i_*g*sum_(k, 1, 2, WL(i, k)*DeltaL(k, j)
                           -DeltaL(i, k)*WL(k, j))
                           -i_*gp*B*DeltaL(i, j);
    id DRDeltaRh(i?, j?) = d(DeltaRh(i, j))
                             -i_*g*sum_(k, 1, 2, WRh(i, k)*DeltaRh(k, j)
                             -DeltaRh(i, k)*WRh(k, j))
                             +i_*gp*B*DeltaRh(i, j);
    id DRDeltaR(i?, j?) = d(DeltaR(i, j))
                           -i_*g*sum_(k, 1, 2, WR(i, k)*DeltaR(k, j)
                           -DeltaR(i, k)*WR(k, j))
                           -i_*gp*B*DeltaR(i, j);
    id DPhih(i?, j?) = d(Phih(i, j))
                         -i_*g*sum_(k, 1, 2, WRh(i, k)*Phih(k, j)
                         -Phih(i, k)*WLh(k, j));
    id DPhi(i?, j?) = d(Phi(i, j))
                       -i_*g*sum_(k, 1, 2, WL(i, k)*Phi(k, j)-Phi(i, k)*WR(k, j));
    
    
    print +s;
    .sort

   exp =
       + d(DeltaLh(1,1))*d(DeltaL(1,1))
       - d(DeltaLh(1,1))*B*DeltaL(1,1)*i_*gp
       - d(DeltaLh(1,1))*WL(1,2)*DeltaL(2,1)*i_*g
       + d(DeltaLh(1,1))*WL(2,1)*DeltaL(1,2)*i_*g
       + d(DeltaLh(1,2))*d(DeltaL(2,1))
       - d(DeltaLh(1,2))*B*DeltaL(2,1)*i_*gp
       + d(DeltaLh(1,2))*WL(1,1)*DeltaL(2,1)*i_*g
       - d(DeltaLh(1,2))*WL(2,1)*DeltaL(1,1)*i_*g
       + d(DeltaLh(1,2))*WL(2,1)*DeltaL(2,2)*i_*g
       - d(DeltaLh(1,2))*WL(2,2)*DeltaL(2,1)*i_*g
       + d(DeltaLh(2,1))*d(DeltaL(1,2))
       - d(DeltaLh(2,1))*B*DeltaL(1,2)*i_*gp
       - d(DeltaLh(2,1))*WL(1,1)*DeltaL(1,2)*i_*g
       + d(DeltaLh(2,1))*WL(1,2)*DeltaL(1,1)*i_*g
       - d(DeltaLh(2,1))*WL(1,2)*DeltaL(2,2)*i_*g
       + d(DeltaLh(2,1))*WL(2,2)*DeltaL(1,2)*i_*g
       + d(DeltaLh(2,2))*d(DeltaL(2,2))
       - d(DeltaLh(2,2))*B*DeltaL(2,2)*i_*gp
       + d(DeltaLh(2,2))*WL(1,2)*DeltaL(2,1)*i_*g
       - d(DeltaLh(2,2))*WL(2,1)*DeltaL(1,2)*i_*g
       + d(DeltaL(1,1))*B*DeltaLh(1,1)*i_*gp
       - d(DeltaL(1,1))*WLh(1,2)*DeltaLh(2,1)*i_*g
       + d(DeltaL(1,1))*WLh(2,1)*DeltaLh(1,2)*i_*g
       + d(DeltaL(1,2))*B*DeltaLh(2,1)*i_*gp
       + d(DeltaL(1,2))*WLh(1,1)*DeltaLh(2,1)*i_*g
       - d(DeltaL(1,2))*WLh(2,1)*DeltaLh(1,1)*i_*g
       + d(DeltaL(1,2))*WLh(2,1)*DeltaLh(2,2)*i_*g
       - d(DeltaL(1,2))*WLh(2,2)*DeltaLh(2,1)*i_*g
       + d(DeltaL(2,1))*B*DeltaLh(1,2)*i_*gp
       - d(DeltaL(2,1))*WLh(1,1)*DeltaLh(1,2)*i_*g
       + d(DeltaL(2,1))*WLh(1,2)*DeltaLh(1,1)*i_*g
       - d(DeltaL(2,1))*WLh(1,2)*DeltaLh(2,2)*i_*g
       + d(DeltaL(2,1))*WLh(2,2)*DeltaLh(1,2)*i_*g
       + d(DeltaL(2,2))*B*DeltaLh(2,2)*i_*gp
       + d(DeltaL(2,2))*WLh(1,2)*DeltaLh(2,1)*i_*g
       - d(DeltaL(2,2))*WLh(2,1)*DeltaLh(1,2)*i_*g
       + d(DeltaRh(1,1))*d(DeltaR(1,1))
       - d(DeltaRh(1,1))*B*DeltaR(1,1)*i_*gp
       - d(DeltaRh(1,1))*WR(1,2)*DeltaR(2,1)*i_*g
       + d(DeltaRh(1,1))*WR(2,1)*DeltaR(1,2)*i_*g
       + d(DeltaRh(1,2))*d(DeltaR(2,1))
       - d(DeltaRh(1,2))*B*DeltaR(2,1)*i_*gp
       + d(DeltaRh(1,2))*WR(1,1)*DeltaR(2,1)*i_*g
       - d(DeltaRh(1,2))*WR(2,1)*DeltaR(1,1)*i_*g
       + d(DeltaRh(1,2))*WR(2,1)*DeltaR(2,2)*i_*g
       - d(DeltaRh(1,2))*WR(2,2)*DeltaR(2,1)*i_*g
       + d(DeltaRh(2,1))*d(DeltaR(1,2))
       - d(DeltaRh(2,1))*B*DeltaR(1,2)*i_*gp
       - d(DeltaRh(2,1))*WR(1,1)*DeltaR(1,2)*i_*g
       + d(DeltaRh(2,1))*WR(1,2)*DeltaR(1,1)*i_*g
       - d(DeltaRh(2,1))*WR(1,2)*DeltaR(2,2)*i_*g
       + d(DeltaRh(2,1))*WR(2,2)*DeltaR(1,2)*i_*g
       + d(DeltaRh(2,2))*d(DeltaR(2,2))
       - d(DeltaRh(2,2))*B*DeltaR(2,2)*i_*gp
       + d(DeltaRh(2,2))*WR(1,2)*DeltaR(2,1)*i_*g
       - d(DeltaRh(2,2))*WR(2,1)*DeltaR(1,2)*i_*g
       + d(DeltaR(1,1))*B*DeltaRh(1,1)*i_*gp
       - d(DeltaR(1,1))*WRh(1,2)*DeltaRh(2,1)*i_*g
       + d(DeltaR(1,1))*WRh(2,1)*DeltaRh(1,2)*i_*g
       + d(DeltaR(1,2))*B*DeltaRh(2,1)*i_*gp
       + d(DeltaR(1,2))*WRh(1,1)*DeltaRh(2,1)*i_*g
       - d(DeltaR(1,2))*WRh(2,1)*DeltaRh(1,1)*i_*g
       + d(DeltaR(1,2))*WRh(2,1)*DeltaRh(2,2)*i_*g
       - d(DeltaR(1,2))*WRh(2,2)*DeltaRh(2,1)*i_*g
       + d(DeltaR(2,1))*B*DeltaRh(1,2)*i_*gp
       - d(DeltaR(2,1))*WRh(1,1)*DeltaRh(1,2)*i_*g
       + d(DeltaR(2,1))*WRh(1,2)*DeltaRh(1,1)*i_*g
       - d(DeltaR(2,1))*WRh(1,2)*DeltaRh(2,2)*i_*g
       + d(DeltaR(2,1))*WRh(2,2)*DeltaRh(1,2)*i_*g
       + d(DeltaR(2,2))*B*DeltaRh(2,2)*i_*gp
       + d(DeltaR(2,2))*WRh(1,2)*DeltaRh(2,1)*i_*g
       - d(DeltaR(2,2))*WRh(2,1)*DeltaRh(1,2)*i_*g
       + d(Phih(1,1))*d(Phi(1,1))
       - d(Phih(1,1))*WL(1,1)*Phi(1,1)*i_*g
       - d(Phih(1,1))*WL(1,2)*Phi(2,1)*i_*g
       + d(Phih(1,1))*WR(1,1)*Phi(1,1)*i_*g
       + d(Phih(1,1))*WR(2,1)*Phi(1,2)*i_*g
       + d(Phih(1,2))*d(Phi(2,1))
       - d(Phih(1,2))*WL(2,1)*Phi(1,1)*i_*g
       - d(Phih(1,2))*WL(2,2)*Phi(2,1)*i_*g
       + d(Phih(1,2))*WR(1,1)*Phi(2,1)*i_*g
       + d(Phih(1,2))*WR(2,1)*Phi(2,2)*i_*g
       + d(Phih(2,1))*d(Phi(1,2))
       - d(Phih(2,1))*WL(1,1)*Phi(1,2)*i_*g
       - d(Phih(2,1))*WL(1,2)*Phi(2,2)*i_*g
       + d(Phih(2,1))*WR(1,2)*Phi(1,1)*i_*g
       + d(Phih(2,1))*WR(2,2)*Phi(1,2)*i_*g
       + d(Phih(2,2))*d(Phi(2,2))
       - d(Phih(2,2))*WL(2,1)*Phi(1,2)*i_*g
       - d(Phih(2,2))*WL(2,2)*Phi(2,2)*i_*g
       + d(Phih(2,2))*WR(1,2)*Phi(2,1)*i_*g
       + d(Phih(2,2))*WR(2,2)*Phi(2,2)*i_*g
       + d(Phi(1,1))*WLh(1,1)*Phih(1,1)*i_*g
       + d(Phi(1,1))*WLh(2,1)*Phih(1,2)*i_*g
       - d(Phi(1,1))*WRh(1,1)*Phih(1,1)*i_*g
       - d(Phi(1,1))*WRh(1,2)*Phih(2,1)*i_*g
       + d(Phi(1,2))*WLh(1,1)*Phih(2,1)*i_*g
       + d(Phi(1,2))*WLh(2,1)*Phih(2,2)*i_*g
       - d(Phi(1,2))*WRh(2,1)*Phih(1,1)*i_*g
       - d(Phi(1,2))*WRh(2,2)*Phih(2,1)*i_*g
       + d(Phi(2,1))*WLh(1,2)*Phih(1,1)*i_*g
       + d(Phi(2,1))*WLh(2,2)*Phih(1,2)*i_*g
       - d(Phi(2,1))*WRh(1,1)*Phih(1,2)*i_*g
       - d(Phi(2,1))*WRh(1,2)*Phih(2,2)*i_*g
       + d(Phi(2,2))*WLh(1,2)*Phih(2,1)*i_*g
       + d(Phi(2,2))*WLh(2,2)*Phih(2,2)*i_*g
       - d(Phi(2,2))*WRh(2,1)*Phih(1,2)*i_*g
       - d(Phi(2,2))*WRh(2,2)*Phih(2,2)*i_*g
       + B^2*DeltaLh(1,1)*DeltaL(1,1)*gp^2
       + B^2*DeltaLh(1,2)*DeltaL(2,1)*gp^2
       + B^2*DeltaLh(2,1)*DeltaL(1,2)*gp^2
       + B^2*DeltaLh(2,2)*DeltaL(2,2)*gp^2
       + B^2*DeltaRh(1,1)*DeltaR(1,1)*gp^2
       + B^2*DeltaRh(1,2)*DeltaR(2,1)*gp^2
       + B^2*DeltaRh(2,1)*DeltaR(1,2)*gp^2
       + B^2*DeltaRh(2,2)*DeltaR(2,2)*gp^2
       - B*WLh(1,1)*DeltaLh(1,2)*DeltaL(2,1)*g*gp
       + B*WLh(1,1)*DeltaLh(2,1)*DeltaL(1,2)*g*gp
       + B*WLh(1,2)*DeltaLh(1,1)*DeltaL(2,1)*g*gp
       - B*WLh(1,2)*DeltaLh(2,1)*DeltaL(1,1)*g*gp
       + B*WLh(1,2)*DeltaLh(2,1)*DeltaL(2,2)*g*gp
       - B*WLh(1,2)*DeltaLh(2,2)*DeltaL(2,1)*g*gp
       - B*WLh(2,1)*DeltaLh(1,1)*DeltaL(1,2)*g*gp
       + B*WLh(2,1)*DeltaLh(1,2)*DeltaL(1,1)*g*gp
       - B*WLh(2,1)*DeltaLh(1,2)*DeltaL(2,2)*g*gp
       + B*WLh(2,1)*DeltaLh(2,2)*DeltaL(1,2)*g*gp
       + B*WLh(2,2)*DeltaLh(1,2)*DeltaL(2,1)*g*gp
       - B*WLh(2,2)*DeltaLh(2,1)*DeltaL(1,2)*g*gp
       - B*WL(1,1)*DeltaLh(1,2)*DeltaL(2,1)*g*gp
       + B*WL(1,1)*DeltaLh(2,1)*DeltaL(1,2)*g*gp
       + B*WL(1,2)*DeltaLh(1,1)*DeltaL(2,1)*g*gp
       - B*WL(1,2)*DeltaLh(2,1)*DeltaL(1,1)*g*gp
       + B*WL(1,2)*DeltaLh(2,1)*DeltaL(2,2)*g*gp
       - B*WL(1,2)*DeltaLh(2,2)*DeltaL(2,1)*g*gp
       - B*WL(2,1)*DeltaLh(1,1)*DeltaL(1,2)*g*gp
       + B*WL(2,1)*DeltaLh(1,2)*DeltaL(1,1)*g*gp
       - B*WL(2,1)*DeltaLh(1,2)*DeltaL(2,2)*g*gp
       + B*WL(2,1)*DeltaLh(2,2)*DeltaL(1,2)*g*gp
       + B*WL(2,2)*DeltaLh(1,2)*DeltaL(2,1)*g*gp
       - B*WL(2,2)*DeltaLh(2,1)*DeltaL(1,2)*g*gp
       - B*WRh(1,1)*DeltaRh(1,2)*DeltaR(2,1)*g*gp
       + B*WRh(1,1)*DeltaRh(2,1)*DeltaR(1,2)*g*gp
       + B*WRh(1,2)*DeltaRh(1,1)*DeltaR(2,1)*g*gp
       - B*WRh(1,2)*DeltaRh(2,1)*DeltaR(1,1)*g*gp
       + B*WRh(1,2)*DeltaRh(2,1)*DeltaR(2,2)*g*gp
       - B*WRh(1,2)*DeltaRh(2,2)*DeltaR(2,1)*g*gp
       - B*WRh(2,1)*DeltaRh(1,1)*DeltaR(1,2)*g*gp
       + B*WRh(2,1)*DeltaRh(1,2)*DeltaR(1,1)*g*gp
       - B*WRh(2,1)*DeltaRh(1,2)*DeltaR(2,2)*g*gp
       + B*WRh(2,1)*DeltaRh(2,2)*DeltaR(1,2)*g*gp
       + B*WRh(2,2)*DeltaRh(1,2)*DeltaR(2,1)*g*gp
       - B*WRh(2,2)*DeltaRh(2,1)*DeltaR(1,2)*g*gp
       - B*WR(1,1)*DeltaRh(1,2)*DeltaR(2,1)*g*gp
       + B*WR(1,1)*DeltaRh(2,1)*DeltaR(1,2)*g*gp
       + B*WR(1,2)*DeltaRh(1,1)*DeltaR(2,1)*g*gp
       - B*WR(1,2)*DeltaRh(2,1)*DeltaR(1,1)*g*gp
       + B*WR(1,2)*DeltaRh(2,1)*DeltaR(2,2)*g*gp
       - B*WR(1,2)*DeltaRh(2,2)*DeltaR(2,1)*g*gp
       - B*WR(2,1)*DeltaRh(1,1)*DeltaR(1,2)*g*gp
       + B*WR(2,1)*DeltaRh(1,2)*DeltaR(1,1)*g*gp
       - B*WR(2,1)*DeltaRh(1,2)*DeltaR(2,2)*g*gp
       + B*WR(2,1)*DeltaRh(2,2)*DeltaR(1,2)*g*gp
       + B*WR(2,2)*DeltaRh(1,2)*DeltaR(2,1)*g*gp
       - B*WR(2,2)*DeltaRh(2,1)*DeltaR(1,2)*g*gp
       + WLh(1,1)*WL(1,1)*DeltaLh(1,2)*DeltaL(2,1)*g^2
       + WLh(1,1)*WL(1,1)*DeltaLh(2,1)*DeltaL(1,2)*g^2
       + WLh(1,1)*WL(1,1)*Phih(1,1)*Phi(1,1)*g^2
       + WLh(1,1)*WL(1,1)*Phih(2,1)*Phi(1,2)*g^2
       - WLh(1,1)*WL(1,2)*DeltaLh(2,1)*DeltaL(1,1)*g^2
       + WLh(1,1)*WL(1,2)*DeltaLh(2,1)*DeltaL(2,2)*g^2
       + WLh(1,1)*WL(1,2)*Phih(1,1)*Phi(2,1)*g^2
       + WLh(1,1)*WL(1,2)*Phih(2,1)*Phi(2,2)*g^2
       - WLh(1,1)*WL(2,1)*DeltaLh(1,2)*DeltaL(1,1)*g^2
       + WLh(1,1)*WL(2,1)*DeltaLh(1,2)*DeltaL(2,2)*g^2
       - WLh(1,1)*WL(2,2)*DeltaLh(1,2)*DeltaL(2,1)*g^2
       - WLh(1,1)*WL(2,2)*DeltaLh(2,1)*DeltaL(1,2)*g^2
       - WLh(1,1)*WR(1,1)*Phih(1,1)*Phi(1,1)*g^2
       - WLh(1,1)*WR(1,2)*Phih(2,1)*Phi(1,1)*g^2
       - WLh(1,1)*WR(2,1)*Phih(1,1)*Phi(1,2)*g^2
       - WLh(1,1)*WR(2,2)*Phih(2,1)*Phi(1,2)*g^2
       - WLh(1,2)*WL(1,1)*DeltaLh(1,1)*DeltaL(2,1)*g^2
       + WLh(1,2)*WL(1,1)*DeltaLh(2,2)*DeltaL(2,1)*g^2
       - 2*WLh(1,2)*WL(1,2)*DeltaLh(2,1)*DeltaL(2,1)*g^2
       + WLh(1,2)*WL(2,1)*DeltaLh(1,1)*DeltaL(1,1)*g^2
       - WLh(1,2)*WL(2,1)*DeltaLh(1,1)*DeltaL(2,2)*g^2
       + 2*WLh(1,2)*WL(2,1)*DeltaLh(2,1)*DeltaL(1,2)*g^2
       - WLh(1,2)*WL(2,1)*DeltaLh(2,2)*DeltaL(1,1)*g^2
       + WLh(1,2)*WL(2,1)*DeltaLh(2,2)*DeltaL(2,2)*g^2
       + WLh(1,2)*WL(2,1)*Phih(1,1)*Phi(1,1)*g^2
       + WLh(1,2)*WL(2,1)*Phih(2,1)*Phi(1,2)*g^2
       + WLh(1,2)*WL(2,2)*DeltaLh(1,1)*DeltaL(2,1)*g^2
       - WLh(1,2)*WL(2,2)*DeltaLh(2,2)*DeltaL(2,1)*g^2
       + WLh(1,2)*WL(2,2)*Phih(1,1)*Phi(2,1)*g^2
       + WLh(1,2)*WL(2,2)*Phih(2,1)*Phi(2,2)*g^2
       - WLh(1,2)*WR(1,1)*Phih(1,1)*Phi(2,1)*g^2
       - WLh(1,2)*WR(1,2)*Phih(2,1)*Phi(2,1)*g^2
       - WLh(1,2)*WR(2,1)*Phih(1,1)*Phi(2,2)*g^2
       - WLh(1,2)*WR(2,2)*Phih(2,1)*Phi(2,2)*g^2
       - WLh(2,1)*WL(1,1)*DeltaLh(1,1)*DeltaL(1,2)*g^2
       + WLh(2,1)*WL(1,1)*DeltaLh(2,2)*DeltaL(1,2)*g^2
       + WLh(2,1)*WL(1,1)*Phih(1,2)*Phi(1,1)*g^2
       + WLh(2,1)*WL(1,1)*Phih(2,2)*Phi(1,2)*g^2
       + WLh(2,1)*WL(1,2)*DeltaLh(1,1)*DeltaL(1,1)*g^2
       - WLh(2,1)*WL(1,2)*DeltaLh(1,1)*DeltaL(2,2)*g^2
       + 2*WLh(2,1)*WL(1,2)*DeltaLh(1,2)*DeltaL(2,1)*g^2
       - WLh(2,1)*WL(1,2)*DeltaLh(2,2)*DeltaL(1,1)*g^2
       + WLh(2,1)*WL(1,2)*DeltaLh(2,2)*DeltaL(2,2)*g^2
       + WLh(2,1)*WL(1,2)*Phih(1,2)*Phi(2,1)*g^2
       + WLh(2,1)*WL(1,2)*Phih(2,2)*Phi(2,2)*g^2
       - 2*WLh(2,1)*WL(2,1)*DeltaLh(1,2)*DeltaL(1,2)*g^2
       + WLh(2,1)*WL(2,2)*DeltaLh(1,1)*DeltaL(1,2)*g^2
       - WLh(2,1)*WL(2,2)*DeltaLh(2,2)*DeltaL(1,2)*g^2
       - WLh(2,1)*WR(1,1)*Phih(1,2)*Phi(1,1)*g^2
       - WLh(2,1)*WR(1,2)*Phih(2,2)*Phi(1,1)*g^2
       - WLh(2,1)*WR(2,1)*Phih(1,2)*Phi(1,2)*g^2
       - WLh(2,1)*WR(2,2)*Phih(2,2)*Phi(1,2)*g^2
       - WLh(2,2)*WL(1,1)*DeltaLh(1,2)*DeltaL(2,1)*g^2
       - WLh(2,2)*WL(1,1)*DeltaLh(2,1)*DeltaL(1,2)*g^2
       + WLh(2,2)*WL(1,2)*DeltaLh(2,1)*DeltaL(1,1)*g^2
       - WLh(2,2)*WL(1,2)*DeltaLh(2,1)*DeltaL(2,2)*g^2
       + WLh(2,2)*WL(2,1)*DeltaLh(1,2)*DeltaL(1,1)*g^2
       - WLh(2,2)*WL(2,1)*DeltaLh(1,2)*DeltaL(2,2)*g^2
       + WLh(2,2)*WL(2,1)*Phih(1,2)*Phi(1,1)*g^2
       + WLh(2,2)*WL(2,1)*Phih(2,2)*Phi(1,2)*g^2
       + WLh(2,2)*WL(2,2)*DeltaLh(1,2)*DeltaL(2,1)*g^2
       + WLh(2,2)*WL(2,2)*DeltaLh(2,1)*DeltaL(1,2)*g^2
       + WLh(2,2)*WL(2,2)*Phih(1,2)*Phi(2,1)*g^2
       + WLh(2,2)*WL(2,2)*Phih(2,2)*Phi(2,2)*g^2
       - WLh(2,2)*WR(1,1)*Phih(1,2)*Phi(2,1)*g^2
       - WLh(2,2)*WR(1,2)*Phih(2,2)*Phi(2,1)*g^2
       - WLh(2,2)*WR(2,1)*Phih(1,2)*Phi(2,2)*g^2
       - WLh(2,2)*WR(2,2)*Phih(2,2)*Phi(2,2)*g^2
       - WL(1,1)*WRh(1,1)*Phih(1,1)*Phi(1,1)*g^2
       - WL(1,1)*WRh(1,2)*Phih(2,1)*Phi(1,1)*g^2
       - WL(1,1)*WRh(2,1)*Phih(1,1)*Phi(1,2)*g^2
       - WL(1,1)*WRh(2,2)*Phih(2,1)*Phi(1,2)*g^2
       - WL(1,2)*WRh(1,1)*Phih(1,1)*Phi(2,1)*g^2
       - WL(1,2)*WRh(1,2)*Phih(2,1)*Phi(2,1)*g^2
       - WL(1,2)*WRh(2,1)*Phih(1,1)*Phi(2,2)*g^2
       - WL(1,2)*WRh(2,2)*Phih(2,1)*Phi(2,2)*g^2
       - WL(2,1)*WRh(1,1)*Phih(1,2)*Phi(1,1)*g^2
       - WL(2,1)*WRh(1,2)*Phih(2,2)*Phi(1,1)*g^2
       - WL(2,1)*WRh(2,1)*Phih(1,2)*Phi(1,2)*g^2
       - WL(2,1)*WRh(2,2)*Phih(2,2)*Phi(1,2)*g^2
       - WL(2,2)*WRh(1,1)*Phih(1,2)*Phi(2,1)*g^2
       - WL(2,2)*WRh(1,2)*Phih(2,2)*Phi(2,1)*g^2
       - WL(2,2)*WRh(2,1)*Phih(1,2)*Phi(2,2)*g^2
       - WL(2,2)*WRh(2,2)*Phih(2,2)*Phi(2,2)*g^2
       + WRh(1,1)*WR(1,1)*DeltaRh(1,2)*DeltaR(2,1)*g^2
       + WRh(1,1)*WR(1,1)*DeltaRh(2,1)*DeltaR(1,2)*g^2
       + WRh(1,1)*WR(1,1)*Phih(1,1)*Phi(1,1)*g^2
       + WRh(1,1)*WR(1,1)*Phih(1,2)*Phi(2,1)*g^2
       - WRh(1,1)*WR(1,2)*DeltaRh(2,1)*DeltaR(1,1)*g^2
       + WRh(1,1)*WR(1,2)*DeltaRh(2,1)*DeltaR(2,2)*g^2
       - WRh(1,1)*WR(2,1)*DeltaRh(1,2)*DeltaR(1,1)*g^2
       + WRh(1,1)*WR(2,1)*DeltaRh(1,2)*DeltaR(2,2)*g^2
       + WRh(1,1)*WR(2,1)*Phih(1,1)*Phi(1,2)*g^2
       + WRh(1,1)*WR(2,1)*Phih(1,2)*Phi(2,2)*g^2
       - WRh(1,1)*WR(2,2)*DeltaRh(1,2)*DeltaR(2,1)*g^2
       - WRh(1,1)*WR(2,2)*DeltaRh(2,1)*DeltaR(1,2)*g^2
       - WRh(1,2)*WR(1,1)*DeltaRh(1,1)*DeltaR(2,1)*g^2
       + WRh(1,2)*WR(1,1)*DeltaRh(2,2)*DeltaR(2,1)*g^2
       + WRh(1,2)*WR(1,1)*Phih(2,1)*Phi(1,1)*g^2
       + WRh(1,2)*WR(1,1)*Phih(2,2)*Phi(2,1)*g^2
       - 2*WRh(1,2)*WR(1,2)*DeltaRh(2,1)*DeltaR(2,1)*g^2
       + WRh(1,2)*WR(2,1)*DeltaRh(1,1)*DeltaR(1,1)*g^2
       - WRh(1,2)*WR(2,1)*DeltaRh(1,1)*DeltaR(2,2)*g^2
       + 2*WRh(1,2)*WR(2,1)*DeltaRh(2,1)*DeltaR(1,2)*g^2
       - WRh(1,2)*WR(2,1)*DeltaRh(2,2)*DeltaR(1,1)*g^2
       + WRh(1,2)*WR(2,1)*DeltaRh(2,2)*DeltaR(2,2)*g^2
       + WRh(1,2)*WR(2,1)*Phih(2,1)*Phi(1,2)*g^2
       + WRh(1,2)*WR(2,1)*Phih(2,2)*Phi(2,2)*g^2
       + WRh(1,2)*WR(2,2)*DeltaRh(1,1)*DeltaR(2,1)*g^2
       - WRh(1,2)*WR(2,2)*DeltaRh(2,2)*DeltaR(2,1)*g^2
       - WRh(2,1)*WR(1,1)*DeltaRh(1,1)*DeltaR(1,2)*g^2
       + WRh(2,1)*WR(1,1)*DeltaRh(2,2)*DeltaR(1,2)*g^2
       + WRh(2,1)*WR(1,2)*DeltaRh(1,1)*DeltaR(1,1)*g^2
       - WRh(2,1)*WR(1,2)*DeltaRh(1,1)*DeltaR(2,2)*g^2
       + 2*WRh(2,1)*WR(1,2)*DeltaRh(1,2)*DeltaR(2,1)*g^2
       - WRh(2,1)*WR(1,2)*DeltaRh(2,2)*DeltaR(1,1)*g^2
       + WRh(2,1)*WR(1,2)*DeltaRh(2,2)*DeltaR(2,2)*g^2
       + WRh(2,1)*WR(1,2)*Phih(1,1)*Phi(1,1)*g^2
       + WRh(2,1)*WR(1,2)*Phih(1,2)*Phi(2,1)*g^2
       - 2*WRh(2,1)*WR(2,1)*DeltaRh(1,2)*DeltaR(1,2)*g^2
       + WRh(2,1)*WR(2,2)*DeltaRh(1,1)*DeltaR(1,2)*g^2
       - WRh(2,1)*WR(2,2)*DeltaRh(2,2)*DeltaR(1,2)*g^2
       + WRh(2,1)*WR(2,2)*Phih(1,1)*Phi(1,2)*g^2
       + WRh(2,1)*WR(2,2)*Phih(1,2)*Phi(2,2)*g^2
       - WRh(2,2)*WR(1,1)*DeltaRh(1,2)*DeltaR(2,1)*g^2
       - WRh(2,2)*WR(1,1)*DeltaRh(2,1)*DeltaR(1,2)*g^2
       + WRh(2,2)*WR(1,2)*DeltaRh(2,1)*DeltaR(1,1)*g^2
       - WRh(2,2)*WR(1,2)*DeltaRh(2,1)*DeltaR(2,2)*g^2
       + WRh(2,2)*WR(1,2)*Phih(2,1)*Phi(1,1)*g^2
       + WRh(2,2)*WR(1,2)*Phih(2,2)*Phi(2,1)*g^2
       + WRh(2,2)*WR(2,1)*DeltaRh(1,2)*DeltaR(1,1)*g^2
       - WRh(2,2)*WR(2,1)*DeltaRh(1,2)*DeltaR(2,2)*g^2
       + WRh(2,2)*WR(2,2)*DeltaRh(1,2)*DeltaR(2,1)*g^2
       + WRh(2,2)*WR(2,2)*DeltaRh(2,1)*DeltaR(1,2)*g^2
       + WRh(2,2)*WR(2,2)*Phih(2,1)*Phi(1,2)*g^2
       + WRh(2,2)*WR(2,2)*Phih(2,2)*Phi(2,2)*g^2
      ;

    ************************************************************
    
    s sq2;
    s ka1,ka2, vr,kap;
    cf wlp, wlm, wrp, wrm, wl3, wr3;
    cf deltal0, deltalp, deltalpp;
    cf deltal0h, deltalm, deltalmm;
    cf deltar0, deltarp, deltarpp;
    cf deltar0h, deltarm, deltarmm;
    cf phi1p, phi1m, phi2p, phi2m, phi10, phi10h, phi20, phi20h;
    
    id WLh(i?, j?) = WL(i, j);
    id WRh(i?, j?) = WR(i, j);
    
    id WL(1, 1) = wl3/2;
    id WL(1, 2) = wlp/sq2;
    id WL(2, 1) = wlm/sq2;
    id WL(2, 2) = -wl3/2;
    
    id WR(1, 1) = wr3/2;
    id WR(1, 2) = wrp/sq2;
    id WR(2, 1) = wrm/sq2;
    id WR(2, 2) = -wr3/2;
    
    id DeltaL(1, 1) = deltalp/sq2;
    id DeltaL(1, 2) = deltalpp;
    id DeltaL(2, 1) = deltal0;
    id DeltaL(2, 2) = -deltalp/sq2;
    
    id DeltaLh(1, 1) = deltalm/sq2;
    id DeltaLh(1, 2) = deltal0h;
    id DeltaLh(2, 1) = deltalmm;
    id DeltaLh(2, 2) = -deltalm/sq2;
    
    id DeltaR(1, 1) = deltarp/sq2;
    id DeltaR(1, 2) = deltarpp;
    id DeltaR(2, 1) = vr/sq2+deltar0;
    id DeltaR(2, 2) = -deltarp/sq2;
    
    id DeltaRh(1, 1) = deltarm/sq2;
    id DeltaRh(1, 2) = vr/sq2+deltar0h;
    id DeltaRh(2, 1) = deltarmm;
    id DeltaRh(2, 2) = -deltarm/sq2;
    
    id Phi(1, 1) = phi10+ka1/sq2;
    id Phi(1, 2) = phi1p;
    id Phi(2, 1) = phi2m;
    id Phi(2, 2) = phi20+ka2/sq2;
    
    id Phih(1, 1) = phi10h+ka1/sq2;
    id Phih(1, 2) = phi2p;
    id Phih(2, 1) = phi1m;
    id Phih(2, 2) = phi20h+ka2/sq2;
    
    id d(DeltaL(1, 1)) = d(deltalp)/sq2;
    id d(DeltaL(1, 2)) = d(deltalpp);
    id d(DeltaL(2, 1)) = d(deltal0);
    id d(DeltaL(2, 2)) = -d(deltalp)/sq2;
    
    id d(DeltaLh(1, 1)) = d(deltalm)/sq2;
    id d(DeltaLh(1, 2)) = d(deltal0h);
    id d(DeltaLh(2, 1)) = d(deltalmm);
    id d(DeltaLh(2, 2)) = -d(deltalm)/sq2;
    
    id d(DeltaR(1, 1)) = d(deltarp)/sq2;
    id d(DeltaR(1, 2)) = d(deltarpp);
    id d(DeltaR(2, 1)) = d(deltar0);
    id d(DeltaR(2, 2)) = -d(deltarp)/sq2;
    
    id d(DeltaRh(1, 1)) = d(deltarm)/sq2;
    id d(DeltaRh(1, 2)) = d(deltar0h);
    id d(DeltaRh(2, 1)) = d(deltarmm);
    id d(DeltaRh(2, 2)) = -d(deltarm)/sq2;
    
    id d(Phi(1, 1)) = d(phi10);
    id d(Phi(1, 2)) = d(phi1p);
    id d(Phi(2, 1)) = d(phi2m);
    id d(Phi(2, 2)) = d(phi20);
    
    id d(Phih(1, 1)) = d(phi10h);
    id d(Phih(1, 2)) = d(phi2p);
    id d(Phih(2, 1)) = d(phi1m);
    id d(Phih(2, 2)) = d(phi20h);
    
    id 1/sq2^2 = 1/2;
    
    *print +s;
    *.end
    .sort
    
    
    set higgsy: deltal0, deltalp, deltalpp,deltal0h, deltalm, deltalmm,deltar0, delt
    arp, deltarpp,
    deltar0h, deltarm, deltarmm,phi1p, phi1m, phi2p, phi2m, phi10, phi10h, phi20, ph
    i20h;
    
    #ifdef `masy'
    id deltalpp?higgsy =0;
    if (match(d(deltalp?higgsy)) !=0);
    discard;
    endif;
    b wr3,wl3,B,wrp,wrm,wlp,wlm;
    print +s;
    *.end
    #else
    #endif
    .sort

   exp =

       + B^2 * (
          + 1/2*gp^2*vr^2
          )

       + B*wr3 * (
          - g*gp*vr^2
          )

       + wlp*wlm * (
          + 1/4*g^2*ka2^2
          + 1/4*g^2*ka1^2
          )

       + wlp*wrm * (
          - 1/2*g^2*ka1*ka2
          )

       + wlm*wrp * (
          - 1/2*g^2*ka1*ka2
          )

       + wrp*wrm * (
          + 1/2*g^2*vr^2
          + 1/4*g^2*ka2^2
          + 1/4*g^2*ka1^2
          )

       + wl3^2 * (
          + 1/8*g^2*ka2^2
          + 1/8*g^2*ka1^2
          )

       + wl3*wr3 * (
          - 1/4*g^2*ka2^2
          - 1/4*g^2*ka1^2
          )

       + wr3^2 * (
          + 1/2*g^2*vr^2
          + 1/8*g^2*ka2^2
          + 1/8*g^2*ka1^2
          );

    .end