i i,j;
i ab1,ab2;
i mu,nu,al,be;
i b1,b2,b3,b4;
i a1,a2,a3,a4,a5,a6;
i aa1,aa2,aa3,aa4,aa5,aa6,aa7;
i bb1,bb2,bb3,bb4,bb5,bb6,bb7;
*
v pp,qp,Qp,Qm,k,kb,BQ;
*
s n;
s zi;
s I0;
s PZ;
s XX;
s PZs;
s Com;
s gexp;
s GT,T;
s selv;
s Xpi;
s n1,n2;
s [2-3];
s [5-6];
s [2+3];
s [5+6];
s eq,egt;
s Qu,Qd,srt;
s denk,denkb;
s one,aux,int;
s tau,x1,x2,z;
s rl1,rl2,rl3;
s sw,swC,Gw,GwC;
s gramd,restore;
s l1,l2,l,lp,lv;
s vv,vvs,vvc,vvq;
s sel0,selu,seluu;
s O,Os,Oc,Oq,Op,Om;
s fws,fwt,fwsC,fwtC;
s um,dm,ums,dms,QmS,QpS;
s omx1,x12,x2z,tx;
s Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2;
s Pws,Pwt,PwsC,PwtC,PwsM2,PwtM2;
s [fws-fwt],[fwsC-fwtC],[Sigmas-Sigmat];
s Sigmas,SigmasC,Sigmat,SigmatC,SigmasM2;
s [1+X],[x2-z],[tau+X+x1],[x1-x2],[1-x1]:
s Xps,Xpis,yc,ys,omx1s,omx1c,omx1q,x12s,x12c,x2zq,x2zc,x2zs,txs;
s y,s,rs,sp,t,tp,u,up,kp,km,zp,zm,tM,uM,X,zX,Xp,Xph,tuM,tdM,Y,zXi;
s cst,cqt,vet,vdt,vut,aet,adt,aut,vets,vdts,vuts,aets,adts,auts,sst;
*
f sl,sq,slC,sqC;
f sum1_gg,sum1_zz,sum2;
f wsub_gg,wsub_gz,wsub_zz;
*
cf es;
cf fun;
cf Rep;
cf redf;
cf kf,qpf;
cf Mods,Re,Im;
cf ae,au,ad;
cf tv1g,tv1gC;
cf tv1z,tv1zC;
cf tv0,tv1,tv0C,tv1C;
cf ve,vu,vd,ss,cs,csi;
cf vvvg,vvvgC,wg,wgC,wgM2;
cf vvvz,vvvzC,wz,wzC,wzM2;
cf I_uu_21,I_uu_22,I_uu_23,I_uu_24;
cf I_dd_21,I_dd_22,I_dd_23,I_dd_24;
cf I_nn_21,I_nn_22,I_nn_23,I_nn_24;
cf rs1wg,rs2wg,rt1wg,rt2wg,rwwg;
cf rs1wz,rs2wz,rt1wz,rt2wz,rwwz,rwwgz;
cf I_un_21,I_un_22,I_un_23,I_un_24,I_un_25;
cf I_dn_21,I_dn_22,I_dn_23,I_dn_24,I_dn_25;
cf I_u_11,I_u_12,I_d_11,I_d_12,I_n_11,I_n_12;
cf C00,B00t,B00s,B00g,C00C,B00tC,B00sC,B00gC;
*
.global
*
*----
*
#procedure wtensor(compo)
*
* First W_{mu,nu} is constructed, gg or gz or zz
*
g w_'compo'= wsub_'compo'(mu,nu)*int;
*
id wsub_gg(mu?,nu?)= sum_(l1,1,4,sum1_gg(l1,mu,nu));
*
*id wsub_zz(mu?,nu?)= sum_(l1,5,9,sum1_zz(l1,mu,nu));
id wsub_zz(mu?,nu?)= sum2(7,7,mu,nu)+sum2(7,9,mu,nu)+
                     sum2(9,7,mu,nu)+sum2(9,9,mu,nu);
*
id wsub_gz(mu?,nu?)= sum_(l1,1,4,sum1_zz(l1,mu,nu))+
                     sum_(l1,5,9,sum1_gg(l1,mu,nu));
*
id sum1_gg(l?,mu?,nu?)= sum_(l2,1,4,sum2(l,l2,mu,nu));
id sum1_zz(l?,mu?,nu?)= sum_(l2,5,9,sum2(l,l2,mu,nu));
*
id sum2(l?,lp?,mu?,nu?)= -sqC(i,l,nu)*(-i_*g_(i,k)+um*gi_(i))*
                          sq(i,lp,mu)*(-i_*g_(i,kb)-dm*gi_(i))*
                          slC(j,l,nu)*g_(j,pp)*sl(j,lp,mu)*
                          (g_(j,pp)-g_(j,Qp));
*
#call diagrams{dummy}
*
trace4,i;
trace4,j;
contract;
*
.store
*
#call kinematics{dummy}
#call lin{'compo'}
#call integration{'compo'}
*
#endprocedure
*
*---------------------------------------------------------------
*
#procedure diagrams(dummy)
*
* The diagrams are constructed
*
id sl(i?,1,mu?)= i_/8*Gt*
                (tv0(mu,a1,a5)+Gs*(tv1g(mu,a1,a5)-Sigmas*tv0(mu,a1,a5)))*
                 es(T)*g_(i,a1)*g6_(i)*vvs/Pws/Pwt/s^2;
id slC(i?,1,nu?)= i_/8*GtC*
                (tv0C(nu,b1,a6)+GsC*(tv1gC(nu,b1,a6)-SigmasC*tv0C(nu,b1,a6)))*
                 es(T)*g_(i,b1)*g6_(i)*vvs/PwsC/PwtC/s^2;
*
id sl(i?,5,mu?)= i_/32*Gt*
          (tv0(mu,aa1,aa5)+Gs*(tv1z(mu,aa1,aa5)-Sigmas*tv0(mu,aa1,aa5)))*
           GT/PZ*g_(i,aa1)*g6_(i)*vvs/Pws/Pwt/s^3;
id slC(i?,5,nu?)= i_/32*GtC*
          (tv0C(nu,bb1,aa6)+GsC*(tv1zC(nu,bb1,aa6)-SigmasC*tv0C(nu,bb1,aa6)))*
           GT/PZ*g_(i,bb1)*g6_(i)*vvs/PwsC/PwtC/s^3;
*
al sl(i?,2,mu?)= i_/(8*y*s)*Gs*es(T)*
                 g_(i,mu)*(g_(i,pp)+g_(i,Qm))*g_(i,a2)*g6_(i)*vv/Pws/s;
al slC(i?,2,nu?)= i_/(8*y*s)*GsC*es(T)*
                 g_(i,b2)*g6_(i)*(g_(i,pp)+g_(i,Qm))*g_(i,nu)*vv/PwsC/s;
*
al sl(i?,6,mu?)= -i_/(8*y*s)*Gs*GT*csi(T)/16/PZ*
                 g_(i,mu)*(ve(T)*g6_(i)+ae(T)*g7_(i))*
                 (g_(i,pp)+g_(i,Qm))*g_(i,aa2)*g6_(i)*vv/Pws/s^2;
al slC(i?,6,nu?)= -i_/(8*y*s)*GsC*GT*csi(T)/16/PZ*
                  g_(i,bb2)*g6_(i)*(g_(i,pp)+g_(i,Qm))*
                  g_(i,nu)*(ve(T)*g6_(i)+ae(T)*g7_(i))*vv/PwsC/s^2;
*
al sl(i?,3,mu?)= i_/8*Qu*Gt*es(T)*
                 g_(i,a3)*g6_(i)*vv/Pwt/s/denk;
al slC(i?,3,nu?)= i_/8*Qu*GtC*es(T)*
                 g_(i,b3)*g6_(i)*vv/PwtC/s/denk;
*
al sl(i?,7,mu?)= -i_/8*Gt*GT*csi(T)/16/PZ*
                 g_(i,aa3)*g6_(i)*vv/Pwt/s^2/denk;
al slC(i?,7,nu?)= -i_/8*GtC*GT*csi(T)/16/PZ*
                  g_(i,bb3)*g6_(i)*vv/PwtC/s^2/denk;
*
al sl(i?,4,mu?)= -i_/8*Qd*Gt*es(T)*
                 g_(i,a4)*g6_(i)*vv/Pwt/s/denkb;
al slC(i?,4,nu?)= -i_/8*Qd*GtC*es(T)*
                  g_(i,b4)*g6_(i)*vv/PwtC/s/denkb;
*
al sl(i?,8,mu?)= -i_/8*Gt*GT*csi(T)/16/PZ*
                 g_(i,aa4)*g6_(i)*vv/Pwt/s^2/denkb;
al slC(i?,8,nu?)= -i_/8*GtC*GT*csi(T)/16/PZ*
                  g_(i,bb4)*g6_(i)*vv/PwtC/s^2/denkb;
*
al sl(i?,9,mu?)= i_/32/km*Gs*GT*csi(T)/PZ*
                 g_(i,aa7)*g6_(i)*(g_(i,Qm)-g_(i,qp))*
                 g_(i,mu)*g6_(i)*vv/Pws/s^2/4;
al slC(i?,9,nu?)= i_/32/km*GsC*GT*csi(T)/PZ*
                  g_(i,nu)*g6_(i)*(g_(i,Qm)-g_(i,qp))*
                  g_(i,bb7)*g6_(i)*vv/PwsC/s^2/4;
*
id sq(i?,1,mu?)= g_(i,a5)*g6_(i);
id sqC(i?,1,nu?)= -g_(i,a6)*g6_(i);
*
al sq(i?,5,mu?)= g_(i,aa5)*g6_(i);
al sqC(i?,5,nu?)= -g_(i,aa6)*g6_(i);
*
al sq(i?,2,mu?)= g_(i,a2)*g6_(i);
al sqC(i?,2,nu?)= -g_(i,b2)*g6_(i);
*
al sq(i?,6,mu?)= g_(i,aa2)*g6_(i);
al sqC(i?,6,nu?)= -g_(i,bb2)*g6_(i);
*
al sq(i?,3,mu?)= g_(i,mu)*(g_(i,Qm)-g_(i,k)-i_*um*gi_(i))*g_(i,a3)*g6_(i);
al sqC(i?,3,nu?)= g_(i,b3)*g6_(i)*(-(g_(i,Qm)-g_(i,k))+i_*um*gi_(i))*g_(i,nu);
*
al sq(i?,7,mu?)= g_(i,mu)*(vu(T)*g6_(i)+au(T)*g7_(i))*
                 (g_(i,Qm)-g_(i,k)-i_*um*gi_(i))*g_(i,aa3)*g6_(i);
al sqC(i?,7,nu?)= g_(i,bb3)*g6_(i)*(-(g_(i,Qm)-g_(i,k))+i_*um*gi_(i))*
                  g_(i,nu)*(vu(T)*g6_(i)+au(T)*g7_(i));
*
al sq(i?,4,mu?)= g_(i,a4)*g6_(i)*(g_(i,Qm)-g_(i,kb)+i_*dm*gi_(i))*g_(i,mu);
al sqC(i?,4,nu?)= g_(i,nu)*(-(g_(i,Qm)-g_(i,kb))-i_*dm*gi_(i))*g_(i,b4)*g6_(i);
*
al sq(i?,8,mu?)= g_(i,aa4)*g6_(i)*(g_(i,Qm)-g_(i,kb)+i_*dm*gi_(i))*
                 g_(i,mu)*(vd(T)*g6_(i)+ad(T)*g7_(i)) ;
al sqC(i?,8,nu?)= g_(i,nu)*(vd(T)*g6_(i)+ad(T)*g7_(i))*
                  (-(g_(i,Qm)-g_(i,kb))-i_*dm*gi_(i))*g_(i,bb4)*g6_(i);
*
al sq(i?,9,mu?)= g_(i,aa7)*g6_(i);
al sqC(i?,9,nu?)= -g_(i,bb7)*g6_(i);
*
#endprocedure
*
.store
*
#procedure kinematics(dummy)
*
* Mostly kinematics
*
g ww_'compo'= w_'compo';
*nwrite statistics;
*
id Qu= 2/3;
al Qd= -1/3;
id srt^-2= 1/2;
*
id pp.pp= 0;
al qp.qp= 0;
al k.k= -ums;
al kb.kb= -dms;
al Qm.Qm= QmS;
id um^2= ums;
al dm^2= dms;
*
*---------------------------------------------------------------------
id Gs= gexp*Gs;
al GsC= gexp*GsC;
id Gt= gexp*Gt;
al GtC= gexp*GtC;
*---------------------------------------------------------------------
*
id Gs*GsC= GsM2;
al Gt*GtC= GtM2;
*
id Pws^-1*PwsC^-1= 1/PwsM2;
id Pwt^-1*PwtC^-1= 1/PwtM2;
id Pws^-1= PwsC/PwsM2;
id PwsC^-1= Pws/PwsM2;
id Pwt^-1= PwtC/PwtM2;
id PwtC^-1= Pwt/PwtM2;
*
* Triple-vertex is substituted
*
id tv0(mu?,al?,be?)= d_(mu,al)*(Qm(be)-Qp(be))+
                     d_(al,be)*(2*Qp(mu)+Qm(mu))-
                     d_(mu,be)*(2*Qm(al)+Qp(al));
al tv0C(mu?,al?,be?)= d_(mu,al)*(Qm(be)-Qp(be))+
                     d_(al,be)*(2*Qp(mu)+Qm(mu))-
                     d_(mu,be)*(2*Qm(al)+Qp(al));
id tv1g(mu?,al?,be?)= sum_(lv,1,10,(vvvg(lv,mu,al,be)));
id tv1gC(mu?,al?,be?)= sum_(lv,1,10,(vvvgC(lv,mu,al,be)));
id tv1z(mu?,al?,be?)= sum_(lv,1,10,(vvvz(lv,mu,al,be)));
id tv1zC(mu?,al?,be?)= sum_(lv,1,10,(vvvzC(lv,mu,al,be)));
*
* Form factors are introduced
*
id vvvg(1,mu?,al?,be?)= 4*wg(1)/rs*(d_(al,be)*Qp(mu)+
                       d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al)));
id vvvg(2,mu?,al?,be?)= 2*wg(2)/rs*(d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al)));
id vvvg(3,mu?,al?,be?)= 2*wg(3)/rs*(d_(mu,al)*Qm(be)+d_(mu,be)*(Qm(al)+Qp(al)));
id vvvg(4,mu?,al?,be?)= 2*wg(4)/rs/(X*y*s)*Qm(be)*(Qp(al)+Qm(al))*Qp(mu);
id vvvg(5,mu?,al?,be?)= wg(5)/rs*e_(Qm,mu,al,be);
id vvvg(6,mu?,al?,be?)= wg(6)/rs*(e_(Qm,mu,al,be)+2*e_(Qp,mu,al,be));
id vvvg(7,mu?,al?,be?)= wg(7)/rs/(X*y*s)*e_(Qp,Qm,al,be)*(Qp(mu)+Qm(mu));
id vvvg(8,mu?,al?,be?)= 2*wg(8)/rs/(X*y*s)*Qm(be)*Qm(mu)*(Qp(al)+Qm(al));
id vvvg(9,mu?,al?,be?)= 2*wg(9)/rs*d_(al,be)*Qm(mu);
id vvvg(10,mu?,al?,be?)= -wg(10)/rs/(X*y*s)*e_(Qp,Qm,al,be)*Qp(mu);
*
id vvvgC(1,mu?,al?,be?)= 4*wgC(1)/rs*(d_(al,be)*Qp(mu)+
                        d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al)));
id vvvgC(2,mu?,al?,be?)= 2*wgC(2)/rs*(d_(mu,al)*Qm(be)-d_(mu,be)*
                        (Qm(al)+Qp(al)));
id vvvgC(3,mu?,al?,be?)= 2*wgC(3)/rs*(d_(mu,al)*Qm(be)+d_(mu,be)*
                        (Qm(al)+Qp(al)));
id vvvgC(4,mu?,al?,be?)= 2*wgC(4)/rs/(X*y*s)*Qm(be)*(Qp(al)+Qm(al))*Qp(mu);
id vvvgC(5,mu?,al?,be?)= -wgC(5)/rs*e_(Qm,mu,al,be);
id vvvgC(6,mu?,al?,be?)= -wgC(6)/rs*(e_(Qm,mu,al,be)+2*e_(Qp,mu,al,be));
id vvvgC(7,mu?,al?,be?)= -wgC(7)/rs/(X*y*s)*e_(Qp,Qm,al,be)*(Qp(mu)+Qm(mu));
id vvvgC(8,mu?,al?,be?)= 2*wgC(8)/rs/(X*y*s)*Qm(be)*Qm(mu)*(Qp(al)+Qm(al));
id vvvgC(9,mu?,al?,be?)= 2*wgC(9)/rs*d_(al,be)*Qm(mu);
id vvvgC(10,mu?,al?,be?)= wgC(10)/rs/(X*y*s)*e_(Qp,Qm,al,be)*Qp(mu);
*
id vvvz(1,mu?,al?,be?)= 4*wz(1)/rs*(d_(al,be)*Qp(mu)+
                       d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al)));
id vvvz(2,mu?,al?,be?)= 2*wz(2)/rs*(d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al)));
id vvvz(3,mu?,al?,be?)= 2*wz(3)/rs*(d_(mu,al)*Qm(be)+d_(mu,be)*(Qm(al)+Qp(al)));
id vvvz(4,mu?,al?,be?)= 2*wz(4)/rs/(X*y*s)*Qm(be)*(Qp(al)+Qm(al))*Qp(mu);
id vvvz(5,mu?,al?,be?)= wz(5)/rs*e_(Qm,mu,al,be);
id vvvz(6,mu?,al?,be?)= wz(6)/rs*(e_(Qm,mu,al,be)+2*e_(Qp,mu,al,be));
id vvvz(7,mu?,al?,be?)= wz(7)/rs/(X*y*s)*e_(Qp,Qm,al,be)*(Qp(mu)+Qm(mu));
id vvvz(8,mu?,al?,be?)= 2*wz(8)/rs/(X*y*s)*Qm(be)*Qm(mu)*(Qp(al)+Qm(al));
id vvvz(9,mu?,al?,be?)= 2*wz(9)/rs*d_(al,be)*Qm(mu);
id vvvz(10,mu?,al?,be?)= -wz(10)/rs/(X*y*s)*e_(Qp,Qm,al,be)*Qp(mu);
*
id vvvzC(1,mu?,al?,be?)= 4*wzC(1)/rs*(d_(al,be)*Qp(mu)+
                        d_(mu,al)*Qm(be)-d_(mu,be)*(Qm(al)+Qp(al)));
id vvvzC(2,mu?,al?,be?)= 2*wzC(2)/rs*(d_(mu,al)*Qm(be)-d_(mu,be)*
                        (Qm(al)+Qp(al)));
