*
* Rechnernutzung in der Physik
* WS 12/13
*
* Gamma-Matrizen
*

index mu,nu,rho,si;
index mu1,...,mu20;
vector p,pp;
vector p1,p2,p3,p4;
symbol m,beta,E,[sin(the/2)];

local t1  = g_(1,mu)*g_(1,nu);
local t2  = g_(1,mu,nu,rho,si);
local t2a = g_(1,p1,p2,p3,p4);
***local t3  = g_(1,mu1,...,mu20);
local t4  = 1/24 * e_(mu,nu,rho,si) * g_(1,mu,nu,rho,si,p1,p2,p3,p4);
local t4a = g_(1,5_, p1,p2,p3,p4);
local t5  = 1/24 * e_(mu,nu,rho,si) * g_(1,mu,nu,rho,si, p1,p2);

* Mott-Streuung:

local mott = 1/(2*m)^2 * (g_(1,pp)+gi_(1)*m) * g_(1,0)
                       * (g_(1,p) +gi_(1)*m) * g_(1,0);
print;
.sort

trace4,1;

***b m;
***print[];
***.end

print +s;
.sort

id p.pp  = m^2 + beta^2*E^2*2*[sin(the/2)]^2;
id p(0)  = E;
id pp(0) = E;

bracket E,m;
print;
.end