Schoonschip, 68000 version of June 27, 1991. Public version.
Date: Fri Aug  9 1991 19:58:46.  Memory: start 0001B84C, length 476348.

Command line: Titard.e Ti_out

C Takes about 22 sec.

C A test of N-traces (Stephan Titard).

C Counter-terms in CDR (dim. regularization)
	C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
C Compute d3,d2,dm,dg,d2c
C Z3=1+d3*g^2 ...
C counterterms(in p-space): w/ rescaling and Zksi=Z3
C d3*g^2*(D(mu,nu)*pDp-p(mu)*p(nu)) --> from gluon self-energy
C d1F*g^2*(i*g*PSIb*gam(mu)*PSI*G(mu)) -->get d1F from quark-quark-gluon vertex
C                                        and dg using  d1F=d2+1/2*d3+dg  
C d2*g^2*(PSIb*gam(p)*PSI)  
C (d2+dm)*g^2*(PSIb*m*PSI) --> from quark self-energy
C  for Nf flavours  T(R) --> Nf*T(R)
	CCCCC
> P ni

Begin. Time 2 sec.
> P input
	Z vertg=PTF3("g,"g,"g)
> TOPO1L{}
	C
	V P1,P2,P3,K,U1,U2,U3
L 1	Id,SELF(I1~,I2~)=
	    DS(I1;-J1;J2;Sym;-J1;J2;(
	    DC("C,TC,I2,-J2,J1)*DC("F,TFE,-1,J1,J2)*DC("F,TGH,-1,J1,J2)*
	    DIA(I1,I2,-J1,J1,-J2,J2) ))
	   +DS(I1;-J3;J3;I2;Sym;-J3;J3;(
	    DIB(I1;I2;-J3,J3) ))
L 3	Id,Anti,TAP
C             
L 4	Id,DIA(I1~,I2~,K4~,K5~,K6~,K7~)=
	   VE3(I1,K4,K7,*,MU4,s1,c1,a1,P1,*,L4,s4,c4,a4,-U1,*,L7,s7,c7,a7,K)*
	   VE3(I2,K6,K5,*,NU4,s2,c2,a2,-P1,*,L6,s6,c6,a6,-K,*,L5,s5,c5,a5,U1)*
	   PROP(K5,K4,*,L5,s5,c5,a5,U1,*,L4,s4,c4,a4,-U1)*
	   PROP(K7,K6,*,L7,s7,c7,a7,K,*,L6,s6,c6,a6,-K)
L 4	Al,DIB(I1~,I2~,K4~,K5~)=
	   VE4(I1,K4,K5,I2,*,MU4,a1,*,L4,a4,*,L5,a5,*,NU4,a2)*
	   PROP(K5,K4,*,L5,s5,c5,a5,K,*,L4,s4,c4,a4,-K)
	C
L 5	Id,PTF3(I1~,I2~,I3~)=
	    DS(I1;-J1;J3;(
	    DS(I3;-J2;J1;(
	    DC("C,TC,I2,-J3,J2)*DC("F,TFE,-1,J1,J2,J3)*DC("F,TGH,-1,J1,J2,J3)*
	    DIC(I1,I2,I3,-J1,J1,-J2,J2,-J3,J3) )) ))
	   +DS(I3;I2;-J5;J4;Sym;-J5;J4;(
	    DC("C,TC,I1,-J4,J5)*DC("F,TGH,-1,J4,J5)*
	    VIR1(I1,I2,I3,-J4,J4,-J5,J5) ))
	   +DS(I1;I3;-J7;J6;Sym;-J7;J6;(
	    DC("C,TC,I2,-J6,J7)*DC("F,TGH,-1,J6,J7)*
	    VIR2(I1,I2,I3,-J6,J6,-J7,J7) ))
	   +DS(I2;I1;-J8;J9;Sym;-J8;J9;(
	    DC("C,TC,I3,-J9,J8)*DC("F,TGH,-1,J8,J9)*
	    VIR3(I1,I2,I3,-J8,J8,-J9,J9) ))
L 8	Id,Anti,TAP
L 9	Id,DIC(I1~,I2~,I3~,K4~,K5~,K6~,K7~,K8~,K9~)=
	   VE3(I1,K4,K9,*,MU4,s1,c1,a1,P1,*,L4,s4,c4,a4,K,*,L9,s9,c9,a9,-U1)*
	   VE3(I2,K8,K7,*,NU4,s2,c2,a2,P2,*,L8,s8,c8,a8,U1,*,L7,s7,c7,a7,-U3)*
	   VE3(I3,K6,K5,*,LA4,s3,c3,a3,P3,*,L6,s6,c6,a6,U3,*,L5,s5,c5,a5,-K)*
	   PROP(K5,K4,*,L5,s5,c5,a5,-K,*,L4,s4,c4,a4,K)*     
	   PROP(K7,K6,*,L7,s7,c7,a7,-U3,*,L6,s6,c6,a6,U3)*
	   PROP(K9,K8,*,L9,s9,c9,a9,-U1,*,L8,s8,c8,a8,U1)
C                              -(P)
C we used inverse conv here  -J ---> J  bec/ we change the internal vectors.
L10	Id,VIR1(I1~,I2~,I3~,K4~,K5~,K6~,K7~)=
	   VE4(I3,I2,K6,K5,*,LA4,a3,*,NU4,a2,*,L6,a6,*,L5,a5)*
	   VE3(I1,K4,K7,*,MU4,s1,c1,a1,P1,*,L4,s4,c4,a4,-U1,*,L7,s7,c7,a7,K)*
	   PROP(K5,K4,*,L5,s5,c5,a5,U1,*,L4,s4,c4,a4,-U1)*
	   PROP(K7,K6,*,L7,s7,c7,a7,K,*,L6,s6,c6,a6,-K)
L10	Al,VIR2(I1~,I2~,I3~,K4~,K5~,K6~,K7~)=    
	   VE4(I1,I3,K6,K5,*,MU4,a1,*,LA4,a3,*,L6,a6,*,L5,a5)*
	   VE3(I2,K4,K7,*,NU4,s2,c2,a2,P2,*,L4,s4,c4,a4,-U2,*,L7,s7,c7,a7,K)* 
	   PROP(K5,K4,*,L5,s5,c5,a5,U2,*,L4,s4,c4,a4,-U2)*
	   PROP(K7,K6,*,L7,s7,c7,a7,K,*,L6,s6,c6,a6,-K)
L10	Al,VIR3(I1~,I2~,I3~,K4~,K5~,K6~,K7~)=
	   VE4(I2,I1,K4,K7,*,NU4,a2,*,MU4,a1,*,L4,a4,*,L7,a7)*
	   VE3(I3,K6,K5,*,LA4,s3,c3,a3,P3,*,L6,s6,c6,a6,-K,*,L5,s5,c5,a5,U3)*
	   PROP(K5,K4,*,L5,s5,c5,a5,U3,*,L4,s4,c4,a4,-U3)*
	   PROP(K7,K6,*,L7,s7,c7,a7,K,*,L6,s6,c6,a6,-K)
C U1=K+P1, U2=K+P2, U3=K-P3=K+P1+P2
	C
L11	Id,PTS3(I1~,I2~,I3~)=
	    DS(I1;-J1;J3;(
	    DS(I3;-J2;J1;(
	    DC("C,TC,I2,-J3,J2)*DC("F,TFE,-1,J1,J2,J3)*DC("F,TGH,-1,J1,J2,J3)*
	    DIC(I1,I2,I3,-J1,J1,-J2,J2,-J3,J3) )) ))
	   +DS(I3;I2;-J5;J4;Sym;-J5;J4;(
	    DC("C,TC,I1,-J4,J5)*DC("F,TGH,-1,J4,J5)*
	    VIR1(I1,I2,I3,-J4,J4,-J5,J5) ))
	   +DS(I1;I3;-J7;J6;Sym;-J7;J6;(
	    DC("C,TC,I2,-J6,J7)*DC("F,TGH,-1,J6,J7)*
	    VIR2(I1,I2,I3,-J6,J6,-J7,J7) ))
	   +DS(I2;I1;-J8;J9;Sym;-J8;J9;(
	    DC("C,TC,I3,-J9,J8)*DC("F,TGH,-1,J8,J9)*
	    VIR3(I1,I2,I3,-J8,J8,-J9,J9) ))
L14	Id,Anti,TAP
L15	Id,DIC(I1~,I2~,I3~,K4~,K5~,K6~,K7~,K8~,K9~)=
	   VE3(I1,K4,K9,*,MU4,s1,c1,a1,P1,*,L4,s4,c4,a4,K,*,L9,s9,c9,a9,-U1)*
	   VE3(I2,K8,K7,*,NU4,s2,c2,a2,P2,*,L8,s8,c8,a8,U1,*,L7,s7,c7,a7,-U3)*
	   VE3(I3,K6,K5,*,LA4,s3,c3,a3,P3,*,L6,s6,c6,a6,U3,*,L5,s5,c5,a5,-K)*
	   PROP(K5,K4,*,L5,s5,c5,a5,-K,*,L4,s4,c4,a4,K)*     
	   PROP(K7,K6,*,L7,s7,c7,a7,-U3,*,L6,s6,c6,a6,U3)*
	   PROP(K9,K8,*,L9,s9,c9,a9,-U1,*,L8,s8,c8,a8,U1)
C                              -(P)
C we used inverse conv here  -J ---> J  bec/ we change the internal vectors.
L16	Id,VIR1(I1~,I2~,I3~,K4~,K5~,K6~,K7~)=
	   VE4p(I3,I2,K6,K5,*,LA4,a3,P3,*,NU4,a2,P2,*,L6,a6,-K,*,L5,a5,U1)*
	   VE3(I1,K4,K7,*,MU4,s1,c1,a1,P1,*,L4,s4,c4,a4,-U1,*,L7,s7,c7,a7,K)*
	   PROP(K5,K4,*,L5,s5,c5,a5,U1,*,L4,s4,c4,a4,-U1)*
	   PROP(K7,K6,*,L7,s7,c7,a7,K,*,L6,s6,c6,a6,-K)
L16	Al,VIR2(I1~,I2~,I3~,K4~,K5~,K6~,K7~)=    
	   VE4p(I1,I3,K6,K5,*,MU4,a1,P1,*,LA4,a3,P3,*,L6,a6,-K,*,L5,a5,U2)*
	   VE3(I2,K4,K7,*,NU4,s2,c2,a2,P2,*,L4,s4,c4,a4,-U2,*,L7,s7,c7,a7,K)* 
	   PROP(K5,K4,*,L5,s5,c5,a5,U2,*,L4,s4,c4,a4,-U2)*
	   PROP(K7,K6,*,L7,s7,c7,a7,K,*,L6,s6,c6,a6,-K)
L16	Al,VIR3(I1~,I2~,I3~,K4~,K5~,K6~,K7~)=
	   VE4p(I2,I1,K4,K7,*,NU4,a2,P2,*,MU4,a1,P1,*,L4,a4,-U3,*,L7,a7,K)*
	   VE3(I3,K6,K5,*,LA4,s3,c3,a3,P3,*,L6,s6,c6,a6,-K,*,L5,s5,c5,a5,U3)*
	   PROP(K5,K4,*,L5,s5,c5,a5,U3,*,L4,s4,c4,a4,-U3)*
	   PROP(K7,K6,*,L7,s7,c7,a7,K,*,L6,s6,c6,a6,-K)
C U1=K+P1, U2=K+P2, U3=K-P3=K+P1+P2
	C
> COMPO{}
	*yep
L 1	Id,Compo,<X>,VE4,VE4p,SVE4,VE3,PROP
L 2	Id,VE4(gg~,L4~,a4~,L5~,a5~,L6~,a6~,L7~,a7~)=gg(L4,a4,L5,a5,L6,a6,L7,a7)
L 2	Al,VE4p(gg~,L4~,a4~,K~,L5~,a5~,P1~,L6~,a6~,P2~,L7~,a7~,P3~)=
	      gg(L4,a4,K,L5,a5,P1,L6,a6,P2,L7,a7,P3)
L 2	Al,SVE4(gq_q1~,MU4~,s1~,c1~,a1~,P~,L4~,s4~,c4~,a4~,K~,L5~,s5~,c5~,a5~,Q~)=
	        gq_q1(MU4,s1,c1,a1,P,L4,s4,c4,a4,K,L5,s5,c5,a5,Q)
L 2	Al,VE3(gg~,MU4~,s1~,c1~,a1~,K~,L4~,s4~,c4~,a4~,Q~,L5~,s5~,c5~,a5~,P~)=
	       gg(MU4,s1,c1,a1,K,L4,s4,c4,a4,Q,L5,s5,c5,a5,P)
L 2	Al,PROP(gg~,L4~,s4~,c4~,a4~,K~,L5~,s5~,c5~,a5~,P~)=
	        gg(L4,s4,c4,a4,L5,s5,c5,a5,P)
L 6	Id,flagA=1
	I Na1=A,Na2=A
	I Na1,Na2
	DO ii=1,2
L 7	Id,f(a0~,a1~,a2~,a3~)=flagA*f(a0,a1,a2,a3)
L 8	Id,flagA^n~*f(a0~,a1~,a2~,a3~)=f(a0,a1,Na1)*f(a2,a3,Na1)
	ENDDO
	Sum,Na1,Na2
	*yep
L 1	Id,Even,nom,1
L 2	Id,Commu,nom,Rf,f
	*yep
> GROUPF{}
	*yep
C exterior colour indices: a_0-3,c_0-3
L 1	Id,flagA=1
C group factor
L 2	Id,Asymm,f,1,2,3
C first products of 3 fs
L 3	Id,Adiso,f(a1~,a4~,a5~)*f(a2~,a5~,a6~)*f(a3~,a6~,a4~)=Nc/2*f(a1,a2,a3)
L 3	Al,Adiso,f(a1~,a5~,a4~)*f(a2~,a5~,a6~)*f(a3~,a6~,a4~)=-Nc/2*f(a1,a2,a3)
L 3	Al,Adiso,f(a1~,a4~,a5~)*f(a2~,a6~,a5~)*f(a3~,a6~,a4~)=-Nc/2*f(a1,a2,a3)
L 3	Al,Adiso,f(a1~,a4~,a5~)*f(a2~,a5~,a6~)*f(a3~,a4~,a6~)=-Nc/2*f(a1,a2,a3)
L 3	Al,Adiso,f(a1~,a5~,a4~)*f(a2~,a6~,a5~)*f(a3~,a6~,a4~)=Nc/2*f(a1,a2,a3)
L 3	Al,Adiso,f(a1~,a5~,a4~)*f(a2~,a5~,a6~)*f(a3~,a4~,a6~)=Nc/2*f(a1,a2,a3)
L 3	Al,Adiso,f(a1~,a4~,a5~)*f(a2~,a6~,a5~)*f(a3~,a4~,a6~)=Nc/2*f(a1,a2,a3)
L 4	Id,Nc=CAS2(Gr)
L 5	Id,Asymm,f,1,2,3
	C
	DO ii=1,3
L 6	Id,Adiso,Rf(a4~,c4~,c5~)*Rf(a5~,c5~,c4~)=Tg(Rep)*D(a4,a5)
L 6	Al,Adiso,Rf(a4~,c4~,c5~)*Rf(a4~,c5~,c6~)=CAS2(Rep)*D(c4,c6)
L 6	Al,Rf(a4~,c4~,c4~)=0
	ENDDO
	C
L 9	IF Adiso,Rf(a4~,c4~,c5~)*Rf(a5~,c5~,c6~)*Rf(a4~,c6~,c7~)*Rf(a6~,c7~,c4~)
	OR Adiso,Rf(a5~,c4~,c5~)*Rf(a4~,c5~,c6~)*Rf(a6~,c6~,c7~)*Rf(a4~,c7~,c4~)
L 9	Al,Addfa,flagA
	ENDIF
C 
	DO ii=1,3
L10	Id,Adiso,Rf(a4~,c4~,c5~)*Rf(a5~,c5~,c6~)*Rf(a4~,c6~,c7~)=
	         Rf(a4,c4,c5)*Rf(a4,c5,c6)*Rf(a5,c6,c7)+
	         i*f(a5,a4,c6)*Rf(a4,c4,c5)*Rf(c6,c5,c7)
L11	Id,Adiso,Rf(a4~,c4~,c5~)*Rf(a4~,c5~,c6~)=CAS2(Rep)*D(c4,c6)
	ENDDO
C here write RR=1/2{R,R}+1/2[R,R] anticommut.piece --> 0
	DO ii=1,3
L16	Id,Asymm,f,1,2,3
L17	Id,Adiso,f(a4~,a5~,a6~)*Rf(a5~,c4~,c5~)*Rf(a6~,c5~,c6~)=
	         f(a4,a5,a6)*i*f(a5,a6,c5)*Rf(c5,c4,c6)*1/2
L17	Al,Adiso,f(a4~,a6~,a5~)*Rf(a5~,c4~,c5~)*Rf(a6~,c5~,c6~)=
	         - f(a4,a5,a6)*i*f(a5,a6,c5)*Rf(c5,c4,c6)*1/2
L17	Al,Adiso,f(a5~,a6~,a4~)*Rf(a5~,c4~,c5~)*Rf(a6~,c5~,c6~)=
	         f(a4,a5,a6)*i*f(a5,a6,c5)*Rf(c5,c4,c6)*1/2
L17	Al,Adiso,f(a6~,a5~,a4~)*Rf(a5~,c4~,c5~)*Rf(a6~,c5~,c6~)=
	         - f(a4,a5,a6)*i*f(a5,a6,c5)*Rf(c5,c4,c6)*1/2
L17	Al,Adiso,f(a5~,a4~,a6~)*Rf(a5~,c4~,c5~)*Rf(a6~,c5~,c6~)=
	         - f(a4,a5,a6)*i*f(a5,a6,c5)*Rf(c5,c4,c6)*1/2
L17	Al,Adiso,f(a6~,a4~,a5~)*Rf(a5~,c4~,c5~)*Rf(a6~,c5~,c6~)=
	         f(a4,a5,a6)*i*f(a5,a6,c5)*Rf(c5,c4,c6)*1/2
L18	Id,Adiso,Rf(a4~,c4~,c5~)*Rf(a4~,c5~,c6~)=CAS2(Rep)*D(c4,c6)
	ENDDO
L25	Id,Asymm,f,1,2,3
	C
	DO ii=1,3
L26	Id,Adiso,Rf(a4~,c4~,c5~)*Rf(a5~,c5~,c4~)=Tg(Rep)*D(a4,a5)
L26	Al,Adiso,Rf(a4~,c4~,c5~)*Rf(a4~,c5~,c6~)=CAS2(Rep)*D(c4,c6)
L27	Id,Adiso,f(a4~,a6~,a7~)*f(a5~,a6~,a7~)=CAS2(Gr)*D(a4,a5)
L27	Al,Adiso,f(a6~,a4~,a7~)*f(a5~,a6~,a7~)= - CAS2(Gr)*D(a4,a5)
L27	Al,Adiso,f(a6~,a7~,a4~)*f(a5~,a6~,a7~)=CAS2(Gr)*D(a4,a5)
L27	Al,Adiso,f(a4~,a6~,a7~)*f(a6~,a5~,a7~)= - CAS2(Gr)*D(a4,a5)
L27	Al,Adiso,f(a6~,a4~,a7~)*f(a6~,a5~,a7~)=CAS2(Gr)*D(a4,a5)
L27	Al,Adiso,f(a6~,a7~,a4~)*f(a6~,a5~,a7~)= - CAS2(Gr)*D(a4,a5)
L27	Al,Adiso,f(a4~,a6~,a7~)*f(a6~,a7~,a5~)=CAS2(Gr)*D(a4,a5)
L27	Al,Adiso,f(a6~,a4~,a7~)*f(a6~,a7~,a5~)= - CAS2(Gr)*D(a4,a5)
L27	Al,Adiso,f(a6~,a7~,a4~)*f(a6~,a7~,a5~)=CAS2(Gr)*D(a4,a5)
	ENDDO
	C
C trace of 3 diff. Rs ( for fund. repr. of SU(N) put Tg(Rep)=1/2)
L32	Id,Adiso,Rf(a1~,a4~,a5~)*Rf(a2~,a5~,a6~)*Rf(a3~,a6~,a4~)=
	         1/2*Tg(Rep)*i*f(a1,a2,a3)
	        +1/2*Tg(Rep)*dsym(a1,a2,a3)
L33	Id,Symme,dsym,1,2,3
	C
L34	Id,Asymm,f,1,2,3
L35	Id,Rf(a4~,c4~,c4~)=0
L36	Id,Commu,CAS2,Tg
L37	Id,D(a1,a2)=KRO(a1,a2)
	*yep
	C
L 1	Id,Adiso,CAS2(Rep)*Tg(Rep)*KRO(a1,a2)=GRF#1
L 1	Al,Adiso,CAS2(Gr)*Tg(Rep)*KRO(a1,a2)=GRF#2
L 1	Al,Adiso,CAS2(Rep,Gr)*Tg(Rep)*KRO(a1,a2)=GRF#3
L 2	IF NOT flag3
L 2	Al,flagA*GRF#2=-2*GRF#3+2*GRF#1
L 2	Al,CAS2(Rep)*g^3=(CAS2(Rep,Gr)+1/2*CAS2(Gr))*g^3
	ENDIF
L 4	Id,flagA=1
C GRF#3 involves C2(R)-C2(G).1/2 >CAS2(Rep,Gr)
L 5	Id,g^4*GRF#1=GRF1
L 5	Al,g^4*GRF#2=GRF2
L 5	Al,g^4*GRF#3=GRF3
	*yep
	C
> PFF{,,}
	C
L 1	Id,Commu,nom
L 2	Id,nom(K,m~)*nom(U1,m1~)*nom(U3,m2~)=C0(m,m1,m2,P1,P2)
L 2	Al,nom(K,m~)*nom(U1,m1~)=B0(m,m1,P1)
L 2	Al,nom(K,m~)*nom(U2,m1~)=B0(m,m1,P2)
L 2	Al,nom(K,m~)*nom(U3,m1~)=B0(m,m1,-P3)
L 2	Al,nom(K,m~)=A0(m)
L 3	Id,Gammas,"C

C remains U1=K+P1,U2=K+P2,U3=K+P1+P2,P3=-P1-P2
C order of vectors:P1,P2,P3,K,U1,U2,U3
L11	Id,U1(mu~)=K(mu)+P1(mu)
L11	Al,Dotpr,U1(mu~)=K(mu)+P1(mu)
L11	Al,Funct,G,Epf,U1(mu~)=K(mu)+P1(mu)
L12	Id,U2(mu~)=K(mu)+P2(mu)
L12	Al,Dotpr,U2(mu~)=K(mu)+P2(mu)
L12	Al,Funct,G,Epf,U2(mu~)=K(mu)+P2(mu)
L13	Id,U3(mu~)=K(mu)+P1(mu)+P2(mu)
L13	Al,Dotpr,U3(mu~)=K(mu)+P1(mu)+P2(mu)
L13	Al,Funct,G,Epf,U3(mu~)=K(mu)+P1(mu)+P2(mu)
L14	Id,P3(mu~)=-P1(mu)-P2(mu)
L14	Al,Dotpr,P3(mu~)=-P1(mu)-P2(mu)
L14	Al,Funct,G,Epf,P3(mu~)=-P1(mu)-P2(mu)
L15	Id,flagA=1
L16	Id,Adiso,nom(U1,m~)*nom(U3,m1~)=B0(m,m1,P2)*flagA
L17	IF flagA=1
L17	Al,K(mu~)=K(mu)-P1(mu)
L17	Al,Dotpr,K(mu~)=K(mu)-P1(mu)
L17	Al,Funct,G,Epf,K(mu~)=K(mu)-P1(mu)
	ENDIF
> INTEG{}
L18	Id,Count,1,G,1,dsym,-1
L19	Id,Count,3,"F,G,K,1

L20	IF NOT G(s4,"t,"4,MU4,K,NU4,K,LA4,K)
L21	Id,Addfa,0
	ENDIF

> P output
	*yep
 
vertg = + G(s4,"t,"4,MU4,K,NU4,K,LA4,K)*f(a1,a2,a3)*Tg(Rep)*C0(m,m,m,P1,P2)
  * ( 1/2*i*g^3 ) + 0.
 
L 1	Id,All,K,N,FK
> P output

L 2	Id,C0(m~,m1~,m2~,P1~,P2~)*FK(AL~,BE~,GA~)=
	      P1(AL)*P1(BE)*P1(GA)*C31(m,m1,m2,P1,P2)
	     +P2(AL)*P2(BE)*P2(GA)*C32(m,m1,m2,P1,P2)     
	     +P1(AL)*P1(BE)*P2(GA)*C33(m,m1,m2,P1,P2)
	     +P1(AL)*P2(BE)*P1(GA)*C33(m,m1,m2,P1,P2)
	     +P2(AL)*P1(BE)*P1(GA)*C33(m,m1,m2,P1,P2)
	     +P1(AL)*P2(BE)*P2(GA)*C34(m,m1,m2,P1,P2)
	     +P2(AL)*P1(BE)*P2(GA)*C34(m,m1,m2,P1,P2)
	     +P2(AL)*P2(BE)*P1(GA)*C34(m,m1,m2,P1,P2)
	     +P1(AL)*D(BE,GA)*C35(m,m1,m2,P1,P2)
	     +P1(BE)*D(AL,GA)*C35(m,m1,m2,P1,P2)
	     +P1(GA)*D(AL,BE)*C35(m,m1,m2,P1,P2)
	     +P2(AL)*D(BE,GA)*C36(m,m1,m2,P1,P2)
	     +P2(BE)*D(AL,GA)*C36(m,m1,m2,P1,P2)
	     +P2(GA)*D(AL,BE)*C36(m,m1,m2,P1,P2)
L 2	Al,C0(m~,m1~,m2~,P1~,P2~)*FK(AL~,BE~)=
	      P1(AL)*P1(BE)*C21(m,m1,m2,P1,P2)
	     +P2(AL)*P2(BE)*C22(m,m1,m2,P1,P2)
	     +P1(AL)*P2(BE)*C23(m,m1,m2,P1,P2)
	     +P2(AL)*P1(BE)*C23(m,m1,m2,P1,P2)
	     +D(AL,BE)*C24(m,m1,m2,P1,P2)
L 2	Al,C0(m~,m1~,m2~,P1~,P2~)*FK(AL~)=
	      P1(AL)*C11(m,m1,m2,P1,P2)
	     +P2(AL)*C12(m,m1,m2,P1,P2)
L 3	Id,B0(m~,m1~,P1~)*FK(AL~,BE~)=
	      P1(AL)*P1(BE)*B21(m,m1,P1)
	     +D(AL,BE)*B22(m,m1,P1)
L 3	Al,B0(m~,m1~,P1~)*FK(AL~)=P1(AL)*B1(m,m1,P1)
L 4	Id,A0(m~)*FK(AL~)=0
L 5	Id,Even,B0,3,B1,3,B21,3,B22,3
	C

L 6	IF NOT C35(m,m,m,P1,P2)
L 7	Id,Addfa,0
	ENDIF

L 8	IF G(s4,"t,"4,MU4,Naa~,NU4,Naa~,LA4,P1)
L 9	Id,Addfa,Xx1
	ENDIF

L10	IF G(s4,"t,"4,MU4,Naa~,NU4,P1,LA4,Naa~)
L11	Id,Addfa,Xx2
	ENDIF
L12	IF G(s4,"t,"4,MU4,P1,NU4,Nac~,LA4,Nac~)
L13	Id,Addfa,Xx3
	ENDIF
L14	Id,Numer,Xx3,0
> P output
	*yep
 
vertg = 
 + G(s4,"t,"4,MU4,Nac,NU4,Nac,LA4,P1)*f(a1,a2,a3)*Tg(Rep)*C35(m,m,m,P1,P2)
  * ( 1/2*i*g^3*Xx1 )
 
 + G(s4,"t,"4,MU4,Nac,NU4,P1,LA4,Nac)*f(a1,a2,a3)*Tg(Rep)*C35(m,m,m,P1,P2)
  * ( 1/2*i*g^3*Xx2 ) + 0.
 
> P output
C *end
C Id,Xx2=0
L 1	Id,G(s4,"t,"4,MU4,Naa~,NU4,Naa~,LA4,P1)=D(Naa,Naa)
L 1	Al,G(s4,"t,"4,MU4,Naa~,NU4,P1,LA4,Naa~)=D(Naa,Naa)
	*end
 
vertg = + f(a1,a2,a3)*Tg(Rep)*C35(m,m,m,P1,P2)
  * ( 1/2*i*N*g^3*Xx1 + 1/2*i*N*g^3*Xx2 ) + 0.
 

End run. Time 22 sec.