id vvvzC(3,mu?,al?,be?)= 2*wzC(3)/rs*(d_(mu,al)*Qm(be)+d_(mu,be)*
                        (Qm(al)+Qp(al)));
id vvvzC(4,mu?,al?,be?)= 2*wzC(4)/rs/(X*y*s)*Qm(be)*(Qp(al)+Qm(al))*Qp(mu);
id vvvzC(5,mu?,al?,be?)= -wzC(5)/rs*e_(Qm,mu,al,be);
id vvvzC(6,mu?,al?,be?)= -wzC(6)/rs*(e_(Qm,mu,al,be)+2*e_(Qp,mu,al,be));
id vvvzC(7,mu?,al?,be?)= -wzC(7)/rs/(X*y*s)*e_(Qp,Qm,al,be)*(Qp(mu)+Qm(mu));
id vvvzC(8,mu?,al?,be?)= 2*wzC(8)/rs/(X*y*s)*Qm(be)*Qm(mu)*(Qp(al)+Qm(al));
id vvvzC(9,mu?,al?,be?)= 2*wzC(9)/rs*d_(al,be)*Qm(mu);
id vvvzC(10,mu?,al?,be?)= wzC(10)/rs/(X*y*s)*e_(Qp,Qm,al,be)*Qp(mu);
*
contract;
*
id Sigmas= 0;
al SigmasC= 0;
*id Sigmas*SigmasC= SigmasM2;
*
id BQ(mu?)= pp(mu)-qp(mu)+Qm(mu);
id Qp(mu?)= pp(mu)-qp(mu);
*
id e_(a1?,Qp,a2?,a3?)= e_(a1,pp,a2,a3)-e_(a1,qp,a2,a3);
id e_(a1?,a2?,Qp,a3?)= e_(a1,a2,pp,a3)-e_(a1,a2,qp,a3);
id e_(a1?,a2?,a3?,Qp)= e_(a1,a2,a3,pp)-e_(a1,a2,a3,qp);
*
id Qp.Qm= pp.Qm-qp.Qm;
id kb.Qp= kb.pp-kb.qp;
id k.Qp= k.pp-k.qp;
id pp.Qp= -pp.qp;
id Qp.Qp= -2*pp.qp;
*
.store
*
#endprocedure
*
#procedure lin(compo)
*
* Linear independent invariants are introduced
*
g www_'compo'= ww_'compo';
*nwrite statistics;
*
contract;
*
id k.k= -ums;
al kb.kb= -dms;
*
id denk^-1*denkb^-1= -1/2/(Qp.Qm)*(1/denk+1/denkb);
*
id pp.pp= 0;
al qp.qp= 0;
al k.k= -ums;
al kb.kb= -dms;
al Qm.Qm= QmS;
*
id denk^-1= -1/tM;
id denkb^-1= 0;
*
id e_(kb,a1?,a2?,a3?)= e_(pp,a1,a2,a3)-e_(qp,a1,a2,a3)+e_(Qm,a1,a2,a3)-
                       e_(k,a1,a2,a3);
id kb(mu?)= pp(mu)-qp(mu)+Qm(mu)-k(mu);
id QpS= -2*pp.qp;
id Qp.Qm= pp.Qm-qp.Qm;
id kb.Qp= kb.pp-kb.qp;
id k.Qp= k.pp-k.qp;
id pp.Qp= -pp.qp;
id Qp.Qp= -2*pp.qp;
id Qp.qp= pp.qp;
*
id pp.Qm= -1/2*(QmS+y*s);
al k.kb= 1/2*(-sp+ums+dms);
al Qm.k= 1/2*(QmS+t-ums);
al pp.kb= 1/2*(tp-dms);
al Qm.kb= 1/2*(QmS+u-dms);
al pp.k= 1/2*(up-ums);
al pp.qp= 1/2*kp;
al Qm.qp= 1/2*(QmS+km);
al kb.qp= 1/2*(dms-zp);
al k.qp= 1/2*(ums-zm);
*
id t= tau*y*s;
al sp= x2*y*s;
al km= z*y*s;
al km^-1= 1/(z*y*s);
al zm= -(x2-x1)*y*s;
id QmS= X*y*s;
*
ratio,tM,t,ums;
ratio,uM,u,dms;
*
id ums*tM^-2= aux;
id ums= 0;
id aux= ums*tM^-2;
id dms= 0;
*
.store
*
g wwww_'compo'= www_'compo';
*nwrite statistics;
*
id kp= y*s * (  - 1 - X + x2 - z );
al u= y*s * (  - 1 - tau - 2*X - z )
       + ums + dms;
al up= y*s * (  - tau - X - x1 )
       + 2*ums + dms;
al tp= y*s * ( tau + X + x1 - x2 + z )
       - ums;
al zp= y*s * ( 1 - x1 )
       + ums + dms;
*
id Qp.Qm^-1= -2/(zX*y*s);
ratio,z,zX,aux;
id aux= 1+2*X;
id aux^-1= 1/2/Xph;
ratio,X,Xph,aux;
id aux= 1/2;
al aux^-1= 2;
*
id k(mu?)= kf(mu);
id qp(mu?)= qpf(mu);
id e_(k,mu?,nu?,al?)= kf(ab1)*e_(ab1,mu,nu,al);
id e_(qp,mu?,nu?,al?)= qpf(ab2)*e_(ab2,mu,nu,al);
*
id ums*tM^-2= aux;
id ums= 0;
id aux= ums*tM^-2;
id dms= 0;
*
b i_,Gs,GsC,Gt,GtC,GT,GsM2,GtM2,Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,ve,ae,vu,au,
  vd,ad,csi,gexp,d_,e_,pp,Qm,PZ,kf,qpf;
print;
.store
*
#endprocedure
*
#procedure integration(compo)
*
* Integration starts
*
g WT1_'compo'= wwww_'compo';
*nwrite statistics;
*
id kf(mu?)*kf(nu?)*int= I_uu_21*pp(mu)*pp(nu)+I_uu_22*Qm(mu)*Qm(nu)+
                        I_uu_23*(pp(mu)*Qm(nu)+pp(nu)*Qm(mu))+
                        I_uu_24*d_(mu,nu);
*
.store
*
g WT2_'compo'= WT1_'compo';
*nwrite statistics;
*
id qpf(mu?)*qpf(nu?)*int= I_nn_21*pp(mu)*pp(nu)+I_nn_22*Qm(mu)*Qm(nu)+
                          I_nn_23*(pp(mu)*Qm(nu)+pp(nu)*Qm(mu))+
                          I_nn_24*d_(mu,nu);
*
.store
*
g WT3_'compo'= WT2_'compo';
*nwrite statistics;
*
id kf(mu?)*qpf(nu?)*int= I_un_21*pp(mu)*pp(nu)+I_un_22*Qm(mu)*Qm(nu)+
                         I_un_23*pp(mu)*Qm(nu)+I_un_24*pp(nu)*Qm(mu)+
                         I_un_25*d_(mu,nu);
*
.store
*
g WT4_'compo'= WT3_'compo';
*nwrite statistics;
*
id kf(mu?)*int= I_u_11*pp(mu)+I_u_12*Qm(mu);
id qpf(mu?)*int= I_n_11*pp(mu)+I_n_12*Qm(mu);
id int= I0;
*
*
.store
*
g W_'compo'_int1= WT4_'compo';
*nwrite statistics;
*
id I_u_11 =
       + I0*ums*Xp^-2*y^-1*s^-1 * ( 2 )
       + I0*ums*Xp^-1*y^-1*s^-1 * (  - 1 )
       + I0*dms*Xp^-2*y^-1*s^-1 * ( 2 )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 2 )
       + I0*Xp^-2 * ( 2 - 2*tau - 2*x1 )
       + I0*Xp^-1 * (  - 3 + tau + 2*x1 )
       + I0 * ( 1 );
al I_n_11 =
       + I0*Xp^-2 * ( 2*x2 - 2*z )
       + I0*Xp^-1 * (  - 1 - 2*x2 + z )
       + I0 * ( 1 );
al I_u_12 =
       + I0*ums*Xp^-1*y^-1*s^-1 * (  - 1 )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 1 )
       + I0*Xp^-1 * (  - 1 + tau + x1 )
       + I0 * ( 1 );
al I_n_12 =
       + I0*Xp^-1 * (  - x2 + z )
       + I0 * ( 1 );
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
id ums*tM^-2= aux;
id ums= 0;
id aux= ums*tM^-2;
id dms= 0;
*
.store
*
g W_'compo'_int2= W_'compo'_int1;
*nwrite statistics;
*
id I_uu_21 =
       + I0*ums*Xp^-4*y^-1*s^-1 * ( 12 - 12*tau - 12*x1 )
       + I0*ums*Xp^-3*y^-1*s^-1 * (  - 24 + 12*tau + 18*x1 )
       + I0*ums*Xp^-2*y^-1*s^-1 * ( 16 - 2*tau - 6*x1 )
       + I0*ums*Xp^-1*y^-1*s^-1 * (  - 4 )
       + I0*dms*Xp^-4*y^-1*s^-1 * ( 12 - 12*tau - 12*x1 )
       + I0*dms*Xp^-3*y^-1*s^-1 * (  - 30 + 18*tau + 24*x1 )
       + I0*dms*Xp^-2*y^-1*s^-1 * ( 24 - 6*tau - 12*x1 )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 6 )
       + I0*Xp^-4 * ( 6 + 12*tau*x1 - 12*tau + 6*tau^2 - 12*x1 + 6*x1^2 )
       + I0*Xp^-3 * (  - 18 - 18*tau*x1 + 24*tau - 6*tau^2 + 30*x1 - 12*x1^2 )
       + I0*Xp^-2 * ( 19 + 6*tau*x1 - 14*tau + tau^2 - 24*x1 + 6*x1^2 )
       + I0*Xp^-1 * (  - 8 + 2*tau + 6*x1 )
       + I0 * ( 1 );
al I_nn_21 =
       + I0*Xp^-4 * (  - 12*x2*z + 6*x2^2 + 6*z^2 )
       + I0*Xp^-3 * ( 18*x2*z - 6*x2 - 12*x2^2 + 6*z - 6*z^2 )
       + I0*Xp^-2 * ( 1 - 6*x2*z + 12*x2 + 6*x2^2 - 8*z + z^2 )
       + I0*Xp^-1 * (  - 2 - 6*x2 + 2*z )
       + I0 * ( 1 );
al I_un_21 =
       + I0*ums*Xp^-4*y^-1*s^-1 * ( 6*x2 - 6*z )
       + I0*ums*Xp^-3*y^-1*s^-1 * (  - 3 - 9*x2 + 6*z )
       + I0*ums*Xp^-2*y^-1*s^-1 * ( 3 + 3*x2 - z )
       + I0*dms*Xp^-4*y^-1*s^-1 * ( 6*x2 - 6*z )
       + I0*dms*Xp^-3*y^-1*s^-1 * (  - 3 - 12*x2 + 9*z )
       + I0*dms*Xp^-2*y^-1*s^-1 * ( 6 + 6*x2 - 3*z )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 3 )
       + I0*Xp^-4 * (  - 6*tau*x2 + 6*tau*z - 6*x1*x2 + 6*x1*z + 6*x2 - 6*z )
       + I0*Xp^-3 * (  - 3 + 9*tau*x2 - 6*tau*z + 3*tau + 12*x1*x2 - 9*x1*z
          + 3*x1 - 15*x2 + 12*z )
       + I0*Xp^-2 * ( 7 - 3*tau*x2 + tau*z - 4*tau - 6*x1*x2 + 3*x1*z - 5*x1
          + 11*x2 - 7*z )
       + I0*Xp^-1 * (  - 5 + tau + 2*x1 - 2*x2 + z )
       + I0 * ( 1 );
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
id ums*tM^-2= aux;
id ums= 0;
id aux= ums*tM^-2;
id dms= 0;
*
.store
*
g W_'compo'_int3= W_'compo'_int2;
*nwrite statistics;
*
id I_uu_22 =
       + I0*ums*Xp^-2*y^-1*s^-1 * ( 2 - 2*tau - 2*x1 )
       + I0*ums*Xp^-1*y^-1*s^-1 * (  - 2 )
       + I0*dms*Xp^-2*y^-1*s^-1 * ( 2 - 2*tau - 2*x1 )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 2 )
       + I0*Xp^-2 * ( 1 + 2*tau*x1 - 2*tau + tau^2 - 2*x1 + x1^2 )
       + I0*Xp^-1 * (  - 2 + 2*tau + 2*x1 )
       + I0 * ( 1 );
al I_nn_22 =
       + I0*Xp^-2 * (  - 2*x2*z + x2^2 + z^2 )
       + I0*Xp^-1 * (  - 2*x2 + 2*z )
       + I0 * ( 1 );
al I_un_22 =
       + I0*ums*Xp^-2*y^-1*s^-1 * ( x2 - z )
       + I0*ums*Xp^-1*y^-1*s^-1 * (  - 1 )
       + I0*dms*Xp^-2*y^-1*s^-1 * ( x2 - z )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 1 )
       + I0*Xp^-2 * (  - tau*x2 + tau*z - x1*x2 + x1*z + x2 - z )
       + I0*Xp^-1 * (  - 1 + tau + x1 - x2 + z )
       + I0 * ( 1 );
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
id ums*tM^-2= aux;
id ums= 0;
id aux= ums*tM^-2;
id dms= 0;
*
.store
*
g W_'compo'_int4= W_'compo'_int3;
*nwrite statistics;
*
id I_uu_23 =
       + I0*ums*Xp^-3*y^-1*s^-1 * (  - 6 + 6*tau + 6*x1 )
       + I0*ums*Xp^-2*y^-1*s^-1 * ( 8 - 2*tau - 4*x1 )
       + I0*ums*Xp^-1*y^-1*s^-1 * (  - 3 )
       + I0*dms*Xp^-3*y^-1*s^-1 * (  - 6 + 6*tau + 6*x1 )
       + I0*dms*Xp^-2*y^-1*s^-1 * ( 10 - 4*tau - 6*x1 )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 4 )
       + I0*Xp^-3 * (  - 3 - 6*tau*x1 + 6*tau - 3*tau^2 + 6*x1 - 3*x1^2 )
       + I0*Xp^-2 * ( 7 + 4*tau*x1 - 8*tau + tau^2 - 10*x1 + 3*x1^2 )
       + I0*Xp^-1 * (  - 5 + 2*tau + 4*x1 )
       + I0 * ( 1 );
al I_nn_23 =
       + I0*Xp^-3 * ( 6*x2*z - 3*x2^2 - 3*z^2 )
       + I0*Xp^-2 * (  - 4*x2*z + 4*x2 + 3*x2^2 - 4*z + z^2 )
       + I0*Xp^-1 * (  - 1 - 4*x2 + 2*z )
       + I0 * ( 1 );
al I_un_23 =
       + I0*ums*Xp^-3*y^-1*s^-1 * (  - 3*x2 + 3*z )
       + I0*ums*Xp^-2*y^-1*s^-1 * ( 5/2 + 3/2*x2 - z )
       + I0*ums*Xp^-1*y^-1*s^-1 * (  - 1/2 )
       + I0*dms*Xp^-3*y^-1*s^-1 * (  - 3*x2 + 3*z )
       + I0*dms*Xp^-2*y^-1*s^-1 * ( 5/2 + 3*x2 - 5/2*z )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 5/2 )
       + I0*Xp^-3 * ( 3*tau*x2 - 3*tau*z + 3*x1*x2 - 3*x1*z - 3*x2 + 3*z )
       + I0*Xp^-2 * ( 5/2 - 3/2*tau*x2 + tau*z - 5/2*tau - 3*x1*x2 + 5/2*x1*z
          - 5/2*x1 + 9/2*x2 - 4*z )
       + I0*Xp^-1 * (  - 7/2 + tau + 2*x1 - x2 + z )
       + I0 * ( 1 );
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
id ums*tM^-2= aux;
id ums= 0;
id aux= ums*tM^-2;
id dms= 0;
*
.store
*
g W_'compo'_int5= W_'compo'_int4;
*nwrite statistics;
*
id I_uu_24 =
       + I0*ums*Xp^-2 * (  - 1 + tau + x1 )
       + I0*ums*Xp^-1 * ( 1 - 1/2*x1 )
       + I0*ums * (  - 1/2 )
       + I0*dms*Xp^-2 * (  - 1 + tau + x1 )
       + I0*dms*Xp^-1 * ( 3/2 - 1/2*tau - x1 )
       + I0*dms * (  - 1/2 )
       + I0*Xp^-2*y*s * (  - 1/2 - tau*x1 + tau - 1/2*tau^2 + x1 - 1/2*x1^2 )
       + I0*Xp^-1*y*s * ( 1 + 1/2*tau*x1 - tau - 3/2*x1 + 1/2*x1^2 )
       + I0*y*s * (  - 1/2 + 1/2*x1 );
al I_nn_24 =
       + I0*Xp^-2*y*s * ( x2*z - 1/2*x2^2 - 1/2*z^2 )
       + I0*Xp^-1*y*s * (  - 1/2*x2*z + 1/2*x2 + 1/2*x2^2 - 1/2*z )
       + I0*y*s * (  - 1/2*x2 );
al I_un_24 =
       + I0*ums*Xp^-3*y^-1*s^-1 * (  - 3*x2 + 3*z )
       + I0*ums*Xp^-2*y^-1*s^-1 * ( 3/2 + 5/2*x2 - z )
       + I0*ums*Xp^-1*y^-1*s^-1 * (  - 1/2 )
       + I0*dms*Xp^-3*y^-1*s^-1 * (  - 3*x2 + 3*z )
       + I0*dms*Xp^-2*y^-1*s^-1 * ( 3/2 + 3*x2 - 3/2*z )
       + I0*dms*Xp^-1*y^-1*s^-1 * (  - 3/2 )
       + I0*Xp^-3 * ( 3*tau*x2 - 3*tau*z + 3*x1*x2 - 3*x1*z - 3*x2 + 3*z )
       + I0*Xp^-2 * ( 3/2 - 5/2*tau*x2 + tau*z - 3/2*tau - 3*x1*x2 + 3/2*x1*z
          - 3/2*x1 + 11/2*x2 - 4*z )
       + I0*Xp^-1 * (  - 5/2 + tau + x1 - 2*x2 + z )
       + I0 * ( 1 );
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
id ums*tM^-2= aux;
id ums= 0;
id aux= ums*tM^-2;
id dms= 0;
*
.store
*
g W_'compo'_int6= W_'compo'_int5;
*nwrite statistics;
*
id I_un_25 =
       + I0*ums*Xp^-2 * (  - 1/2*x2 + 1/2*z )
       + I0*ums*Xp^-1 * ( 1/4 + 1/4*x2 )
       + I0*ums * ( 1/4 )
       + I0*dms*Xp^-2 * (  - 1/2*x2 + 1/2*z )
       + I0*dms*Xp^-1 * ( 1/4 + 1/2*x2 - 1/4*z )
       + I0*dms * (  - 1/4 )
       + I0*Xp^-2*y*s * ( 1/2*tau*x2 - 1/2*tau*z + 1/2*x1*x2 - 1/2*x1*z - 1/2
         *x2 + 1/2*z )
       + I0*Xp^-1*y*s * ( 1/4 - 1/4*tau*x2 - 1/4*tau - 1/2*x1*x2 + 1/4*x1*z
          - 1/4*x1 + 3/4*x2 - 1/2*z )
       + I0*y*s * (  - 1/4 );
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
id ums*tM^-2= aux;
id ums= 0;
id aux= ums*tM^-2;
id dms= 0;
*
id tM^-1= 1/(y*s)/tuM;
*
ratio,tuM,tau,aux;
id aux= ums/(y*s);
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
ratio,z,zX,aux;
id aux= 1+2*X;
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
.store
*
*save wrl$grp04:[giampiero]outputampgg_A.sav W_'compo'_int6;
*save wrl$grp04:[giampiero]outputampgz_A.sav W_'compo'_int6;
save wrl$grp04:[giampiero]outputampzz_A.sav W_'compo'_int6;
*
#endprocedure
*
#procedure totalgg(dummy)
*
* Now the full stuff
*
load wrl$grp04:[giampiero]outputampgg_A.sav;
*
g T_gg= PwsM2*PwtM2*W_gg_int6/(eq*I0);
*nwrite statistics;
*
id es(T)*es(T)= eq;
*
id wg(n?)= wg(n)*rs;
al wgC(n?)= wgC(n)*rs;
*
id wg(4)= X*wg(4);
al wg(7)= X*wg(7);
al wg(8)= X*wg(8);
al wg(10)= X*wg(10);
id wgC(4)= X*wgC(4);
al wgC(7)= X*wgC(7);
al wgC(8)= X*wgC(8);
al wgC(10)= X*wgC(10);
*
id wg(n?)*wgC(n?)= wgM2(n);
*
id GtM2*GsM2*wg(2)*wgC(1)= GtM2*GsM2*(-wgC(2)*wg(1)+
                                 2*Re(wg(2)*wgC(1)));
al GtM2*GsM2*wg(3)*wgC(1)= GtM2*GsM2*(-wgC(3)*wg(1)+
                                 2*Re(wg(3)*wgC(1)));
al GtM2*GsM2*wg(4)*wgC(1)= GtM2*GsM2*(-wgC(4)*wg(1)+
                                 2*Re(wg(4)*wgC(1)));
al GtM2*GsM2*wg(5)*wgC(1)= GtM2*GsM2*(-wgC(5)*wg(1)+
                                 2*Re(wg(5)*wgC(1)));
al GtM2*GsM2*wg(6)*wgC(1)= GtM2*GsM2*(-wgC(6)*wg(1)+
                                 2*Re(wg(6)*wgC(1)));
al GtM2*GsM2*wg(7)*wgC(1)= GtM2*GsM2*(-wgC(7)*wg(1)+
                                 2*Re(wg(7)*wgC(1)));
al GtM2*GsM2*wg(8)*wgC(1)= GtM2*GsM2*(-wgC(8)*wg(1)+
                                 2*Re(wg(8)*wgC(1)));
al GtM2*GsM2*wg(9)*wgC(1)= GtM2*GsM2*(-wgC(9)*wg(1)+
                                 2*Re(wg(9)*wgC(1)));
al GtM2*GsM2*wg(10)*wgC(1)= GtM2*GsM2*(-wgC(10)*wg(1)+
                                 2*Re(wg(10)*wgC(1)));
*
al GtM2*GsM2*wg(3)*wgC(2)= GtM2*GsM2*(-wgC(3)*wg(2)+
                                 2*Re(wg(3)*wgC(2)));
al GtM2*GsM2*wg(4)*wgC(2)= GtM2*GsM2*(-wgC(4)*wg(2)+
                                 2*Re(wg(4)*wgC(2)));
al GtM2*GsM2*wg(5)*wgC(2)= GtM2*GsM2*(-wgC(5)*wg(2)+
                                 2*Re(wg(5)*wgC(2)));
al GtM2*GsM2*wg(6)*wgC(2)= GtM2*GsM2*(-wgC(6)*wg(2)+
                                 2*Re(wg(6)*wgC(2)));
al GtM2*GsM2*wg(7)*wgC(2)= GtM2*GsM2*(-wgC(7)*wg(2)+
                                 2*Re(wg(7)*wgC(2)));
al GtM2*GsM2*wg(8)*wgC(2)= GtM2*GsM2*(-wgC(8)*wg(2)+
                                 2*Re(wg(8)*wgC(2)));
al GtM2*GsM2*wg(9)*wgC(2)= GtM2*GsM2*(-wgC(9)*wg(2)+
                                 2*Re(wg(9)*wgC(2)));
al GtM2*GsM2*wg(10)*wgC(2)= GtM2*GsM2*(-wgC(10)*wg(2)+
                                 2*Re(wg(10)*wgC(2)));
*
al GtM2*GsM2*wg(4)*wgC(3)= GtM2*GsM2*(-wgC(4)*wg(3)+
                                 2*Re(wg(4)*wgC(3)));
al GtM2*GsM2*wg(5)*wgC(3)= GtM2*GsM2*(-wgC(5)*wg(3)+
                                 2*Re(wg(5)*wgC(3)));
al GtM2*GsM2*wg(6)*wgC(3)= GtM2*GsM2*(-wgC(6)*wg(3)+
                                 2*Re(wg(6)*wgC(3)));
al GtM2*GsM2*wg(7)*wgC(3)= GtM2*GsM2*(-wgC(7)*wg(3)+
                                 2*Re(wg(7)*wgC(3)));
al GtM2*GsM2*wg(8)*wgC(3)= GtM2*GsM2*(-wgC(8)*wg(3)+
                                 2*Re(wg(8)*wgC(3)));
al GtM2*GsM2*wg(9)*wgC(3)= GtM2*GsM2*(-wgC(9)*wg(3)+
                                 2*Re(wg(9)*wgC(3)));
al GtM2*GsM2*wg(10)*wgC(3)= GtM2*GsM2*(-wgC(10)*wg(3)+
                                 2*Re(wg(10)*wgC(3)));
*
al GtM2*GsM2*wg(5)*wgC(4)= GtM2*GsM2*(-wgC(5)*wg(4)+
                                 2*Re(wg(5)*wgC(4)));
al GtM2*GsM2*wg(6)*wgC(4)= GtM2*GsM2*(-wgC(6)*wg(4)+
                                 2*Re(wg(6)*wgC(4)));
al GtM2*GsM2*wg(7)*wgC(4)= GtM2*GsM2*(-wgC(7)*wg(4)+
                                 2*Re(wg(7)*wgC(4)));
al GtM2*GsM2*wg(8)*wgC(4)= GtM2*GsM2*(-wgC(8)*wg(4)+
                                 2*Re(wg(8)*wgC(4)));
al GtM2*GsM2*wg(9)*wgC(4)= GtM2*GsM2*(-wgC(9)*wg(4)+
                                 2*Re(wg(9)*wgC(4)));
al GtM2*GsM2*wg(10)*wgC(4)= GtM2*GsM2*(-wgC(10)*wg(4)+
                                 2*Re(wg(10)*wgC(4)));
*
al GtM2*GsM2*wg(6)*wgC(5)= GtM2*GsM2*(-wgC(6)*wg(5)+
                                 2*Re(wg(6)*wgC(5)));
al GtM2*GsM2*wg(7)*wgC(5)= GtM2*GsM2*(-wgC(7)*wg(5)+
                                 2*Re(wg(7)*wgC(5)));
al GtM2*GsM2*wg(8)*wgC(5)= GtM2*GsM2*(-wgC(8)*wg(5)+
                                 2*Re(wg(8)*wgC(5)));
al GtM2*GsM2*wg(9)*wgC(5)= GtM2*GsM2*(-wgC(9)*wg(5)+
                                 2*Re(wg(9)*wgC(5)));
al GtM2*GsM2*wg(10)*wgC(5)= GtM2*GsM2*(-wgC(10)*wg(5)+
                                 2*Re(wg(10)*wgC(5)));
*
al GtM2*GsM2*wg(7)*wgC(6)= GtM2*GsM2*(-wgC(7)*wg(6)+
                                 2*Re(wg(7)*wgC(6)));
al GtM2*GsM2*wg(8)*wgC(6)= GtM2*GsM2*(-wgC(8)*wg(6)+
                                 2*Re(wg(8)*wgC(6)));
al GtM2*GsM2*wg(9)*wgC(6)= GtM2*GsM2*(-wgC(9)*wg(6)+
                                 2*Re(wg(9)*wgC(6)));
al GtM2*GsM2*wg(10)*wgC(6)= GtM2*GsM2*(-wgC(10)*wg(6)+
                                 2*Re(wg(10)*wgC(6)));
*
al GtM2*GsM2*wg(8)*wgC(7)= GtM2*GsM2*(-wgC(8)*wg(7)+
                                 2*Re(wg(8)*wgC(7)));
al GtM2*GsM2*wg(9)*wgC(7)= GtM2*GsM2*(-wgC(9)*wg(7)+
                                 2*Re(wg(9)*wgC(7)));
al GtM2*GsM2*wg(10)*wgC(7)= GtM2*GsM2*(-wgC(10)*wg(7)+
                                 2*Re(wg(10)*wgC(7)));
*
al GtM2*GsM2*wg(9)*wgC(8)= GtM2*GsM2*(-wgC(9)*wg(8)+
                                 2*Re(wg(9)*wgC(8)));
al GtM2*GsM2*wg(10)*wgC(8)= GtM2*GsM2*(-wgC(10)*wg(8)+
                                 2*Re(wg(10)*wgC(8)));
*
al GtM2*GsM2*wg(10)*wgC(9)= GtM2*GsM2*(-wgC(10)*wg(9)+
                                 2*Re(wg(10)*wgC(9)));
*
id wg(n?)*Gs*GtM2*PwsC= GtM2*(-wgC(n)*GsC*Pws+2*Re(wg(n)*Gs*PwsC));
id wg(n?)*Gt*GsM2*PwtC= GsM2*(-wgC(n)*GtC*Pwt+2*Re(wg(n)*Gt*PwtC));
id wg(n?)*Gs*GtM2= GtM2*(-GsC*wgC(n)+2*Re(Gs*wg(n)));
id wg(n?)*Gt*GsM2= GsM2*(-GtC*wgC(n)+2*Re(Gt*wg(n)));
*
id GsM2*GtC*PwtC*wgC(n?)= GsM2*(-Gt*Pwt*wg(n)+2*Re(Gt*Pwt*wg(n)));
id GtM2*GsC*PwsC*wgC(n?)= GtM2*(-Gs*Pws*wg(n)+2*Re(Gs*Pws*wg(n)));
*
id GsM2*GtC*PwtC= GsM2*(-Gt*Pwt+2*Re(Gt*Pwt));
id GtM2*GsC*PwsC= GtM2*(-Gs*Pws+2*Re(Gs*Pws));
*
id GsC*PwsC*Gt*Pwt= -Gs*Pws*GtC*PwtC+2*Re(GsC*Gt*PwsC*Pwt);
id GsC*Gt*Pws*PwtC= -Gs*GtC*PwsC*Pwt+2*Re(Gs*GtC*PwsC*Pwt);
*
id Gs*GtC*Pwt= -GsC*Gt*PwtC+2*Re(Gs*GtC*Pwt);
*
id GtM2*PwsC= GtM2*(-Pws+2*Re(Pws));
id GsM2*PwtC= GsM2*(-Pwt+2*Re(Pwt));
*
id Re(wgC(2)*wg(1))= Re(wg(2)*wgC(1));
al Re(wgC(3)*wg(1))= Re(wg(3)*wgC(1));
al Re(wgC(4)*wg(1))= Re(wg(4)*wgC(1));
al Re(wgC(5)*wg(1))= Re(wg(5)*wgC(1));
al Re(wgC(6)*wg(1))= Re(wg(6)*wgC(1));
al Re(wgC(7)*wg(1))= Re(wg(7)*wgC(1));
al Re(wgC(8)*wg(1))= Re(wg(8)*wgC(1));
al Re(wgC(9)*wg(1))= Re(wg(9)*wgC(1));
al Re(wgC(10)*wg(1))= Re(wg(10)*wgC(1));
*
id Re(wgC(3)*wg(2))= Re(wg(3)*wgC(2));
al Re(wgC(4)*wg(2))= Re(wg(4)*wgC(2));
al Re(wgC(5)*wg(2))= Re(wg(5)*wgC(2));
al Re(wgC(6)*wg(2))= Re(wg(6)*wgC(2));
al Re(wgC(7)*wg(2))= Re(wg(7)*wgC(2));
al Re(wgC(8)*wg(2))= Re(wg(8)*wgC(2));
al Re(wgC(9)*wg(2))= Re(wg(9)*wgC(2));
al Re(wgC(10)*wg(2))= Re(wg(10)*wgC(2));
*
id Re(wgC(4)*wg(3))= Re(wg(4)*wgC(3));
al Re(wgC(5)*wg(3))= Re(wg(5)*wgC(3));
al Re(wgC(6)*wg(3))= Re(wg(6)*wgC(3));
al Re(wgC(7)*wg(3))= Re(wg(7)*wgC(3));
al Re(wgC(8)*wg(3))= Re(wg(8)*wgC(3));
al Re(wgC(9)*wg(3))= Re(wg(9)*wgC(3));
al Re(wgC(10)*wg(3))= Re(wg(10)*wgC(3));
*
id Re(wgC(5)*wg(4))= Re(wg(5)*wgC(4));
al Re(wgC(6)*wg(4))= Re(wg(6)*wgC(4));
al Re(wgC(7)*wg(4))= Re(wg(7)*wgC(4));
al Re(wgC(8)*wg(4))= Re(wg(8)*wgC(4));
al Re(wgC(9)*wg(4))= Re(wg(9)*wgC(4));
al Re(wgC(10)*wg(4))= Re(wg(10)*wgC(4));
*
id Re(wgC(6)*wg(5))= Re(wg(6)*wgC(5));
al Re(wgC(7)*wg(5))= Re(wg(7)*wgC(5));
al Re(wgC(8)*wg(5))= Re(wg(8)*wgC(5));
al Re(wgC(9)*wg(5))= Re(wg(9)*wgC(5));
al Re(wgC(10)*wg(5))= Re(wg(10)*wgC(5));
*
id Re(wgC(7)*wg(6))= Re(wg(7)*wgC(6));
al Re(wgC(8)*wg(6))= Re(wg(8)*wgC(6));
al Re(wgC(9)*wg(6))= Re(wg(9)*wgC(6));
al Re(wgC(10)*wg(6))= Re(wg(10)*wgC(6));
*
id Re(wgC(8)*wg(7))= Re(wg(8)*wgC(7));
al Re(wgC(9)*wg(7))= Re(wg(9)*wgC(7));
al Re(wgC(10)*wg(7))= Re(wg(10)*wgC(7));
*
id Re(wgC(9)*wg(8))= Re(wg(9)*wgC(8));
al Re(wgC(10)*wg(8))= Re(wg(10)*wgC(8));
*
id Re(wgC(10)*wg(9))= Re(wg(10)*wgC(9));
*
.store
*
#endprocedure
*
#procedure contractgg(dummy)
*
g Wd_gg= T_gg*d_(mu,nu);
g Wpp_gg= T_gg*pp(mu)*pp(nu);
*
id rs^-2= -1/(X*y*s);
id pp.pp= 0;
al qp.qp= 0;
al Qm.Qm= QmS;
id Qp.Qm= pp.Qm-qp.Qm;
id kb.Qp= kb.pp-kb.qp;
id k.Qp= k.pp-k.qp;
id pp.Qp= -pp.qp;
id Qp.Qp= -2*pp.qp;
*
id pp.Qm= -1/2*(QmS+y*s);
al k.kb= 1/2*(-sp+ums+dms);
al Qm.k= 1/2*(QmS+t-ums);
al pp.kb= 1/2*(tp-dms);
al Qm.kb= 1/2*(QmS+u-dms);
al pp.k= 1/2*(up-ums);
al pp.qp= 1/2*kp;
al Qm.qp= 1/2*(QmS+km);
al kb.qp= 1/2*(dms-zp);
al k.qp= 1/2*(ums-zm);
*
id t= tau*y*s;
al sp= x2*y*s;
al km= z*y*s;
al km^-1= 1/(z*y*s);
al zm= -(x2-x1)*y*s;
id QmS= X*y*s;
*
ratio,tuM,tau,aux;
id aux= ums/(y*s);
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
id kp= y*s * (  - 1 - X + x2 - z );
al u= y*s * (  - 1 - tau - 2*X - z )
       + ums + dms;
al up= y*s * (  - tau - X - x1 )
       + 2*ums + dms;
al tp= y*s * ( tau + X + x1 - x2 + z )
       - ums;
al zp= y*s * ( 1 - x1 )
       + ums + dms;
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
.store
*
g GG10= (-1/2*Wd_gg-2*X/Xp^2*Wpp_gg/(y*s))*sel0*s;
g GG20= (-1/2*Wd_gg-6*X/Xp^2*Wpp_gg/(y*s))*sel0*s;
g GG11= (-1/2*Wd_gg-2*X/Xp^2*Wpp_gg/(y*s))*selu*s;
g GG21= (-1/2*Wd_gg-6*X/Xp^2*Wpp_gg/(y*s))*selu*s;
g GG12= (-1/2*Wd_gg-2*X/Xp^2*Wpp_gg/(y*s))*seluu*s;
g GG22= (-1/2*Wd_gg-6*X/Xp^2*Wpp_gg/(y*s))*seluu*s;
*
id tuM^-2*seluu= 1;
id seluu= 0;
*
id tuM^-2*selu= 0;
id tuM^-1*selu= 1;
id selu= 0;
*
id tuM^-2*sel0= 0;
id tuM^-1*sel0= 0;
*
id sel0= 1;
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
*--
id Re(wg(n?)*wgC(2))= 1/2*(Re(wg(n)*wgC([2+3]))+Re(wg(n)*wgC([2-3])));
id Re(wg(n?)*wgC(3))= 1/2*(Re(wg(n)*wgC([2+3]))-Re(wg(n)*wgC([2-3])));
id Re(wg(2)*wgC(2))= wgM2(2);
id Re(wg(3)*wgC(3))= wgM2(3);
*
id Re(wg(2)*wgC(n?))= 1/2*(Re(wg([2+3])*wgC(n))+Re(wg([2-3])*wgC(n)));
id Re(wg(3)*wgC(n?))= 1/2*(Re(wg([2+3])*wgC(n))-Re(wg([2-3])*wgC(n)));
id Re(wg(2)*wgC(2))= wgM2(2);
id Re(wg(3)*wgC(3))= wgM2(3);
*
id Re(wg(n?)*wgC(5))= 1/2*(Re(wg(n)*wgC([5+6]))+Re(wg(n)*wgC([5-6])));
id Re(wg(n?)*wgC(6))= 1/2*(Re(wg(n)*wgC([5+6]))-Re(wg(n)*wgC([5-6])));
id Re(wg(5)*wgC(5))= wgM2(5);
id Re(wg(6)*wgC(6))= wgM2(6);
*
id Re(wg(5)*wgC(n?))= 1/2*(Re(wg([5+6])*wgC(n))+Re(wg([5-6])*wgC(n)));
id Re(wg(6)*wgC(n?))= 1/2*(Re(wg([5+6])*wgC(n))-Re(wg([5-6])*wgC(n)));
id Re(wg(5)*wgC(5))= wgM2(5);
id Re(wg(6)*wgC(6))= wgM2(6);
*
id Re(wg(2)*X?*Y?)= 1/2*(Re(wg([2+3])*X*Y)+Re(wg([2-3])*X*Y));
id Re(wg(3)*X?*Y?)= 1/2*(Re(wg([2+3])*X*Y)-Re(wg([2-3])*X*Y));
id Re(wg(5)*X?*Y?)= 1/2*(Re(wg([5+6])*X*Y)+Re(wg([5-6])*X*Y));
id Re(wg(6)*X?*Y?)= 1/2*(Re(wg([5+6])*X*Y)-Re(wg([5-6])*X*Y));
*
id Re(wg(2)*X?)= 1/2*(Re(wg([2+3])*X)+Re(wg([2-3])*X));
id Re(wg(3)*X?)= 1/2*(Re(wg([2+3])*X)-Re(wg([2-3])*X));
id Re(wg(5)*X?)= 1/2*(Re(wg([5+6])*X)+Re(wg([5-6])*X));
id Re(wg(6)*X?)= 1/2*(Re(wg([5+6])*X)-Re(wg([5-6])*X));
*
id Re(wg([2+3])*wgC([2-3]))= Re(wg([2-3])*wgC([2+3]));
id Re(wg([2+3])*wgC([5-6]))= Re(wg([5-6])*wgC([2+3]));
id Re(wg([2-3])*wgC([5+6]))= Re(wg([5+6])*wgC([2-3]));
id Re(wg([5+6])*wgC([5-6]))= Re(wg([5-6])*wgC([5+6]));
*
id Re(wg([2+3])*wgC([2+3]))= wgM2(2)+wgM2(3)+2*Re(wg(2)*wgC(3));
id Re(wg([2-3])*wgC([2-3]))= wgM2(2)+wgM2(3)-2*Re(wg(2)*wgC(3));
id Re(wg([5+6])*wgC([5+6]))= wgM2(5)+wgM2(6)+2*Re(wg(5)*wgC(6));
id Re(wg([5-6])*wgC([5-6]))= wgM2(5)+wgM2(6)-2*Re(wg(5)*wgC(6));
*--
id X= [1+X]-1;
id x1= 1-[1-x1];
id x2= 1-[1-x1]-[x1-x2];
id z= 1-[1-x1]-[x1-x2]-[x2-z];
id tau= [tau+X+x1]-[1+X]+[1-x1];
id Xp^-1= 1/[1+X];
ratio,X,[1+X],one;
id one= 1;
al one^-1= 1;
*
id Re(wg(X?)*wgC(Y?))= Rep(X,Y);
id GsM2*GtM2= GstM2;
id [1+X] = Xp;
al [1+X]^-1= Xpi;
id [1-x1]= omx1;
al [x1-x2]= x12;
al [x2-z]= x2z;
al [tau+X+x1]= tx;
al zX^-1= zXi;
id Xpi^2= Xpi-aux;
id aux= XX;
*
id Xp^2= Xps;
id Xpi^2= Xpis;
id y^3= yc;
id y^2= ys;
*
id omx1= O;
id O^4= Oq;
id O^3= Oc;
id O^2= Os;
id x12^3= x12c;
id x12^2= x12s;
id x2z^4= x2zq;
id x2z^3= x2zc;
id x2z^2= x2zs;
id tx^2= txs;
id vv^2= vvs;
id vv^3= vvc;
id vvs*vv= vvc;
id vvs^2= vvq;
id vvc*vv= vvq;
id z^-1= zi;
*
id Re(Gs*GtC*Pwt*PwsC)= rl1;
id Re(Gs*GtC*Pwt)= rl2;
id Re(Pws)= rl3;
id Re(wg(n?)*Gs*PwsC)= rs1wg(n);
id Re(wg(n?)*Gs)= rs2wg(n);
id Re(wg(n?)*Gt*PwtC)= rt1wg(n);
id Re(wg(n?)*Gt)= rt2wg(n);
id Rep(n1?,n2?)= rwwg(n1,n2);
id Re(wg(n1?)*wgC(n2?))= rwwg(n1,n2);
id wgM2(n?)= rwwg(n,n);
*
id fun?([2+3],n?)= fun(11,n);
al fun?([2-3],n?)= fun(12,n);
al fun?([5+6],n?)= fun(13,n);
al fun?([5-6],n?)= fun(14,n);
*
id fun?(n?,[2+3])= fun(n,11);
al fun?(n?,[2-3])= fun(n,12);
al fun?(n?,[5+6])= fun(n,13);
al fun?(n?,[5-6])= fun(n,14);
*
id fun?([2+3])= fun(11);
al fun?([2-3])= fun(12);
al fun?([5+6])= fun(13);
al fun?([5-6])= fun(14);
*
.store
*
*----g^2--------------------
*
g GG10s= GG10;
g GG20s= GG20;
g GG11s= GG11;
g GG21s= GG21;
g GG12s= GG12;
g GG22s= GG22;
*
id gexp^4= 0;
id gexp^3= 0;
id gexp^2= 1;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wg,wgC,wgM2,rs1wg,rs2wg,
  rt1wg,rt2wg,rwwg,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
*  tx,txs,x2z,x2zs,x2zc,x2zq,x12,x12s,x12c,zXi,XX;
  zXi,XX,vv,vvs,vvc,vvq;
format fortran;
print;
*
.store
*
*----g^3--------------------
*
g GG10c= GG10;
g GG20c= GG20;
g GG11c= GG11;
g GG21c= GG21;
g GG12c= GG12;
g GG22c= GG22;
*
id gexp^4= 0;
id gexp^3= 1;
id gexp^2= 0;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wg,wgC,wgM2,rs1wg,rs2wg,
  rt1wg,rt2wg,rwwg,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
*  tx,txs,x2z,x2zs,x2zc,x2zq,x12,x12s,x12c,zXi,XX;
  zXi,XX,vv,vvs,vvc,vvq;
format fortran;
print;
*
*
.store
*
*----g^4--------------------
*
g GG10q= GG10/GstM2;
g GG20q= GG20/GstM2;
g GG11q= GG11/GstM2;
g GG21q= GG21/GstM2;
g GG12q= GG12/GstM2;
g GG22q= GG22/GstM2;
*
id gexp^4= 1;
id gexp^3= 0;
id gexp^2= 0;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wg,wgC,wgM2,rs1wg,rs2wg,
  rt1wg,rt2wg,rwwg,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
*  tx,txs,x2z,x2zs,x2zc,x2zq,x12,x12s,x12c,zXi,XX;
  zXi,XX,vv,vvs,vvc,vvq;
format fortran;
print;
*
#endprocedure
*
#procedure totalgz(dummy)
*
* Now the full stuff
*
load wrl$grp04:[giampiero]outputampgz_A.sav;
*
g T_gz= PwsM2*PwtM2*W_gz_int6/(egt*I0);
*nwrite statistics;
*
id es(T)*GT= egt;
*
id wg(n?)= wg(n)*rs;
al wgC(n?)= wgC(n)*rs;
id wz(n?)= wz(n)*rs;
al wzC(n?)= wzC(n)*rs;
*
id wg(4)= X*wg(4);
al wg(7)= X*wg(7);
al wg(8)= X*wg(8);
al wg(10)= X*wg(10);
id wgC(4)= X*wgC(4);
al wgC(7)= X*wgC(7);
al wgC(8)= X*wgC(8);
al wgC(10)= X*wgC(10);
*
id wz(4)= X*wz(4);
al wz(7)= X*wz(7);
al wz(8)= X*wz(8);
al wz(10)= X*wz(10);
id wzC(4)= X*wzC(4);
al wzC(7)= X*wzC(7);
al wzC(8)= X*wzC(8);
al wzC(10)= X*wzC(10);
*
id wg(n?)*wgC(n?)= wgM2(n);
id wz(n?)*wzC(n?)= wzM2(n);
*
id GtM2*GsM2*wg(2)*wzC(1)= GtM2*GsM2*(-wgC(2)*wz(1)+
                                 2*Re(wg(2)*wzC(1)));
al GtM2*GsM2*wg(3)*wzC(1)= GtM2*GsM2*(-wgC(3)*wz(1)+
                                 2*Re(wg(3)*wzC(1)));
al GtM2*GsM2*wg(4)*wzC(1)= GtM2*GsM2*(-wgC(4)*wz(1)+
                                 2*Re(wg(4)*wzC(1)));
al GtM2*GsM2*wg(5)*wzC(1)= GtM2*GsM2*(-wgC(5)*wz(1)+
                                 2*Re(wg(5)*wzC(1)));
al GtM2*GsM2*wg(6)*wzC(1)= GtM2*GsM2*(-wgC(6)*wz(1)+
                                 2*Re(wg(6)*wzC(1)));
al GtM2*GsM2*wg(7)*wzC(1)= GtM2*GsM2*(-wgC(7)*wz(1)+
                                 2*Re(wg(7)*wzC(1)));
al GtM2*GsM2*wg(8)*wzC(1)= GtM2*GsM2*(-wgC(8)*wz(1)+
                                 2*Re(wg(8)*wzC(1)));
al GtM2*GsM2*wg(9)*wzC(1)= GtM2*GsM2*(-wgC(9)*wz(1)+
                                 2*Re(wg(9)*wzC(1)));
al GtM2*GsM2*wg(10)*wzC(1)= GtM2*GsM2*(-wgC(10)*wz(1)+
                                 2*Re(wg(10)*wzC(1)));
*
al GtM2*GsM2*wg(3)*wzC(2)= GtM2*GsM2*(-wgC(3)*wz(2)+
                                 2*Re(wg(3)*wzC(2)));
al GtM2*GsM2*wg(4)*wzC(2)= GtM2*GsM2*(-wgC(4)*wz(2)+
                                 2*Re(wg(4)*wzC(2)));
al GtM2*GsM2*wg(5)*wzC(2)= GtM2*GsM2*(-wgC(5)*wz(2)+
                                 2*Re(wg(5)*wzC(2)));
al GtM2*GsM2*wg(6)*wzC(2)= GtM2*GsM2*(-wgC(6)*wz(2)+
                                 2*Re(wg(6)*wzC(2)));
al GtM2*GsM2*wg(7)*wzC(2)= GtM2*GsM2*(-wgC(7)*wz(2)+
                                 2*Re(wg(7)*wzC(2)));
al GtM2*GsM2*wg(8)*wzC(2)= GtM2*GsM2*(-wgC(8)*wz(2)+
                                 2*Re(wg(8)*wzC(2)));
al GtM2*GsM2*wg(9)*wzC(2)= GtM2*GsM2*(-wgC(9)*wz(2)+
                                 2*Re(wg(9)*wzC(2)));
al GtM2*GsM2*wg(10)*wzC(2)= GtM2*GsM2*(-wgC(10)*wz(2)+
                                 2*Re(wg(10)*wzC(2)));
*
al GtM2*GsM2*wg(4)*wzC(3)= GtM2*GsM2*(-wgC(4)*wz(3)+
                                 2*Re(wg(4)*wzC(3)));
al GtM2*GsM2*wg(5)*wzC(3)= GtM2*GsM2*(-wgC(5)*wz(3)+
                                 2*Re(wg(5)*wzC(3)));
al GtM2*GsM2*wg(6)*wzC(3)= GtM2*GsM2*(-wgC(6)*wz(3)+
                                 2*Re(wg(6)*wzC(3)));
al GtM2*GsM2*wg(7)*wzC(3)= GtM2*GsM2*(-wgC(7)*wz(3)+
                                 2*Re(wg(7)*wzC(3)));
al GtM2*GsM2*wg(8)*wzC(3)= GtM2*GsM2*(-wgC(8)*wz(3)+
                                 2*Re(wg(8)*wzC(3)));
al GtM2*GsM2*wg(9)*wzC(3)= GtM2*GsM2*(-wgC(9)*wz(3)+
                                 2*Re(wg(9)*wzC(3)));
al GtM2*GsM2*wg(10)*wzC(3)= GtM2*GsM2*(-wgC(10)*wz(3)+
                                 2*Re(wg(10)*wzC(3)));
*
al GtM2*GsM2*wg(5)*wzC(4)= GtM2*GsM2*(-wgC(5)*wz(4)+
                                 2*Re(wg(5)*wzC(4)));
al GtM2*GsM2*wg(6)*wzC(4)= GtM2*GsM2*(-wgC(6)*wz(4)+
                                 2*Re(wg(6)*wzC(4)));
al GtM2*GsM2*wg(7)*wzC(4)= GtM2*GsM2*(-wgC(7)*wz(4)+
                                 2*Re(wg(7)*wzC(4)));
al GtM2*GsM2*wg(8)*wzC(4)= GtM2*GsM2*(-wgC(8)*wz(4)+
                                 2*Re(wg(8)*wzC(4)));
al GtM2*GsM2*wg(9)*wzC(4)= GtM2*GsM2*(-wgC(9)*wz(4)+
                                 2*Re(wg(9)*wzC(4)));
al GtM2*GsM2*wg(10)*wzC(4)= GtM2*GsM2*(-wgC(10)*wz(4)+
                                 2*Re(wg(10)*wzC(4)));
*
al GtM2*GsM2*wg(6)*wzC(5)= GtM2*GsM2*(-wgC(6)*wz(5)+
                                 2*Re(wg(6)*wzC(5)));
al GtM2*GsM2*wg(7)*wzC(5)= GtM2*GsM2*(-wgC(7)*wz(5)+
                                 2*Re(wg(7)*wzC(5)));
al GtM2*GsM2*wg(8)*wzC(5)= GtM2*GsM2*(-wgC(8)*wz(5)+
                                 2*Re(wg(8)*wzC(5)));
al GtM2*GsM2*wg(9)*wzC(5)= GtM2*GsM2*(-wgC(9)*wz(5)+
                                 2*Re(wg(9)*wzC(5)));
al GtM2*GsM2*wg(10)*wzC(5)= GtM2*GsM2*(-wgC(10)*wz(5)+
                                 2*Re(wg(10)*wzC(5)));
*
al GtM2*GsM2*wg(7)*wzC(6)= GtM2*GsM2*(-wgC(7)*wz(6)+
                                 2*Re(wg(7)*wzC(6)));
al GtM2*GsM2*wg(8)*wzC(6)= GtM2*GsM2*(-wgC(8)*wz(6)+
                                 2*Re(wg(8)*wzC(6)));
al GtM2*GsM2*wg(9)*wzC(6)= GtM2*GsM2*(-wgC(9)*wz(6)+
                                 2*Re(wg(9)*wzC(6)));
al GtM2*GsM2*wg(10)*wzC(6)= GtM2*GsM2*(-wgC(10)*wz(6)+
                                 2*Re(wg(10)*wzC(6)));
*
al GtM2*GsM2*wg(8)*wzC(7)= GtM2*GsM2*(-wgC(8)*wz(7)+
                                 2*Re(wg(8)*wzC(7)));
al GtM2*GsM2*wg(9)*wzC(7)= GtM2*GsM2*(-wgC(9)*wz(7)+
                                 2*Re(wg(9)*wzC(7)));
al GtM2*GsM2*wg(10)*wzC(7)= GtM2*GsM2*(-wgC(10)*wz(7)+
                                 2*Re(wg(10)*wzC(7)));
*
al GtM2*GsM2*wg(9)*wzC(8)= GtM2*GsM2*(-wgC(9)*wz(8)+
                                 2*Re(wg(9)*wzC(8)));
al GtM2*GsM2*wg(10)*wzC(8)= GtM2*GsM2*(-wgC(10)*wz(8)+
                                 2*Re(wg(10)*wzC(8)));
*
al GtM2*GsM2*wg(10)*wzC(9)= GtM2*GsM2*(-wgC(10)*wz(9)+
                                 2*Re(wg(10)*wzC(9)));
*
id GtM2*GsM2*wz(2)*wgC(1)= GtM2*GsM2*(-wzC(2)*wg(1)+
                                 2*Re(wz(2)*wgC(1)));
al GtM2*GsM2*wz(3)*wgC(1)= GtM2*GsM2*(-wzC(3)*wg(1)+
                                 2*Re(wz(3)*wgC(1)));
al GtM2*GsM2*wz(4)*wgC(1)= GtM2*GsM2*(-wzC(4)*wg(1)+
                                 2*Re(wz(4)*wgC(1)));
al GtM2*GsM2*wz(5)*wgC(1)= GtM2*GsM2*(-wzC(5)*wg(1)+
                                 2*Re(wz(5)*wgC(1)));
al GtM2*GsM2*wz(6)*wgC(1)= GtM2*GsM2*(-wzC(6)*wg(1)+
                                 2*Re(wz(6)*wgC(1)));
al GtM2*GsM2*wz(7)*wgC(1)= GtM2*GsM2*(-wzC(7)*wg(1)+
                                 2*Re(wz(7)*wgC(1)));
al GtM2*GsM2*wz(8)*wgC(1)= GtM2*GsM2*(-wzC(8)*wg(1)+
                                 2*Re(wz(8)*wgC(1)));
al GtM2*GsM2*wz(9)*wgC(1)= GtM2*GsM2*(-wzC(9)*wg(1)+
                                 2*Re(wz(9)*wgC(1)));
al GtM2*GsM2*wz(10)*wgC(1)= GtM2*GsM2*(-wzC(10)*wg(1)+
                                 2*Re(wz(10)*wgC(1)));
*
al GtM2*GsM2*wz(3)*wgC(2)= GtM2*GsM2*(-wzC(3)*wg(2)+
                                 2*Re(wz(3)*wgC(2)));
al GtM2*GsM2*wz(4)*wgC(2)= GtM2*GsM2*(-wzC(4)*wg(2)+
                                 2*Re(wz(4)*wgC(2)));
al GtM2*GsM2*wz(5)*wgC(2)= GtM2*GsM2*(-wzC(5)*wg(2)+
                                 2*Re(wz(5)*wgC(2)));
al GtM2*GsM2*wz(6)*wgC(2)= GtM2*GsM2*(-wzC(6)*wg(2)+
                                 2*Re(wz(6)*wgC(2)));
al GtM2*GsM2*wz(7)*wgC(2)= GtM2*GsM2*(-wzC(7)*wg(2)+
                                 2*Re(wz(7)*wgC(2)));
al GtM2*GsM2*wz(8)*wgC(2)= GtM2*GsM2*(-wzC(8)*wg(2)+
                                 2*Re(wz(8)*wgC(2)));
al GtM2*GsM2*wz(9)*wgC(2)= GtM2*GsM2*(-wzC(9)*wg(2)+
                                 2*Re(wz(9)*wgC(2)));
al GtM2*GsM2*wz(10)*wgC(2)= GtM2*GsM2*(-wzC(10)*wg(2)+
                                 2*Re(wz(10)*wgC(2)));
*
al GtM2*GsM2*wz(4)*wgC(3)= GtM2*GsM2*(-wzC(4)*wg(3)+
                                 2*Re(wz(4)*wgC(3)));
al GtM2*GsM2*wz(5)*wgC(3)= GtM2*GsM2*(-wzC(5)*wg(3)+
                                 2*Re(wz(5)*wgC(3)));
al GtM2*GsM2*wz(6)*wgC(3)= GtM2*GsM2*(-wzC(6)*wg(3)+
                                 2*Re(wz(6)*wgC(3)));
al GtM2*GsM2*wz(7)*wgC(3)= GtM2*GsM2*(-wzC(7)*wg(3)+
                                 2*Re(wz(7)*wgC(3)));
al GtM2*GsM2*wz(8)*wgC(3)= GtM2*GsM2*(-wzC(8)*wg(3)+
                                 2*Re(wz(8)*wgC(3)));
al GtM2*GsM2*wz(9)*wgC(3)= GtM2*GsM2*(-wzC(9)*wg(3)+
                                 2*Re(wz(9)*wgC(3)));
al GtM2*GsM2*wz(10)*wgC(3)= GtM2*GsM2*(-wzC(10)*wg(3)+
                                 2*Re(wz(10)*wgC(3)));
*
al GtM2*GsM2*wz(5)*wgC(4)= GtM2*GsM2*(-wzC(5)*wg(4)+
                                 2*Re(wz(5)*wgC(4)));
al GtM2*GsM2*wz(6)*wgC(4)= GtM2*GsM2*(-wzC(6)*wg(4)+
                                 2*Re(wz(6)*wgC(4)));
al GtM2*GsM2*wz(7)*wgC(4)= GtM2*GsM2*(-wzC(7)*wg(4)+
                                 2*Re(wz(7)*wgC(4)));
al GtM2*GsM2*wz(8)*wgC(4)= GtM2*GsM2*(-wzC(8)*wg(4)+
                                 2*Re(wz(8)*wgC(4)));
al GtM2*GsM2*wz(9)*wgC(4)= GtM2*GsM2*(-wzC(9)*wg(4)+
                                 2*Re(wz(9)*wgC(4)));
al GtM2*GsM2*wz(10)*wgC(4)= GtM2*GsM2*(-wzC(10)*wg(4)+
                                 2*Re(wz(10)*wgC(4)));
*
al GtM2*GsM2*wz(6)*wgC(5)= GtM2*GsM2*(-wzC(6)*wg(5)+
                                 2*Re(wz(6)*wgC(5)));
al GtM2*GsM2*wz(7)*wgC(5)= GtM2*GsM2*(-wzC(7)*wg(5)+
                                 2*Re(wz(7)*wgC(5)));
al GtM2*GsM2*wz(8)*wgC(5)= GtM2*GsM2*(-wzC(8)*wg(5)+
                                 2*Re(wz(8)*wgC(5)));
al GtM2*GsM2*wz(9)*wgC(5)= GtM2*GsM2*(-wzC(9)*wg(5)+
                                 2*Re(wz(9)*wgC(5)));
al GtM2*GsM2*wz(10)*wgC(5)= GtM2*GsM2*(-wzC(10)*wg(5)+
                                 2*Re(wz(10)*wgC(5)));
*
al GtM2*GsM2*wz(7)*wgC(6)= GtM2*GsM2*(-wzC(7)*wg(6)+
                                 2*Re(wz(7)*wgC(6)));
al GtM2*GsM2*wz(8)*wgC(6)= GtM2*GsM2*(-wzC(8)*wg(6)+
                                 2*Re(wz(8)*wgC(6)));
al GtM2*GsM2*wz(9)*wgC(6)= GtM2*GsM2*(-wzC(9)*wg(6)+
                                 2*Re(wz(9)*wgC(6)));
al GtM2*GsM2*wz(10)*wgC(6)= GtM2*GsM2*(-wzC(10)*wg(6)+
                                 2*Re(wz(10)*wgC(6)));
*
al GtM2*GsM2*wz(8)*wgC(7)= GtM2*GsM2*(-wzC(8)*wg(7)+
                                 2*Re(wz(8)*wgC(7)));
al GtM2*GsM2*wz(9)*wgC(7)= GtM2*GsM2*(-wzC(9)*wg(7)+
                                 2*Re(wz(9)*wgC(7)));
al GtM2*GsM2*wz(10)*wgC(7)= GtM2*GsM2*(-wzC(10)*wg(7)+
                                 2*Re(wz(10)*wgC(7)));
*
al GtM2*GsM2*wz(9)*wgC(8)= GtM2*GsM2*(-wzC(9)*wg(8)+
                                 2*Re(wz(9)*wgC(8)));
al GtM2*GsM2*wz(10)*wgC(8)= GtM2*GsM2*(-wzC(10)*wg(8)+
                                 2*Re(wz(10)*wgC(8)));
*
al GtM2*GsM2*wz(10)*wgC(9)= GtM2*GsM2*(-wzC(10)*wg(9)+
                                 2*Re(wz(10)*wgC(9)));
*
id wg(n?)*Gs*GtM2*PwsC= GtM2*(-wgC(n)*GsC*Pws+2*Re(wg(n)*Gs*PwsC));
id wg(n?)*Gt*GsM2*PwtC= GsM2*(-wgC(n)*GtC*Pwt+2*Re(wg(n)*Gt*PwtC));
id wz(n?)*Gs*GtM2*PwsC= GtM2*(-wzC(n)*GsC*Pws+2*Re(wz(n)*Gs*PwsC));
id wz(n?)*Gt*GsM2*PwtC= GsM2*(-wzC(n)*GtC*Pwt+2*Re(wz(n)*Gt*PwtC));
*
id GsM2*GtC*PwtC*wgC(n?)= GsM2*(-Gt*Pwt*wg(n)+2*Re(Gt*Pwt*wg(n)));
id GtM2*GsC*PwsC*wgC(n?)= GtM2*(-Gs*Pws*wg(n)+2*Re(Gs*Pws*wg(n)));
id GsM2*GtC*PwtC*wzC(n?)= GsM2*(-Gt*Pwt*wz(n)+2*Re(Gt*Pwt*wz(n)));
id GtM2*GsC*PwsC*wzC(n?)= GtM2*(-Gs*Pws*wz(n)+2*Re(Gs*Pws*wz(n)));
*
id wg(n?)*Gs*GtM2= GtM2*(-GsC*wgC(n)+2*Re(Gs*wg(n)));
id wg(n?)*Gt*GsM2= GsM2*(-GtC*wgC(n)+2*Re(Gt*wg(n)));
id wz(n?)*Gs*GtM2= GtM2*(-GsC*wzC(n)+2*Re(Gs*wz(n)));
id wz(n?)*Gt*GsM2= GsM2*(-GtC*wzC(n)+2*Re(Gt*wz(n)));
*
id GsM2*GtC*PwtC= GsM2*(-Gt*Pwt+2*Re(Gt*Pwt));
id GtM2*GsC*PwsC= GtM2*(-Gs*Pws+2*Re(Gs*Pws));
*
id GsC*PwsC*Gt*Pwt= -Gs*Pws*GtC*PwtC+2*Re(GsC*Gt*PwsC*Pwt);
id GsC*Gt*Pws*PwtC= -Gs*GtC*PwsC*Pwt+2*Re(Gs*GtC*PwsC*Pwt);
id GsC*Gt*PwsC*PwtC= -Gs*GtC*Pws*Pwt+2*Re(Gs*GtC*Pws*Pwt);
*
id Gs*GtC*Pwt= -GsC*Gt*PwtC+2*Re(Gs*GtC*Pwt);
*
id GtM2*PwsC= GtM2*(-Pws+2*Re(Pws));
id GsM2*PwtC= GsM2*(-Pwt+2*Re(Pwt));
*
id Re(wzC(2)*wg(1))= Re(wz(2)*wgC(1));
al Re(wzC(3)*wg(1))= Re(wz(3)*wgC(1));
al Re(wzC(4)*wg(1))= Re(wz(4)*wgC(1));
al Re(wzC(5)*wg(1))= Re(wz(5)*wgC(1));
al Re(wzC(6)*wg(1))= Re(wz(6)*wgC(1));
al Re(wzC(7)*wg(1))= Re(wz(7)*wgC(1));
al Re(wzC(8)*wg(1))= Re(wz(8)*wgC(1));
al Re(wzC(9)*wg(1))= Re(wz(9)*wgC(1));
al Re(wzC(10)*wg(1))= Re(wz(10)*wgC(1));
*
id Re(wgC(2)*wz(1))= Re(wg(2)*wzC(1));
al Re(wgC(3)*wz(1))= Re(wg(3)*wzC(1));
al Re(wgC(4)*wz(1))= Re(wg(4)*wzC(1));
al Re(wgC(5)*wz(1))= Re(wg(5)*wzC(1));
al Re(wgC(6)*wz(1))= Re(wg(6)*wzC(1));
al Re(wgC(7)*wz(1))= Re(wg(7)*wzC(1));
al Re(wgC(8)*wz(1))= Re(wg(8)*wzC(1));
al Re(wgC(9)*wz(1))= Re(wg(9)*wzC(1));
al Re(wgC(10)*wz(1))= Re(wg(10)*wzC(1));
*
id Re(wgC(3)*wz(2))= Re(wg(3)*wzC(2));
al Re(wgC(4)*wz(2))= Re(wg(4)*wzC(2));
al Re(wgC(5)*wz(2))= Re(wg(5)*wzC(2));
al Re(wgC(6)*wz(2))= Re(wg(6)*wzC(2));
al Re(wgC(7)*wz(2))= Re(wg(7)*wzC(2));
al Re(wgC(8)*wz(2))= Re(wg(8)*wzC(2));
al Re(wgC(9)*wz(2))= Re(wg(9)*wzC(2));
al Re(wgC(10)*wz(2))= Re(wg(10)*wzC(2));
*
id Re(wzC(3)*wg(2))= Re(wz(3)*wgC(2));
al Re(wzC(4)*wg(2))= Re(wz(4)*wgC(2));
al Re(wzC(5)*wg(2))= Re(wz(5)*wgC(2));
al Re(wzC(6)*wg(2))= Re(wz(6)*wgC(2));
al Re(wzC(7)*wg(2))= Re(wz(7)*wgC(2));
al Re(wzC(8)*wg(2))= Re(wz(8)*wgC(2));
al Re(wzC(9)*wg(2))= Re(wz(9)*wgC(2));
al Re(wzC(10)*wg(2))= Re(wz(10)*wgC(2));
*
id Re(wgC(4)*wz(3))= Re(wg(4)*wzC(3));
al Re(wgC(5)*wz(3))= Re(wg(5)*wzC(3));
al Re(wgC(6)*wz(3))= Re(wg(6)*wzC(3));
al Re(wgC(7)*wz(3))= Re(wg(7)*wzC(3));
al Re(wgC(8)*wz(3))= Re(wg(8)*wzC(3));
al Re(wgC(9)*wz(3))= Re(wg(9)*wzC(3));
al Re(wgC(10)*wz(3))= Re(wg(10)*wzC(3));
*
id Re(wzC(4)*wg(3))= Re(wz(4)*wgC(3));
al Re(wzC(5)*wg(3))= Re(wz(5)*wgC(3));
al Re(wzC(6)*wg(3))= Re(wz(6)*wgC(3));
al Re(wzC(7)*wg(3))= Re(wz(7)*wgC(3));
al Re(wzC(8)*wg(3))= Re(wz(8)*wgC(3));
al Re(wzC(9)*wg(3))= Re(wz(9)*wgC(3));
al Re(wzC(10)*wg(3))= Re(wz(10)*wgC(3));
*
id Re(wgC(5)*wz(4))= Re(wg(5)*wzC(4));
al Re(wgC(6)*wz(4))= Re(wg(6)*wzC(4));
al Re(wgC(7)*wz(4))= Re(wg(7)*wzC(4));
al Re(wgC(8)*wz(4))= Re(wg(8)*wzC(4));
al Re(wgC(9)*wz(4))= Re(wg(9)*wzC(4));
al Re(wgC(10)*wz(4))= Re(wg(10)*wzC(4));
*
id Re(wzC(5)*wg(4))= Re(wg(5)*wzC(4));
al Re(wzC(6)*wg(4))= Re(wg(6)*wzC(4));
al Re(wzC(7)*wg(4))= Re(wg(7)*wzC(4));
al Re(wzC(8)*wg(4))= Re(wg(8)*wzC(4));
al Re(wzC(9)*wg(4))= Re(wg(9)*wzC(4));
al Re(wzC(10)*wg(4))= Re(wg(10)*wzC(4));
*
id Re(wgC(6)*wz(5))= Re(wg(6)*wzC(5));
al Re(wgC(7)*wz(5))= Re(wg(7)*wzC(5));
al Re(wgC(8)*wz(5))= Re(wg(8)*wzC(5));
al Re(wgC(9)*wz(5))= Re(wg(9)*wzC(5));
al Re(wgC(10)*wz(5))= Re(wg(10)*wzC(5));
*
id Re(wzC(6)*wg(5))= Re(wz(6)*wgC(5));
al Re(wzC(7)*wg(5))= Re(wz(7)*wgC(5));
al Re(wzC(8)*wg(5))= Re(wz(8)*wgC(5));
al Re(wzC(9)*wg(5))= Re(wz(9)*wgC(5));
al Re(wzC(10)*wg(5))= Re(wz(10)*wgC(5));
*
id Re(wgC(7)*wz(6))= Re(wg(7)*wzC(6));
al Re(wgC(8)*wz(6))= Re(wg(8)*wzC(6));
al Re(wgC(9)*wz(6))= Re(wg(9)*wzC(6));
al Re(wgC(10)*wz(6))= Re(wg(10)*wzC(6));
*
id Re(wzC(7)*wg(6))= Re(wz(7)*wgC(6));
al Re(wzC(8)*wg(6))= Re(wz(8)*wgC(6));
al Re(wzC(9)*wg(6))= Re(wz(9)*wgC(6));
al Re(wzC(10)*wg(6))= Re(wz(10)*wgC(6));
*
id Re(wgC(8)*wz(7))= Re(wg(8)*wzC(7));
al Re(wgC(9)*wz(7))= Re(wg(9)*wzC(7));
al Re(wgC(10)*wz(7))= Re(wg(10)*wzC(7));
*
id Re(wzC(8)*wg(7))= Re(wz(8)*wgC(7));
al Re(wzC(9)*wg(7))= Re(wz(9)*wgC(7));
al Re(wzC(10)*wg(7))= Re(wz(10)*wgC(7));
*
id Re(wgC(9)*wz(8))= Re(wg(9)*wzC(8));
al Re(wgC(10)*wz(8))= Re(wg(10)*wzC(8));
*
id Re(wzC(9)*wg(8))= Re(wg(9)*wgC(8));
al Re(wzC(10)*wg(8))= Re(wg(10)*wgC(8));
*
id Re(wgC(10)*wz(9))= Re(wg(10)*wzC(9));
*
id Re(wzC(10)*wg(9))= Re(wz(10)*wgC(9));
*
id wg(n?)*wzC(n?)= -wgC(n)*wz(n)+2*Re(wg(n)*wzC(n));
*
.store
*
#endprocedure
*
#procedure contractgz(dummy)
*
g Wd_gz= T_gz*d_(mu,nu);
g Wpp_gz= T_gz*pp(mu)*pp(nu);
g We_gz= T_gz*e_(mu,nu,Qm,pp);
*
contract;
*
id rs^-1= -1/(X*y*s);
id pp.pp= 0;
al qp.qp= 0;
al Qm.Qm= QmS;
id Qp.Qm= pp.Qm-qp.Qm;
id kb.Qp= kb.pp-kb.qp;
id k.Qp= k.pp-k.qp;
id pp.Qp= -pp.qp;
id Qp.Qp= -2*pp.qp;
id qp.Qp= qp.pp;
*
id pp.Qm= -1/2*(QmS+y*s);
al k.kb= 1/2*(-sp+ums+dms);
al Qm.k= 1/2*(QmS+t-ums);
al pp.kb= 1/2*(tp-dms);
al Qm.kb= 1/2*(QmS+u-dms);
al pp.k= 1/2*(up-ums);
al pp.qp= 1/2*kp;
al Qm.qp= 1/2*(QmS+km);
al kb.qp= 1/2*(dms-zp);
al k.qp= 1/2*(ums-zm);
*
id t= tau*y*s;
al sp= x2*y*s;
al km= z*y*s;
al km^-1= 1/(z*y*s);
al zm= -(x2-x1)*y*s;
id QmS= X*y*s;
*
ratio,tuM,tau,aux;
id aux= ums/(y*s);
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
id kp= y*s * (  - 1 - X + x2 - z );
al u= y*s * (  - 1 - tau - 2*X - z )
       + ums + dms;
al up= y*s * (  - tau - X - x1 )
       + 2*ums + dms;
al tp= y*s * ( tau + X + x1 - x2 + z )
       - ums;
al zp= y*s * ( 1 - x1 )
       + ums + dms;
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
.store
*
g GZ10= (-1/2*Wd_gz-2*X/Xp^2*Wpp_gz/(y*s))*sel0*s^2;
g GZ20= (-1/2*Wd_gz-6*X/Xp^2*Wpp_gz/(y*s))*sel0*s^2;
g GZ30= (We_gz/(y*s)/Xp)*sel0*s^2;
g GZ11= (-1/2*Wd_gz-2*X/Xp^2*Wpp_gz/(y*s))*selu*s^2;
g GZ21= (-1/2*Wd_gz-6*X/Xp^2*Wpp_gz/(y*s))*selu*s^2;
g GZ31= (We_gz/(y*s)/Xp)*selu*s^2;
g GZ12= (-1/2*Wd_gz-2*X/Xp^2*Wpp_gz/(y*s))*seluu*s^2;
g GZ22= (-1/2*Wd_gz-6*X/Xp^2*Wpp_gz/(y*s))*seluu*s^2;
g GZ32= (We_gz/(y*s)/Xp)*seluu*s^2;
*
id tuM^-2*seluu= 1;
id seluu= 0;
*
id tuM^-2*selu= 0;
id tuM^-1*selu= 1;
id selu= 0;
*
id tuM^-2*sel0= 0;
id tuM^-1*sel0= 0;
*
id sel0= 1;
*
ratio,X,Xp,one;
id one= 1;
*--
id Re(wg(n?)*wzC(2))= 1/2*(Re(wg(n)*wzC([2+3]))+Re(wg(n)*wzC([2-3])));
id Re(wg(n?)*wzC(3))= 1/2*(Re(wg(n)*wzC([2+3]))-Re(wg(n)*wzC([2-3])));
id Re(wz(n?)*wgC(2))= 1/2*(Re(wz(n)*wgC([2+3]))+Re(wz(n)*wgC([2-3])));
id Re(wz(n?)*wgC(3))= 1/2*(Re(wz(n)*wgC([2+3]))-Re(wz(n)*wgC([2-3])));
*
id Re(wg(2)*wzC(n?))= 1/2*(Re(wg([2+3])*wzC(n))+Re(wg([2-3])*wzC(n)));
id Re(wg(3)*wzC(n?))= 1/2*(Re(wg([2+3])*wzC(n))-Re(wg([2-3])*wzC(n)));
id Re(wz(2)*wgC(n?))= 1/2*(Re(wz([2+3])*wgC(n))+Re(wz([2-3])*wgC(n)));
id Re(wz(3)*wgC(n?))= 1/2*(Re(wz([2+3])*wgC(n))-Re(wz([2-3])*wgC(n)));
*
id Re(wg(n?)*wzC(5))= 1/2*(Re(wg(n)*wzC([5+6]))+Re(wg(n)*wzC([5-6])));
id Re(wg(n?)*wzC(6))= 1/2*(Re(wg(n)*wzC([5+6]))-Re(wg(n)*wzC([5-6])));
id Re(wz(n?)*wgC(5))= 1/2*(Re(wz(n)*wgC([5+6]))+Re(wz(n)*wgC([5-6])));
id Re(wz(n?)*wgC(6))= 1/2*(Re(wz(n)*wgC([5+6]))-Re(wz(n)*wgC([5-6])));
*
id Re(wg(5)*wzC(n?))= 1/2*(Re(wg([5+6])*wzC(n))+Re(wg([5-6])*wzC(n)));
id Re(wg(6)*wzC(n?))= 1/2*(Re(wg([5+6])*wzC(n))-Re(wg([5-6])*wzC(n)));
id Re(wz(5)*wgC(n?))= 1/2*(Re(wz([5+6])*wgC(n))+Re(wz([5-6])*wgC(n)));
id Re(wz(6)*wgC(n?))= 1/2*(Re(wz([5+6])*wgC(n))-Re(wz([5-6])*wgC(n)));
*
id Re(wg(2)*X?*Y?)= 1/2*(Re(wg([2+3])*X*Y)+Re(wg([2-3])*X*Y));
id Re(wg(3)*X?*Y?)= 1/2*(Re(wg([2+3])*X*Y)-Re(wg([2-3])*X*Y));
id Re(wg(5)*X?*Y?)= 1/2*(Re(wg([5+6])*X*Y)+Re(wg([5-6])*X*Y));
id Re(wg(6)*X?*Y?)= 1/2*(Re(wg([5+6])*X*Y)-Re(wg([5-6])*X*Y));
id Re(wz(2)*X?*Y?)= 1/2*(Re(wz([2+3])*X*Y)+Re(wz([2-3])*X*Y));
id Re(wz(3)*X?*Y?)= 1/2*(Re(wz([2+3])*X*Y)-Re(wz([2-3])*X*Y));
id Re(wz(5)*X?*Y?)= 1/2*(Re(wz([5+6])*X*Y)+Re(wz([5-6])*X*Y));
id Re(wz(6)*X?*Y?)= 1/2*(Re(wz([5+6])*X*Y)-Re(wz([5-6])*X*Y));
*
id Re(wg(2)*X?)= 1/2*(Re(wg([2+3])*X)+Re(wg([2-3])*X));
id Re(wg(3)*X?)= 1/2*(Re(wg([2+3])*X)-Re(wg([2-3])*X));
id Re(wg(5)*X?)= 1/2*(Re(wg([5+6])*X)+Re(wg([5-6])*X));
id Re(wg(6)*X?)= 1/2*(Re(wg([5+6])*X)-Re(wg([5-6])*X));
id Re(wz(2)*X?)= 1/2*(Re(wz([2+3])*X)+Re(wz([2-3])*X));
id Re(wz(3)*X?)= 1/2*(Re(wz([2+3])*X)-Re(wz([2-3])*X));
id Re(wz(5)*X?)= 1/2*(Re(wz([5+6])*X)+Re(wz([5-6])*X));
id Re(wz(6)*X?)= 1/2*(Re(wz([5+6])*X)-Re(wz([5-6])*X));
*
id Re(wg([2+3])*wzC([2-3]))= Re(wz([2-3])*wgC([2+3]));
id Re(wg([2+3])*wzC([5-6]))= Re(wz([5-6])*wgC([2+3]));
id Re(wg([2-3])*wzC([5+6]))= Re(wz([5+6])*wgC([2-3]));
id Re(wg([5+6])*wzC([5-6]))= Re(wz([5-6])*wgC([5+6]));
id Re(wz([2+3])*wgC([2-3]))= Re(wg([2-3])*wzC([2+3]));
id Re(wz([2+3])*wgC([5-6]))= Re(wg([5-6])*wzC([2+3]));
id Re(wz([2-3])*wgC([5+6]))= Re(wg([5+6])*wzC([2-3]));
id Re(wz([5+6])*wgC([5-6]))= Re(wg([5-6])*wzC([5+6]));
*
id X= [1+X]-1;
id x1= 1-[1-x1];
id x2= 1-[1-x1]-[x1-x2];
id z= 1-[1-x1]-[x1-x2]-[x2-z];
id tau= [tau+X+x1]-[1+X]+[1-x1];
id Xp^-1= 1/[1+X];
ratio,X,[1+X],one;
id one= 1;
al one^-1= 1;
*
id Re(wg(X?)*wzC(Y?))= Rep(X,Y);
id Re(wz(X?)*wgC(Y?))= Rep(Y,X);
id GsM2*GtM2= GstM2;
id [1+X] = Xp;
al [1+X]^-1= Xpi;
id [1-x1]= omx1;
al [x1-x2]= x12;
al [x2-z]= x2z;
al [tau+X+x1]= tx;
al zX^-1= zXi;
id Xpi^2= Xpi-aux;
id aux= XX;
*
id Xp^2= Xps;
id Xpi^2= Xpis;
id y^3= yc;
id y^2= ys;
*
id omx1= O;
id O^4= Oq;
id O^3= Oc;
id O^2= Os;
id x12^3= x12c;
id x12^2= x12s;
id x2z^4= x2zq;
id x2z^3= x2zc;
id x2z^2= x2zs;
id tx^2= txs;
id vv^2= vvs;
id vv^3= vvc;
id vvs*vv= vvc;
id vvs^2= vvq;
id vvc*vv= vvq;
id z^-1= zi;
*
id Re(Gs*GtC*Pwt*PwsC)= rl1;
id Re(Gs*GtC*Pwt)= rl2;
id Re(Pws)= rl3;
id Re(wg(n?)*Gs*PwsC)= rs1wg(n);
id Re(wz(n?)*Gs*PwsC)= rs1wz(n);
id Re(wg(n?)*Gs)= rs2wg(n);
id Re(wz(n?)*Gs)= rs2wz(n);
id Re(wg(n?)*Gt*PwtC)= rt1wg(n);
id Re(wz(n?)*Gt*PwtC)= rt1wz(n);
id Re(wg(n?)*Gt)= rt2wg(n);
id Re(wz(n?)*Gt)= rt2wz(n);
id Rep(n1?,n2?)= rwwgz(n1,n2);
id Re(wg(n1?)*wzC(n2?))= rwwgz(n1,n2);
id Re(wz(n1?)*wgC(n2?))= rwwgz(n2,n1);
*
id fun?([2+3],n?)= fun(11,n);
al fun?([2-3],n?)= fun(12,n);
al fun?([5+6],n?)= fun(13,n);
al fun?([5-6],n?)= fun(14,n);
*
id fun?(n?,[2+3])= fun(n,11);
al fun?(n?,[2-3])= fun(n,12);
al fun?(n?,[5+6])= fun(n,13);
al fun?(n?,[5-6])= fun(n,14);
*
id fun?([2+3])= fun(11);
al fun?([2-3])= fun(12);
al fun?([5+6])= fun(13);
al fun?([5-6])= fun(14);
*--
id ve(T)= -1/2+ss(T);
al vu(T)= 1/2-2/3*ss(T);
al vd(T)= -1/2+1/3*ss(T);
*
id ss(T)*csi(T)= csi(T)-1;
id csi(T)= 1/cst;
al ss(T)= sst;
*
.store
*
*----g^2--------------------
*
g GZ10s= GZ10;
g GZ20s= GZ20;
g GZ30s= GZ30;
g GZ11s= GZ11;
g GZ21s= GZ21;
g GZ31s= GZ31;
g GZ12s= GZ12;
g GZ22s= GZ22;
g GZ32s= GZ32;
*
id gexp^4= 0;
id gexp^3= 0;
id gexp^2= 1;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,rs1wg,rs2wg,rt1wg,rt2wg,rwwg,
  rs1wz,rs2wz,rt1wz,rt2wz,rwwz,rwwgz,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
  zXi,XX,vv,vvs,vvc,vvq,cst,sst,PZ;
format fortran;
print;
*
.store
*
*----g^3--------------------
*
g GZ10c= GZ10;
g GZ20c= GZ20;
g GZ30c= GZ30;
g GZ11c= GZ11;
g GZ21c= GZ21;
g GZ31c= GZ31;
g GZ12c= GZ12;
g GZ22c= GZ22;
g GZ32c= GZ32;
*
id gexp^4= 0;
id gexp^3= 1;
id gexp^2= 0;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,rs1wg,rs2wg,rt1wg,rt2wg,rwwg,
  rs1wz,rs2wz,rt1wz,rt2wz,rwwz,rwwgz,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
  zXi,XX,vv,vvs,vvc,vvq,cst,sst,PZ;
format fortran;
print;
*
*
.store
*
*----g^4--------------------
*
g GZ10q= GZ10;
g GZ20q= GZ20;
g GZ30q= GZ30;
g GZ11q= GZ11;
g GZ21q= GZ21;
g GZ31q= GZ31;
g GZ12q= GZ12;
g GZ22q= GZ22;
g GZ32q= GZ32;
*
id gexp^4= 1;
id gexp^3= 0;
id gexp^2= 0;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,rs1wg,rs2wg,rt1wg,rt2wg,rwwg,
  rs1wz,rs2wz,rt1wz,rt2wz,rwwz,rwwgz,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
  zXi,XX,vv,vvs,vvc,vvq,cst,sst,PZ;
format fortran;
print;
*
#endprocedure
*
#procedure totalzze(dummy)
*
* Now the full stuff
*
load wrl$grp04:[giampiero]outputampzz_A.sav;
*
g T_zz= PwsM2*PwtM2*W_zz_int6/(GT^2*I0);
*nwrite statistics;
*
id e_(a1?,a2?,a3?,a4?)*sel0= e_(a1,a2,a3,a4);
id sel0= 0;
*
id wz(n?)= wz(n)*rs;
al wzC(n?)= wzC(n)*rs;
*
id wz(4)= X*wz(4);
al wz(7)= X*wz(7);
al wz(8)= X*wz(8);
al wz(10)= X*wz(10);
id wzC(4)= X*wzC(4);
al wzC(7)= X*wzC(7);
al wzC(8)= X*wzC(8);
al wzC(10)= X*wzC(10);
*
id wz(n?)*wzC(n?)= wzM2(n);
*
id pp.pp= 0;
al qp.qp= 0;
al Qm.Qm= QmS;
id Qp.Qm= pp.Qm-qp.Qm;
id kb.Qp= kb.pp-kb.qp;
id k.Qp= k.pp-k.qp;
id pp.Qp= -pp.qp;
id Qp.Qp= -2*pp.qp;
id qp.Qp= qp.pp;
*
id pp.Qm= -1/2*(QmS+y*s);
al k.kb= 1/2*(-sp+ums+dms);
al Qm.k= 1/2*(QmS+t-ums);
al pp.kb= 1/2*(tp-dms);
al Qm.kb= 1/2*(QmS+u-dms);
al pp.k= 1/2*(up-ums);
al pp.qp= 1/2*kp;
al Qm.qp= 1/2*(QmS+km);
al kb.qp= 1/2*(dms-zp);
al k.qp= 1/2*(ums-zm);
*
id t= tau*y*s;
al sp= x2*y*s;
al km= z*y*s;
al km^-1= 1/(z*y*s);
al zm= -(x2-x1)*y*s;
id QmS= X*y*s;
*
ratio,tuM,tau,aux;
id aux= ums/(y*s);
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
id kp= y*s * (  - 1 - X + x2 - z );
al u= y*s * (  - 1 - tau - 2*X - z )
       + ums + dms;
al up= y*s * (  - tau - X - x1 )
       + 2*ums + dms;
al tp= y*s * ( tau + X + x1 - x2 + z )
       - ums;
al zp= y*s * ( 1 - x1 )
       + ums + dms;
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
id GtM2*GsM2*wz(2)*wzC(1)= GtM2*GsM2*(-wzC(2)*wz(1)+
                                 2*Re(wz(2)*wzC(1)));
al GtM2*GsM2*wz(3)*wzC(1)= GtM2*GsM2*(-wzC(3)*wz(1)+
                                 2*Re(wz(3)*wzC(1)));
al GtM2*GsM2*wz(4)*wzC(1)= GtM2*GsM2*(-wzC(4)*wz(1)+
                                 2*Re(wz(4)*wzC(1)));
al GtM2*GsM2*wz(5)*wzC(1)= GtM2*GsM2*(-wzC(5)*wz(1)+
                                 2*Re(wz(5)*wzC(1)));
al GtM2*GsM2*wz(6)*wzC(1)= GtM2*GsM2*(-wzC(6)*wz(1)+
                                 2*Re(wz(6)*wzC(1)));
al GtM2*GsM2*wz(7)*wzC(1)= GtM2*GsM2*(-wzC(7)*wz(1)+
                                 2*Re(wz(7)*wzC(1)));
al GtM2*GsM2*wz(8)*wzC(1)= GtM2*GsM2*(-wzC(8)*wz(1)+
                                 2*Re(wz(8)*wzC(1)));
al GtM2*GsM2*wz(9)*wzC(1)= GtM2*GsM2*(-wzC(9)*wz(1)+
                                 2*Re(wz(9)*wzC(1)));
al GtM2*GsM2*wz(10)*wzC(1)= GtM2*GsM2*(-wzC(10)*wz(1)+
                                 2*Re(wz(10)*wzC(1)));
*
al GtM2*GsM2*wz(3)*wzC(2)= GtM2*GsM2*(-wzC(3)*wz(2)+
                                 2*Re(wz(3)*wzC(2)));
al GtM2*GsM2*wz(4)*wzC(2)= GtM2*GsM2*(-wzC(4)*wz(2)+
                                 2*Re(wz(4)*wzC(2)));
al GtM2*GsM2*wz(5)*wzC(2)= GtM2*GsM2*(-wzC(5)*wz(2)+
                                 2*Re(wz(5)*wzC(2)));
al GtM2*GsM2*wz(6)*wzC(2)= GtM2*GsM2*(-wzC(6)*wz(2)+
                                 2*Re(wz(6)*wzC(2)));
al GtM2*GsM2*wz(7)*wzC(2)= GtM2*GsM2*(-wzC(7)*wz(2)+
                                 2*Re(wz(7)*wzC(2)));
al GtM2*GsM2*wz(8)*wzC(2)= GtM2*GsM2*(-wzC(8)*wz(2)+
                                 2*Re(wz(8)*wzC(2)));
al GtM2*GsM2*wz(9)*wzC(2)= GtM2*GsM2*(-wzC(9)*wz(2)+
                                 2*Re(wz(9)*wzC(2)));
al GtM2*GsM2*wz(10)*wzC(2)= GtM2*GsM2*(-wzC(10)*wz(2)+
                                 2*Re(wz(10)*wzC(2)));
*
al GtM2*GsM2*wz(4)*wzC(3)= GtM2*GsM2*(-wzC(4)*wz(3)+
                                 2*Re(wz(4)*wzC(3)));
al GtM2*GsM2*wz(5)*wzC(3)= GtM2*GsM2*(-wzC(5)*wz(3)+
                                 2*Re(wz(5)*wzC(3)));
al GtM2*GsM2*wz(6)*wzC(3)= GtM2*GsM2*(-wzC(6)*wz(3)+
                                 2*Re(wz(6)*wzC(3)));
al GtM2*GsM2*wz(7)*wzC(3)= GtM2*GsM2*(-wzC(7)*wz(3)+
                                 2*Re(wz(7)*wzC(3)));
al GtM2*GsM2*wz(8)*wzC(3)= GtM2*GsM2*(-wzC(8)*wz(3)+
                                 2*Re(wz(8)*wzC(3)));
al GtM2*GsM2*wz(9)*wzC(3)= GtM2*GsM2*(-wzC(9)*wz(3)+
                                 2*Re(wz(9)*wzC(3)));
al GtM2*GsM2*wz(10)*wzC(3)= GtM2*GsM2*(-wzC(10)*wz(3)+
                                 2*Re(wz(10)*wzC(3)));
*
al GtM2*GsM2*wz(5)*wzC(4)= GtM2*GsM2*(-wzC(5)*wz(4)+
                                 2*Re(wz(5)*wzC(4)));
al GtM2*GsM2*wz(6)*wzC(4)= GtM2*GsM2*(-wzC(6)*wz(4)+
                                 2*Re(wz(6)*wzC(4)));
al GtM2*GsM2*wz(7)*wzC(4)= GtM2*GsM2*(-wzC(7)*wz(4)+
                                 2*Re(wz(7)*wzC(4)));
al GtM2*GsM2*wz(8)*wzC(4)= GtM2*GsM2*(-wzC(8)*wz(4)+
                                 2*Re(wz(8)*wzC(4)));
al GtM2*GsM2*wz(9)*wzC(4)= GtM2*GsM2*(-wzC(9)*wz(4)+
                                 2*Re(wz(9)*wzC(4)));
al GtM2*GsM2*wz(10)*wzC(4)= GtM2*GsM2*(-wzC(10)*wz(4)+
                                 2*Re(wz(10)*wzC(4)));
*
al GtM2*GsM2*wz(6)*wzC(5)= GtM2*GsM2*(-wzC(6)*wz(5)+
                                 2*Re(wz(6)*wzC(5)));
al GtM2*GsM2*wz(7)*wzC(5)= GtM2*GsM2*(-wzC(7)*wz(5)+
                                 2*Re(wz(7)*wzC(5)));
al GtM2*GsM2*wz(8)*wzC(5)= GtM2*GsM2*(-wzC(8)*wz(5)+
                                 2*Re(wz(8)*wzC(5)));
al GtM2*GsM2*wz(9)*wzC(5)= GtM2*GsM2*(-wzC(9)*wz(5)+
                                 2*Re(wz(9)*wzC(5)));
al GtM2*GsM2*wz(10)*wzC(5)= GtM2*GsM2*(-wzC(10)*wz(5)+
                                 2*Re(wz(10)*wzC(5)));
*
al GtM2*GsM2*wz(7)*wzC(6)= GtM2*GsM2*(-wzC(7)*wz(6)+
                                 2*Re(wz(7)*wzC(6)));
al GtM2*GsM2*wz(8)*wzC(6)= GtM2*GsM2*(-wzC(8)*wz(6)+
                                 2*Re(wz(8)*wzC(6)));
al GtM2*GsM2*wz(9)*wzC(6)= GtM2*GsM2*(-wzC(9)*wz(6)+
                                 2*Re(wz(9)*wzC(6)));
al GtM2*GsM2*wz(10)*wzC(6)= GtM2*GsM2*(-wzC(10)*wz(6)+
                                 2*Re(wz(10)*wzC(6)));
*
al GtM2*GsM2*wz(8)*wzC(7)= GtM2*GsM2*(-wzC(8)*wz(7)+
                                 2*Re(wz(8)*wzC(7)));
al GtM2*GsM2*wz(9)*wzC(7)= GtM2*GsM2*(-wzC(9)*wz(7)+
                                 2*Re(wz(9)*wzC(7)));
al GtM2*GsM2*wz(10)*wzC(7)= GtM2*GsM2*(-wzC(10)*wz(7)+
                                 2*Re(wz(10)*wzC(7)));
*
al GtM2*GsM2*wz(9)*wzC(8)= GtM2*GsM2*(-wzC(9)*wz(8)+
                                 2*Re(wz(9)*wzC(8)));
al GtM2*GsM2*wz(10)*wzC(8)= GtM2*GsM2*(-wzC(10)*wz(8)+
                                 2*Re(wz(10)*wzC(8)));
*
al GtM2*GsM2*wz(10)*wzC(9)= GtM2*GsM2*(-wzC(10)*wz(9)+
                                 2*Re(wz(10)*wzC(9)));
*
id wz(n?)*Gs*GtM2*PwsC= GtM2*(-wzC(n)*GsC*Pws+2*Re(wz(n)*Gs*PwsC));
id wz(n?)*Gt*GsM2*PwtC= GsM2*(-wzC(n)*GtC*Pwt+2*Re(wz(n)*Gt*PwtC));
id wz(n?)*Gs*GtM2= GtM2*(-GsC*wzC(n)+2*Re(Gs*wz(n)));
id wz(n?)*Gt*GsM2= GsM2*(-GtC*wzC(n)+2*Re(Gt*wz(n)));
*
id GsM2*GtC*PwtC*wzC(n?)= GsM2*(-Gt*Pwt*wz(n)+2*Re(Gt*Pwt*wz(n)));
id GtM2*GsC*PwsC*wzC(n?)= GtM2*(-Gs*Pws*wz(n)+2*Re(Gs*Pws*wz(n)));
*
id GsM2*GtC*PwtC= GsM2*(-Gt*Pwt+2*Re(Gt*Pwt));
id GtM2*GsC*PwsC= GtM2*(-Gs*Pws+2*Re(Gs*Pws));
*
id GsC*PwsC*Gt*Pwt= -Gs*Pws*GtC*PwtC+2*Re(GsC*Gt*PwsC*Pwt);
*id GsC*Gt*Pws*PwtC= -Gs*GtC*PwsC*Pwt+2*Re(Gs*GtC*PwsC*Pwt);
id GsC*Gt*Pws*PwtC= Re(Gs*GtC*PwsC*Pwt)+i_*Im(Gs*GtC*PwsC*Pwt);
id Gs*GtC*PwsC*Pwt= Re(Gs*GtC*PwsC*Pwt)-i_*Im(Gs*GtC*PwsC*Pwt);
id GsC*Gt*PwsC*PwtC= -Gs*GtC*Pws*Pwt+2*Re(Gs*GtC*Pws*Pwt);
*
id Gs*GtC*Pwt= -GsC*Gt*PwtC+2*Re(Gs*GtC*Pwt);
*
id GtM2*PwsC= GtM2*(-Pws+2*Re(Pws));
id GsM2*PwtC= GsM2*(-Pwt+2*Re(Pwt));
*
id Re(wzC(2)*wz(1))= Re(wz(2)*wzC(1));
al Re(wzC(3)*wz(1))= Re(wz(3)*wzC(1));
al Re(wzC(4)*wz(1))= Re(wz(4)*wzC(1));
al Re(wzC(5)*wz(1))= Re(wz(5)*wzC(1));
al Re(wzC(6)*wz(1))= Re(wz(6)*wzC(1));
al Re(wzC(7)*wz(1))= Re(wz(7)*wzC(1));
al Re(wzC(8)*wz(1))= Re(wz(8)*wzC(1));
al Re(wzC(9)*wz(1))= Re(wz(9)*wzC(1));
al Re(wzC(10)*wz(1))= Re(wz(10)*wzC(1));
*
id Re(wzC(3)*wz(2))= Re(wz(3)*wzC(2));
al Re(wzC(4)*wz(2))= Re(wz(4)*wzC(2));
al Re(wzC(5)*wz(2))= Re(wz(5)*wzC(2));
al Re(wzC(6)*wz(2))= Re(wz(6)*wzC(2));
al Re(wzC(7)*wz(2))= Re(wz(7)*wzC(2));
al Re(wzC(8)*wz(2))= Re(wz(8)*wzC(2));
al Re(wzC(9)*wz(2))= Re(wz(9)*wzC(2));
al Re(wzC(10)*wz(2))= Re(wz(10)*wzC(2));
*
id Re(wzC(4)*wz(3))= Re(wz(4)*wzC(3));
al Re(wzC(5)*wz(3))= Re(wz(5)*wzC(3));
al Re(wzC(6)*wz(3))= Re(wz(6)*wzC(3));
al Re(wzC(7)*wz(3))= Re(wz(7)*wzC(3));
al Re(wzC(8)*wz(3))= Re(wz(8)*wzC(3));
al Re(wzC(9)*wz(3))= Re(wz(9)*wzC(3));
al Re(wzC(10)*wz(3))= Re(wz(10)*wzC(3));
*
id Re(wzC(5)*wz(4))= Re(wz(5)*wzC(4));
al Re(wzC(6)*wz(4))= Re(wz(6)*wzC(4));
al Re(wzC(7)*wz(4))= Re(wz(7)*wzC(4));
al Re(wzC(8)*wz(4))= Re(wz(8)*wzC(4));
al Re(wzC(9)*wz(4))= Re(wz(9)*wzC(4));
al Re(wzC(10)*wz(4))= Re(wz(10)*wzC(4));
*
id Re(wzC(6)*wz(5))= Re(wz(6)*wzC(5));
al Re(wzC(7)*wz(5))= Re(wz(7)*wzC(5));
al Re(wzC(8)*wz(5))= Re(wz(8)*wzC(5));
al Re(wzC(9)*wz(5))= Re(wz(9)*wzC(5));
al Re(wzC(10)*wz(5))= Re(wz(10)*wzC(5));
*
id Re(wzC(7)*wz(6))= Re(wz(7)*wzC(6));
al Re(wzC(8)*wz(6))= Re(wz(8)*wzC(6));
al Re(wzC(9)*wz(6))= Re(wz(9)*wzC(6));
al Re(wzC(10)*wz(6))= Re(wz(10)*wzC(6));
*
id Re(wzC(8)*wz(7))= Re(wz(8)*wzC(7));
al Re(wzC(9)*wz(7))= Re(wz(9)*wzC(7));
al Re(wzC(10)*wz(7))= Re(wz(10)*wzC(7));
*
id Re(wzC(9)*wz(8))= Re(wz(9)*wzC(8));
al Re(wzC(10)*wz(8))= Re(wz(10)*wzC(8));
*
id Re(wzC(10)*wz(9))= Re(wz(10)*wzC(9));
*
id vv^2= vvs;
id vv^3= vvc;
id vvs*vv= vvc;
id vvs^2= vvq;
id vvc*vv= vvq;
*
b i_,pp,Qm,d_,e_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,ve,vd,vu,ae,ad,au,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,X,Xp,Xps,Xpi,Xpis,
  zXi,XX,vv,vvs,vvc,vvq,cst,sst,PZ,csi,gexp,wz,wzC,wzM2,Re,Im,e_;
print;
*
#endprocedure
*
#procedure totalzz(dummy)
*
* Now the full stuff
*
load wrl$grp04:[giampiero]outputampzz_A.sav;
*
g T_zz= PwsM2*PwtM2*W_zz_int6/(GT^2*I0);
*nwrite statistics;
*
id wz(n?)= wz(n)*rs;
al wzC(n?)= wzC(n)*rs;
*
id wz(4)= X*wz(4);
al wz(7)= X*wz(7);
al wz(8)= X*wz(8);
al wz(10)= X*wz(10);
id wzC(4)= X*wzC(4);
al wzC(7)= X*wzC(7);
al wzC(8)= X*wzC(8);
al wzC(10)= X*wzC(10);
*
id wz(n?)*wzC(n?)= wzM2(n);
*
id GtM2*GsM2*wz(2)*wzC(1)= GtM2*GsM2*(-wzC(2)*wz(1)+
                                 2*Re(wz(2)*wzC(1)));
al GtM2*GsM2*wz(3)*wzC(1)= GtM2*GsM2*(-wzC(3)*wz(1)+
                                 2*Re(wz(3)*wzC(1)));
al GtM2*GsM2*wz(4)*wzC(1)= GtM2*GsM2*(-wzC(4)*wz(1)+
                                 2*Re(wz(4)*wzC(1)));
al GtM2*GsM2*wz(5)*wzC(1)= GtM2*GsM2*(-wzC(5)*wz(1)+
                                 2*Re(wz(5)*wzC(1)));
al GtM2*GsM2*wz(6)*wzC(1)= GtM2*GsM2*(-wzC(6)*wz(1)+
                                 2*Re(wz(6)*wzC(1)));
al GtM2*GsM2*wz(7)*wzC(1)= GtM2*GsM2*(-wzC(7)*wz(1)+
                                 2*Re(wz(7)*wzC(1)));
al GtM2*GsM2*wz(8)*wzC(1)= GtM2*GsM2*(-wzC(8)*wz(1)+
                                 2*Re(wz(8)*wzC(1)));
al GtM2*GsM2*wz(9)*wzC(1)= GtM2*GsM2*(-wzC(9)*wz(1)+
                                 2*Re(wz(9)*wzC(1)));
al GtM2*GsM2*wz(10)*wzC(1)= GtM2*GsM2*(-wzC(10)*wz(1)+
                                 2*Re(wz(10)*wzC(1)));
*
al GtM2*GsM2*wz(3)*wzC(2)= GtM2*GsM2*(-wzC(3)*wz(2)+
                                 2*Re(wz(3)*wzC(2)));
al GtM2*GsM2*wz(4)*wzC(2)= GtM2*GsM2*(-wzC(4)*wz(2)+
                                 2*Re(wz(4)*wzC(2)));
al GtM2*GsM2*wz(5)*wzC(2)= GtM2*GsM2*(-wzC(5)*wz(2)+
                                 2*Re(wz(5)*wzC(2)));
al GtM2*GsM2*wz(6)*wzC(2)= GtM2*GsM2*(-wzC(6)*wz(2)+
                                 2*Re(wz(6)*wzC(2)));
al GtM2*GsM2*wz(7)*wzC(2)= GtM2*GsM2*(-wzC(7)*wz(2)+
                                 2*Re(wz(7)*wzC(2)));
al GtM2*GsM2*wz(8)*wzC(2)= GtM2*GsM2*(-wzC(8)*wz(2)+
                                 2*Re(wz(8)*wzC(2)));
al GtM2*GsM2*wz(9)*wzC(2)= GtM2*GsM2*(-wzC(9)*wz(2)+
                                 2*Re(wz(9)*wzC(2)));
al GtM2*GsM2*wz(10)*wzC(2)= GtM2*GsM2*(-wzC(10)*wz(2)+
                                 2*Re(wz(10)*wzC(2)));
*
al GtM2*GsM2*wz(4)*wzC(3)= GtM2*GsM2*(-wzC(4)*wz(3)+
                                 2*Re(wz(4)*wzC(3)));
al GtM2*GsM2*wz(5)*wzC(3)= GtM2*GsM2*(-wzC(5)*wz(3)+
                                 2*Re(wz(5)*wzC(3)));
al GtM2*GsM2*wz(6)*wzC(3)= GtM2*GsM2*(-wzC(6)*wz(3)+
                                 2*Re(wz(6)*wzC(3)));
al GtM2*GsM2*wz(7)*wzC(3)= GtM2*GsM2*(-wzC(7)*wz(3)+
                                 2*Re(wz(7)*wzC(3)));
al GtM2*GsM2*wz(8)*wzC(3)= GtM2*GsM2*(-wzC(8)*wz(3)+
                                 2*Re(wz(8)*wzC(3)));
al GtM2*GsM2*wz(9)*wzC(3)= GtM2*GsM2*(-wzC(9)*wz(3)+
                                 2*Re(wz(9)*wzC(3)));
al GtM2*GsM2*wz(10)*wzC(3)= GtM2*GsM2*(-wzC(10)*wz(3)+
                                 2*Re(wz(10)*wzC(3)));
*
al GtM2*GsM2*wz(5)*wzC(4)= GtM2*GsM2*(-wzC(5)*wz(4)+
                                 2*Re(wz(5)*wzC(4)));
al GtM2*GsM2*wz(6)*wzC(4)= GtM2*GsM2*(-wzC(6)*wz(4)+
                                 2*Re(wz(6)*wzC(4)));
al GtM2*GsM2*wz(7)*wzC(4)= GtM2*GsM2*(-wzC(7)*wz(4)+
                                 2*Re(wz(7)*wzC(4)));
al GtM2*GsM2*wz(8)*wzC(4)= GtM2*GsM2*(-wzC(8)*wz(4)+
                                 2*Re(wz(8)*wzC(4)));
al GtM2*GsM2*wz(9)*wzC(4)= GtM2*GsM2*(-wzC(9)*wz(4)+
                                 2*Re(wz(9)*wzC(4)));
al GtM2*GsM2*wz(10)*wzC(4)= GtM2*GsM2*(-wzC(10)*wz(4)+
                                 2*Re(wz(10)*wzC(4)));
*
al GtM2*GsM2*wz(6)*wzC(5)= GtM2*GsM2*(-wzC(6)*wz(5)+
                                 2*Re(wz(6)*wzC(5)));
al GtM2*GsM2*wz(7)*wzC(5)= GtM2*GsM2*(-wzC(7)*wz(5)+
                                 2*Re(wz(7)*wzC(5)));
al GtM2*GsM2*wz(8)*wzC(5)= GtM2*GsM2*(-wzC(8)*wz(5)+
                                 2*Re(wz(8)*wzC(5)));
al GtM2*GsM2*wz(9)*wzC(5)= GtM2*GsM2*(-wzC(9)*wz(5)+
                                 2*Re(wz(9)*wzC(5)));
al GtM2*GsM2*wz(10)*wzC(5)= GtM2*GsM2*(-wzC(10)*wz(5)+
                                 2*Re(wz(10)*wzC(5)));
*
al GtM2*GsM2*wz(7)*wzC(6)= GtM2*GsM2*(-wzC(7)*wz(6)+
                                 2*Re(wz(7)*wzC(6)));
al GtM2*GsM2*wz(8)*wzC(6)= GtM2*GsM2*(-wzC(8)*wz(6)+
                                 2*Re(wz(8)*wzC(6)));
al GtM2*GsM2*wz(9)*wzC(6)= GtM2*GsM2*(-wzC(9)*wz(6)+
                                 2*Re(wz(9)*wzC(6)));
al GtM2*GsM2*wz(10)*wzC(6)= GtM2*GsM2*(-wzC(10)*wz(6)+
                                 2*Re(wz(10)*wzC(6)));
*
al GtM2*GsM2*wz(8)*wzC(7)= GtM2*GsM2*(-wzC(8)*wz(7)+
                                 2*Re(wz(8)*wzC(7)));
al GtM2*GsM2*wz(9)*wzC(7)= GtM2*GsM2*(-wzC(9)*wz(7)+
                                 2*Re(wz(9)*wzC(7)));
al GtM2*GsM2*wz(10)*wzC(7)= GtM2*GsM2*(-wzC(10)*wz(7)+
                                 2*Re(wz(10)*wzC(7)));
*
al GtM2*GsM2*wz(9)*wzC(8)= GtM2*GsM2*(-wzC(9)*wz(8)+
                                 2*Re(wz(9)*wzC(8)));
al GtM2*GsM2*wz(10)*wzC(8)= GtM2*GsM2*(-wzC(10)*wz(8)+
                                 2*Re(wz(10)*wzC(8)));
*
al GtM2*GsM2*wz(10)*wzC(9)= GtM2*GsM2*(-wzC(10)*wz(9)+
                                 2*Re(wz(10)*wzC(9)));
*
id wz(n?)*Gs*GtM2*PwsC= GtM2*(-wzC(n)*GsC*Pws+2*Re(wz(n)*Gs*PwsC));
id wz(n?)*Gt*GsM2*PwtC= GsM2*(-wzC(n)*GtC*Pwt+2*Re(wz(n)*Gt*PwtC));
id wz(n?)*Gs*GtM2= GtM2*(-GsC*wzC(n)+2*Re(Gs*wz(n)));
id wz(n?)*Gt*GsM2= GsM2*(-GtC*wzC(n)+2*Re(Gt*wz(n)));
*
id GsM2*GtC*PwtC*wzC(n?)= GsM2*(-Gt*Pwt*wz(n)+2*Re(Gt*Pwt*wz(n)));
id GtM2*GsC*PwsC*wzC(n?)= GtM2*(-Gs*Pws*wz(n)+2*Re(Gs*Pws*wz(n)));
*
id GsM2*GtC*PwtC= GsM2*(-Gt*Pwt+2*Re(Gt*Pwt));
id GtM2*GsC*PwsC= GtM2*(-Gs*Pws+2*Re(Gs*Pws));
*
id GsC*PwsC*Gt*Pwt= -Gs*Pws*GtC*PwtC+2*Re(GsC*Gt*PwsC*Pwt);
*id GsC*Gt*Pws*PwtC= -Gs*GtC*PwsC*Pwt+2*Re(Gs*GtC*PwsC*Pwt);
id GsC*Gt*Pws*PwtC= Re(Gs*GtC*PwsC*Pwt)+i_*Im(Gs*GtC*PwsC*Pwt);
id Gs*GtC*PwsC*Pwt= Re(Gs*GtC*PwsC*Pwt)-i_*Im(Gs*GtC*PwsC*Pwt);
id GsC*Gt*PwsC*PwtC= -Gs*GtC*Pws*Pwt+2*Re(Gs*GtC*Pws*Pwt);
*
id Gs*GtC*Pwt= -GsC*Gt*PwtC+2*Re(Gs*GtC*Pwt);
*
id GtM2*PwsC= GtM2*(-Pws+2*Re(Pws));
id GsM2*PwtC= GsM2*(-Pwt+2*Re(Pwt));
*
id Re(wzC(2)*wz(1))= Re(wz(2)*wzC(1));
al Re(wzC(3)*wz(1))= Re(wz(3)*wzC(1));
al Re(wzC(4)*wz(1))= Re(wz(4)*wzC(1));
al Re(wzC(5)*wz(1))= Re(wz(5)*wzC(1));
al Re(wzC(6)*wz(1))= Re(wz(6)*wzC(1));
al Re(wzC(7)*wz(1))= Re(wz(7)*wzC(1));
al Re(wzC(8)*wz(1))= Re(wz(8)*wzC(1));
al Re(wzC(9)*wz(1))= Re(wz(9)*wzC(1));
al Re(wzC(10)*wz(1))= Re(wz(10)*wzC(1));
*
id Re(wzC(3)*wz(2))= Re(wz(3)*wzC(2));
al Re(wzC(4)*wz(2))= Re(wz(4)*wzC(2));
al Re(wzC(5)*wz(2))= Re(wz(5)*wzC(2));
al Re(wzC(6)*wz(2))= Re(wz(6)*wzC(2));
al Re(wzC(7)*wz(2))= Re(wz(7)*wzC(2));
al Re(wzC(8)*wz(2))= Re(wz(8)*wzC(2));
al Re(wzC(9)*wz(2))= Re(wz(9)*wzC(2));
al Re(wzC(10)*wz(2))= Re(wz(10)*wzC(2));
*
id Re(wzC(4)*wz(3))= Re(wz(4)*wzC(3));
al Re(wzC(5)*wz(3))= Re(wz(5)*wzC(3));
al Re(wzC(6)*wz(3))= Re(wz(6)*wzC(3));
al Re(wzC(7)*wz(3))= Re(wz(7)*wzC(3));
al Re(wzC(8)*wz(3))= Re(wz(8)*wzC(3));
al Re(wzC(9)*wz(3))= Re(wz(9)*wzC(3));
al Re(wzC(10)*wz(3))= Re(wz(10)*wzC(3));
*
id Re(wzC(5)*wz(4))= Re(wz(5)*wzC(4));
al Re(wzC(6)*wz(4))= Re(wz(6)*wzC(4));
al Re(wzC(7)*wz(4))= Re(wz(7)*wzC(4));
al Re(wzC(8)*wz(4))= Re(wz(8)*wzC(4));
al Re(wzC(9)*wz(4))= Re(wz(9)*wzC(4));
al Re(wzC(10)*wz(4))= Re(wz(10)*wzC(4));
*
id Re(wzC(6)*wz(5))= Re(wz(6)*wzC(5));
al Re(wzC(7)*wz(5))= Re(wz(7)*wzC(5));
al Re(wzC(8)*wz(5))= Re(wz(8)*wzC(5));
al Re(wzC(9)*wz(5))= Re(wz(9)*wzC(5));
al Re(wzC(10)*wz(5))= Re(wz(10)*wzC(5));
*
id Re(wzC(7)*wz(6))= Re(wz(7)*wzC(6));
al Re(wzC(8)*wz(6))= Re(wz(8)*wzC(6));
al Re(wzC(9)*wz(6))= Re(wz(9)*wzC(6));
al Re(wzC(10)*wz(6))= Re(wz(10)*wzC(6));
*
id Re(wzC(8)*wz(7))= Re(wz(8)*wzC(7));
al Re(wzC(9)*wz(7))= Re(wz(9)*wzC(7));
al Re(wzC(10)*wz(7))= Re(wz(10)*wzC(7));
*
id Re(wzC(9)*wz(8))= Re(wz(9)*wzC(8));
al Re(wzC(10)*wz(8))= Re(wz(10)*wzC(8));
*
id Re(wzC(10)*wz(9))= Re(wz(10)*wzC(9));
*
.store
*
#endprocedure
*
#procedure contractzz(dummy)
*
g Wd_zz= T_zz*d_(mu,nu);
g Wpp_zz= T_zz*pp(mu)*pp(nu);
g Wqq_zz= T_zz*Qm(mu)*Qm(nu);
g Wpq_zz= T_zz*(pp(mu)*Qm(nu)+pp(nu)*Qm(mu));
g We_zz= T_zz*e_(mu,nu,Qm,pp);
*
contract;
*
id rs^-1= -1/(X*y*s);
id pp.pp= 0;
al qp.qp= 0;
al Qm.Qm= QmS;
id Qp.Qm= pp.Qm-qp.Qm;
id kb.Qp= kb.pp-kb.qp;
id k.Qp= k.pp-k.qp;
id pp.Qp= -pp.qp;
id Qp.Qp= -2*pp.qp;
id qp.Qp= qp.pp;
*
id pp.Qm= -1/2*(QmS+y*s);
al k.kb= 1/2*(-sp+ums+dms);
al Qm.k= 1/2*(QmS+t-ums);
al pp.kb= 1/2*(tp-dms);
al Qm.kb= 1/2*(QmS+u-dms);
al pp.k= 1/2*(up-ums);
al pp.qp= 1/2*kp;
al Qm.qp= 1/2*(QmS+km);
al kb.qp= 1/2*(dms-zp);
al k.qp= 1/2*(ums-zm);
*
id t= tau*y*s;
al sp= x2*y*s;
al km= z*y*s;
al km^-1= 1/(z*y*s);
al zm= -(x2-x1)*y*s;
id QmS= X*y*s;
*
ratio,tuM,tau,aux;
id aux= ums/(y*s);
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
id kp= y*s * (  - 1 - X + x2 - z );
al u= y*s * (  - 1 - tau - 2*X - z )
       + ums + dms;
al up= y*s * (  - tau - X - x1 )
       + 2*ums + dms;
al tp= y*s * ( tau + X + x1 - x2 + z )
       - ums;
al zp= y*s * ( 1 - x1 )
       + ums + dms;
*
id ums*tuM^-2= aux;
id ums= 0;
id aux= ums*tuM^-2;
id dms= 0;
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
.store
*
g ZZ10= (-1/2*Wd_zz-2*X/Xp^2*Wpp_zz/(y*s)-Wpq_zz/(y*s)/Xp)*sel0*s^3;
g ZZ20= (-1/2*Wd_zz-6*X/Xp^2*Wpp_zz/(y*s)-3*Wpq_zz/(y*s)/Xp-
         Wqq_zz/(y*s)/Xp)*sel0*s^3;
g ZZ30= (We_zz/(y*s)/Xp)*sel0*s^3;
g ZZ40= (-2/Xp*Wqq_zz/(y*s)-2*X/Xp^2*Wpq_zz/(y*s))*sel0*s^3;
*g W5_zz_0u= (-Wqq_zz/(y*s)/X-2/Xp*Wpq_zz/(y*s))*sel0*s^3;
g ZZ11= (-1/2*Wd_zz-2*X/Xp^2*Wpp_zz/(y*s)-Wpq_zz/(y*s)/Xp)*selu*s^3;
g ZZ21= (-1/2*Wd_zz-6*X/Xp^2*Wpp_zz/(y*s)-3*Wpq_zz/(y*s)/Xp-
         Wqq_zz/(y*s)/Xp)*selu*s^3;
g ZZ31= (We_zz/(y*s)/Xp)*selu*s^3;
g ZZ41= (-2/Xp*Wqq_zz/(y*s)-2*X/Xp^2*Wpq_zz/(y*s))*selu*s^3;
*g W5_zz_u= (-Wqq_zz/(y*s)/X-2/Xp*Wpq_zz/(y*s))*selu*s^3;
g ZZ12= (-1/2*Wd_zz-2*X/Xp^2*Wpp_zz/(y*s)-Wpq_zz/(y*s)/Xp)*seluu*s^3;
g ZZ22= (-1/2*Wd_zz-6*X/Xp^2*Wpp_zz/(y*s)-3*Wpq_zz/(y*s)/Xp-
         Wqq_zz/(y*s)/Xp)*seluu*s^3;
g ZZ32= (We_zz/(y*s)/Xp)*seluu*s^3;
g ZZ42= (-2/Xp*Wqq_zz/(y*s)-2*X/Xp^2*Wpq_zz/(y*s))*seluu*s^3;
*g W5_zz_uu= (-Wqq_zz/(y*s)/X-2/Xp*Wpq_zz/(y*s))*seluu*s^3;
*
id tuM^-2*seluu= 1;
id seluu= 0;
*
id tuM^-2*selu= 0;
id tuM^-1*selu= 1;
id selu= 0;
*
id tuM^-2*sel0= 0;
id tuM^-1*sel0= 0;
*
id sel0= 1;
*
ratio,X,Xp,one;
id one= 1;
al one^-1= 1;
*
id cs(T)*csi(T)= 1;
*
*--
id Re(wz(n?)*wzC(2))= 1/2*(Re(wz(n)*wzC([2+3]))+Re(wz(n)*wzC([2-3])));
id Re(wz(n?)*wzC(3))= 1/2*(Re(wz(n)*wzC([2+3]))-Re(wz(n)*wzC([2-3])));
id Re(wz(2)*wzC(2))= wzM2(2);
id Re(wz(3)*wzC(3))= wzM2(3);
*
id Re(wz(2)*wzC(n?))= 1/2*(Re(wz([2+3])*wzC(n))+Re(wz([2-3])*wzC(n)));
id Re(wz(3)*wzC(n?))= 1/2*(Re(wz([2+3])*wzC(n))-Re(wz([2-3])*wzC(n)));
id Re(wz(2)*wzC(2))= wzM2(2);
id Re(wz(3)*wzC(3))= wzM2(3);
*
id Re(wz(n?)*wzC(5))= 1/2*(Re(wz(n)*wzC([5+6]))+Re(wz(n)*wzC([5-6])));
id Re(wz(n?)*wzC(6))= 1/2*(Re(wz(n)*wzC([5+6]))-Re(wz(n)*wzC([5-6])));
id Re(wz(5)*wzC(5))= wzM2(5);
id Re(wz(6)*wzC(6))= wzM2(6);
*
id Re(wz(5)*wzC(n?))= 1/2*(Re(wz([5+6])*wzC(n))+Re(wz([5-6])*wzC(n)));
id Re(wz(6)*wzC(n?))= 1/2*(Re(wz([5+6])*wzC(n))-Re(wz([5-6])*wzC(n)));
id Re(wz(5)*wzC(5))= wzM2(5);
id Re(wz(6)*wzC(6))= wzM2(6);
*
id Re(wz(2)*X?*Y?)= 1/2*(Re(wz([2+3])*X*Y)+Re(wz([2-3])*X*Y));
id Re(wz(3)*X?*Y?)= 1/2*(Re(wz([2+3])*X*Y)-Re(wz([2-3])*X*Y));
id Re(wz(5)*X?*Y?)= 1/2*(Re(wz([5+6])*X*Y)+Re(wz([5-6])*X*Y));
id Re(wz(6)*X?*Y?)= 1/2*(Re(wz([5+6])*X*Y)-Re(wz([5-6])*X*Y));
*
id Re(wz(2)*X?)= 1/2*(Re(wz([2+3])*X)+Re(wz([2-3])*X));
id Re(wz(3)*X?)= 1/2*(Re(wz([2+3])*X)-Re(wz([2-3])*X));
id Re(wz(5)*X?)= 1/2*(Re(wz([5+6])*X)+Re(wz([5-6])*X));
id Re(wz(6)*X?)= 1/2*(Re(wz([5+6])*X)-Re(wz([5-6])*X));
*
id Re(wz([2+3])*wzC([2-3]))= Re(wz([2-3])*wzC([2+3]));
id Re(wz([2+3])*wzC([5-6]))= Re(wz([5-6])*wzC([2+3]));
id Re(wz([2-3])*wzC([5+6]))= Re(wz([5+6])*wzC([2-3]));
id Re(wz([5+6])*wzC([5-6]))= Re(wz([5-6])*wzC([5+6]));
*
id Re(wz([2+3])*wzC([2+3]))= wzM2(2)+wzM2(3)+2*Re(wz(2)*wzC(3));
id Re(wz([2-3])*wzC([2-3]))= wzM2(2)+wzM2(3)-2*Re(wz(2)*wzC(3));
id Re(wz([5+6])*wzC([5+6]))= wzM2(5)+wzM2(6)+2*Re(wz(5)*wzC(6));
id Re(wz([5-6])*wzC([5-6]))= wzM2(5)+wzM2(6)-2*Re(wz(5)*wzC(6));
*--
id X= [1+X]-1;
id x1= 1-[1-x1];
id x2= 1-[1-x1]-[x1-x2];
id z= 1-[1-x1]-[x1-x2]-[x2-z];
id tau= [tau+X+x1]-[1+X]+[1-x1];
id Xp^-1= 1/[1+X];
ratio,X,[1+X],one;
id one= 1;
al one^-1= 1;
*
id Re(wz(X?)*wzC(Y?))= Rep(X,Y);
id GsM2*GtM2= GstM2;
id [1+X] = Xp;
al [1+X]^-1= Xpi;
id [1-x1]= omx1;
al [x1-x2]= x12;
al [x2-z]= x2z;
al [tau+X+x1]= tx;
al zX^-1= zXi;
id Xpi^2= Xpi-aux;
id aux= XX;
*
id Xp^2= Xps;
id Xpi^2= Xpis;
id y^3= yc;
id y^2= ys;
id PZ^-2= 1/PZs;
*
id omx1= O;
id O^4= Oq;
id O^3= Oc;
id O^2= Os;
id x12^3= x12c;
id x12^2= x12s;
id x2z^4= x2zq;
id x2z^3= x2zc;
id x2z^2= x2zs;
id tx^2= txs;
id vv^2= vvs;
id vv^3= vvc;
id vvs*vv= vvc;
id vvs^2= vvq;
id vvc*vv= vvq;
id csi(T)= 1/cst;
id cst^-2= 1/cqt;
id ve(T)= vet;
al vd(T)= vdt;
al vu(T)= vut;
id ae(T)= aet;
al ad(T)= adt;
al au(T)= aut;
*
id vet^2= vets;
al vdt^2= vdts;
al vut^2= vuts;
id aet^2= aets;
al adt^2= adts;
al aut^2= auts;
*
id Re(Gs*GtC*Pwt*PwsC)= rl1;
id Re(Gs*GtC*Pwt)= rl2;
id Re(Pws)= rl3;
id Re(wz(n?)*Gs*PwsC)= rs1wz(n);
id Re(wz(n?)*Gs)= rs2wz(n);
id Re(wz(n?)*Gt*PwtC)= rt1wz(n);
id Re(wz(n?)*Gt)= rt2wz(n);
id Rep(n1?,n2?)= rwwz(n1,n2);
id Re(wz(n1?)*wzC(n2?))= rwwz(n1,n2);
id wzM2(n?)= rwwz(n,n);
*
id fun?([2+3],n?)= fun(11,n);
al fun?([2-3],n?)= fun(12,n);
al fun?([5+6],n?)= fun(13,n);
al fun?([5-6],n?)= fun(14,n);
*
id fun?(n?,[2+3])= fun(n,11);
al fun?(n?,[2-3])= fun(n,12);
al fun?(n?,[5+6])= fun(n,13);
al fun?(n?,[5-6])= fun(n,14);
*
id fun?([2+3])= fun(11);
al fun?([2-3])= fun(12);
al fun?([5+6])= fun(13);
al fun?([5-6])= fun(14);
*
.store
*
*----g^2--------------------
*
g ZZ10s= ZZ10;
g ZZ20s= ZZ20;
g ZZ30s= ZZ30;
g ZZ40s= ZZ40;
g ZZ11s= ZZ11;
g ZZ21s= ZZ21;
g ZZ31s= ZZ31;
g ZZ41s= ZZ41;
g ZZ12s= ZZ12;
g ZZ22s= ZZ22;
g ZZ32s= ZZ32;
g ZZ42s= ZZ32;
*
id gexp^4= 0;
id gexp^3= 0;
id gexp^2= 1;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep,PZs,
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wz,wzC,wzM2,rs1wz,rs2wz,
  rt1wz,rt2wz,rwwz,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
  zXi,XX,vv,vvs,vvc,vvq,
  cst,cqt,vet,vdt,vut,aet,adt,aut,vets,vdts,vuts,aets,adts,auts;
*
format fortran;
print;
*
.store
*
*----g^3--------------------
*
g ZZ10c= ZZ10;
g ZZ20c= ZZ20;
g ZZ30c= ZZ30;
g ZZ40c= ZZ40;
g ZZ11c= ZZ11;
g ZZ21c= ZZ21;
g ZZ31c= ZZ31;
g ZZ41c= ZZ41;
g ZZ12c= ZZ12;
g ZZ22c= ZZ22;
g ZZ32c= ZZ32;
g ZZ42c= ZZ42;
*
id gexp^4= 0;
id gexp^3= 1;
id gexp^2= 0;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep,PZs,
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wz,wzC,wzM2,rs1wz,rs2wz,
  rt1wz,rt2wz,rwwz,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
  zXi,XX,vv,vvs,vvc,vvq,
  cst,cqt,vet,vdt,vut,aet,adt,aut,vets,vdts,vuts,aets,adts,auts;
format fortran;
print;
*
.store
*
*----g^4--------------------
*
g ZZ10q= ZZ10/GstM2;
g ZZ20q= ZZ20/GstM2;
g ZZ30q= ZZ30/GstM2;
g ZZ40q= ZZ40/GstM2;
g ZZ11q= ZZ11/GstM2;
g ZZ21q= ZZ21/GstM2;
g ZZ31q= ZZ31/GstM2;
g ZZ41q= ZZ41/GstM2;
g ZZ12q= ZZ12/GstM2;
g ZZ22q= ZZ22/GstM2;
g ZZ32q= ZZ32/GstM2;
g ZZ42q= ZZ42/GstM2;
*
id gexp^4= 1;
id gexp^3= 0;
id gexp^2= 0;
*
b i_,Gs,GsC,Gt,GtC,GsM2,GtM2,GstM2,Re,Mods,Rep,PZs,
  zX,ums,dms,y,ys,yc,rs,s,tuM,tdM,wz,wzC,wzM2,rs1wz,rs2wz,
  rt1wz,rt2wz,rwwz,rl1,rl2,rl3,
  Pws,PwsC,Pwt,PwtC,PwsM2,PwtM2,Sigmas,SigmasC,SigmasM2,X,Xp,Xps,Xpi,Xpis,
  zXi,XX,vv,vvs,vvc,vvq,
  cst,cqt,vet,vdt,vut,aet,adt,aut,vets,vdts,vuts,aets,adts,auts;
format fortran;
print;
*
#endprocedure
*
*///////////////////////////////////////////////////////////////////////
*
*#call wtensor{gg}
*#call totalgg{dummy}
*#call contractgg{dummy}
*----------------------------------------------------------------
*#call wtensor{gz}
*#call totalgz{dummy}
*#call contractgz{dummy}
*----------------------------------------------------------------
#call wtensor{zz}
*#call totalzz{dummy}
*#call contractzz{dummy}
*#call totalzze{dummy}
*
.end