add_fortran_list:='(af1);
 $
vpdlist:='((15 -1) (4 -1) (2 -1));
 $
funclist:='(v_ckm cg tp dmas);
 $
fermionprog:='(sf)
 $
bosonprog:='(s11 s1 s0)
 $
propagatorl:='(sf3 s11 s1 s0 sf)
 $
ss_dmas_list:='(s1)
 $
new_const_list:='nil
 $
bound_state_vertex_list:='nil
 $
color_vertices:='nil
 $
fcolor_vertices:='nil
 $
splist:='((s11 gs) (s1 jpsi2 jpsi w z p) (s0 pis gsg h0) (sf
qt qc qu tau mu ef qb qs qd nut numu nue))
 $
spflist:='(((sf3 l k m v va) i (quotient (plus (times 8 (
cons k v) (cons k va) (dkk k 1) m) (minus (times 9 (cons v 
va) m)) (times 3 (g l k v) (cons k va) (dkk k 1) m) (minus (
times 2 (g l k va) (cons k v) (dkk k 1) m)) (times 10 (g l k
) (cons k v) (cons k va) (dkk k 1)) (minus (times 9 (g l k)
(cons v va))) (minus (times 2 (g l v k va))) (minus (times 2
(g l v k) (cons k va) (dkk k 1) m)) (times 2 (g l v va) m) (
minus (times (g l v) (cons k va))) (times 3 (g l va k v)) (
times 3 (g l va k) (cons k v) (dkk k 1) m) (minus (times 3 (
g l va v) m)) (minus (times (g l va) (cons k v)))) 15)) ((
s11 r1 m v va) i (minus (cons v va))) ((s1 r1 m v va) i (
plus (times (dmas m 2) (cons r1 v) (cons r1 va)) (minus (
cons v va)))) ((s0 r1 m) i 1) ((sf l r1 m) i (plus (minus (g
l r1)) (minus m))))
 $
singlist:='(pis jpsi2 jpsi gs w z p gsg h0 qb qs qd qt qc qu
tau mu ef nut numu nue)
 $
wavefunction:='(u ub)
 $
hpsumrange:='((pa4 ef mu tau nue numu nut qd qs qb qu qc qt)
(pa3 qt qc qu qb qs qd tau mu ef) (pa6 nut numu nue qt qc qu
) (pa5 tau mu ef qb qs qd) (pa1 qd qs qb qu qc qt) (pa7 qu 
qc qt) (pa8 qd qs qb) (pa2 z p))
 $
feynmanrulelist:='(((20 va v) 0 (cjpsi2 (cons v va)) nil nil
) ((19 r4 r3 r1 va v) 0 (1 (plus (times (expt (cons r1 r3) 2
) (cons v va) cf3 dhjpmp4) (minus (times 2 (cons r1 r3) (
cons r1 r4) (cons v va) cf3 dhjpmp4)) (times (cons r1 r3) (
cons v va) cf1 dhjpmp2) (times (cons r1 r3) (eps r1 r3 v va)
cf11 dhjpmp4) (times (cons r1 r3) (eps r1 r3 v va) cf9 
dhjpmp4) (minus (times (cons r1 r3) (eps r1 r4 v va) cf11 
dhjpmp4)) (times (cons r1 r3) (eps r1 r4 v va) cf9 dhjpmp4)
(times (expt (cons r1 r4) 2) (cons v va) cf3 dhjpmp4) (minus
(times (cons r1 r4) (cons v va) cf1 dhjpmp2)) (minus (times
(cons r1 r4) (eps r1 r3 v va) cf11 dhjpmp4)) (minus (times (
cons r1 r4) (eps r1 r3 v va) cf9 dhjpmp4)) (times (cons r1 
r4) (eps r1 r4 v va) cf11 dhjpmp4) (minus (times (cons r1 r4
) (eps r1 r4 v va) cf9 dhjpmp4)) (times (cons r3 r4) (cons v
va) cf2 dhjpmp2) (times (cons r3 v) (cons r3 va) cf4 dhjpmp2
) (times (cons r3 v) (cons r3 va) cf5 dhjpmp2) (times (cons
r3 v) (cons r4 va) cf6 dhjpmp2) (times (cons r3 v) (cons r4
va) cf7 dhjpmp2) (times (cons r3 va) (cons r4 v) cf6 dhjpmp2
) (minus (times (cons r3 va) (cons r4 v) cf7 dhjpmp2)) (
times (cons r4 v) (cons r4 va) cf4 dhjpmp2) (minus (times (
cons r4 v) (cons r4 va) cf5 dhjpmp2)) (times (cons v va) cf0
) (times (eps r1 r3 v va) cf10 dhjpmp2) (times (eps r1 r3 v
va) cf8 dhjpmp2) (minus (times (eps r1 r4 v va) cf10 dhjpmp2
)) (times (eps r1 r4 v va) cf8 dhjpmp2) (times (eps r3 r4 v
va) cf12 dhjpmp2 i))) nil nil) ((1 v va vb vc) 0 (i (plus (
times 2 (cons v va) (cons vb vc) cas1) (times (cons v va) (
cons vb vc) cas2) (minus (times (cons v vb) (cons va vc) 
cas1)) (minus (times 2 (cons v vb) (cons va vc) cas2)) (
minus (times (cons v vc) (cons va vb) cas1)) (times (cons v
vc) (cons va vb) cas2))) nil ((1))) ((2 v va vb r1 r2 r3) 0
(1 (plus (times (cons r1 v) (cons va vb)) (minus (times (
cons r1 vb) (cons v va))) (minus (times (cons r2 v) (cons va
vb))) (times (cons r2 va) (cons v vb)) (minus (times (cons 
r3 va) (cons v vb))) (times (cons r3 vb) (cons v va)))) nil
((25))) ((3 l pa1 v) 1 ((quotient i 2) (times (g l v) (tp 1
3 3 ic1 id1 id2))) ((qt 51) (qc 50) (qu 49) (qb 39) (qs 38)
(qd 37)) ((37) (38) (39) (49) (50) (51))) ((4 v r3) 0 (-1 (
cons r3 v)) ((gs gsg 52)) ((52))) ((5 v va) 0 ((quotient i (
times 2 (expt (cos theta) 2))) (cons v va)) ((h0 z 3)) ((3))
) ((6 v va) 0 ((quotient i 2) (cons v va)) ((h0 w 2)) ((2)))
((7 v va) 0 ((quotient (times i wm) (expt (cos theta) 2)) (
cons v va)) ((h0 z 5)) ((5))) ((8 v va) 0 ((times i wm) (
cons v va)) ((h0 w 4)) ((4))) ((9 pa2 v va vb vc) 0 (cc9 (
plus (minus (times 2 (cons v va) (cons vb vc))) (times (cons
v vb) (cons va vc)) (times (cons v vc) (cons va vb)))) ((p 7
) (z 9)) ((7 (cc9 times (expt (sin theta) 2) i)) (9 (cc9 
times (expt (cos theta) 2) i)))) ((10 v va vb vc) 0 ((times
(cos theta) (sin theta) i) (plus (minus (times 2 (cons v va)
(cons vb vc))) (times (cons v vb) (cons va vc)) (times (cons
v vc) (cons va vb)))) nil ((8))) ((11 v va vb vc) 0 (i (plus
(times 2 (cons v va) (cons vb vc)) (minus (times (cons v vb)
(cons va vc))) (minus (times (cons v vc) (cons va vb))))) 
nil ((6))) ((12 l pa2 pa4 v) 4 (2 (cc12c0 (plus (times (g l
a v) cc12c2) (times (g l v) cc12c1)) (plus (times (g l a v)
cc12c1) (minus (times (g l a v) cc12c2)) (times (g l v) 
cc12c1) (times (g l v) cc12c2))) ((z qt 48) (z qc 47) (z qu
46) (z qb 33) (z qs 32) (z qd 31) (z nut 24) (z numu 23) (z
nue 22) (z tau 15) (z mu 14) (z ef 13) (p ef 10) (p mu 11) (
p tau 12) (p qd 28) (p qs 29) (p qb 30) (p qu 43) (p qc 44)
(p qt 45)) ((10 (cc12c1 . 1) (cc12c2 . 0) (cc12c0 minus (
times (sin theta) i))) (11 (cc12c1 . 1) (cc12c2 . 0) (cc12c0
minus (times (sin theta) i))) (12 (cc12c1 . 1) (cc12c2 . 0)
(cc12c0 minus (times (sin theta) i))) (13 (cc12c1 plus (
times 4 (expt (sin theta) 2)) (minus 1)) (cc12c2 . -1) (
cc12c0 quotient i (times 4 (cos theta)))) (14 (cc12c1 plus (
times 4 (expt (sin theta) 2)) (minus 1)) (cc12c2 . -1) (
cc12c0 quotient i (times 4 (cos theta)))) (15 (cc12c1 plus (
times 4 (expt (sin theta) 2)) (minus 1)) (cc12c2 . -1) (
cc12c0 quotient i (times 4 (cos theta)))) (22 (cc12c1 . 1) (
cc12c2 . 1) (cc12c0 quotient i (times 4 (cos theta)))) (23 (
cc12c1 . 1) (cc12c2 . 1) (cc12c0 quotient i (times 4 (cos 
theta)))) (24 (cc12c1 . 1) (cc12c2 . 1) (cc12c0 quotient i (
times 4 (cos theta)))) (28 (cc12c1 . 1) (cc12c2 . 0) (cc12c0
quotient (minus (times (sin theta) i)) 3)) (29 (cc12c1 . 1)
(cc12c2 . 0) (cc12c0 quotient (minus (times (sin theta) i))
3)) (30 (cc12c1 . 1) (cc12c2 . 0) (cc12c0 quotient (minus (
times (sin theta) i)) 3)) (31 (cc12c1 plus (times 4 (expt (
sin theta) 2)) (minus 3)) (cc12c2 . -3) (cc12c0 quotient i (
times 12 (cos theta)))) (32 (cc12c1 plus (times 4 (expt (sin
theta) 2)) (minus 3)) (cc12c2 . -3) (cc12c0 quotient i (
times 12 (cos theta)))) (33 (cc12c1 plus (times 4 (expt (sin
theta) 2)) (minus 3)) (cc12c2 . -3) (cc12c0 quotient i (
times 12 (cos theta)))) (43 (cc12c1 . 1) (cc12c2 . 0) (
cc12c0 quotient (times 2 (sin theta) i) 3)) (44 (cc12c1 . 1)
(cc12c2 . 0) (cc12c0 quotient (times 2 (sin theta) i) 3)) (
45 (cc12c1 . 1) (cc12c2 . 0) (cc12c0 quotient (times 2 (sin
theta) i) 3)) (46 (cc12c1 plus (times 8 (expt (cos theta) 2)
) (minus 5)) (cc12c2 . 3) (cc12c0 quotient i (times 12 (cos
theta)))) (47 (cc12c1 plus (times 8 (expt (cos theta) 2)) (
minus 5)) (cc12c2 . 3) (cc12c0 quotient i (times 12 (cos 
theta)))) (48 (cc12c1 plus (times 8 (expt (cos theta) 2)) (
minus 5)) (cc12c2 . 3) (cc12c0 quotient i (times 12 (cos 
theta))))) ((10 (cc12c1 quotient 1 2) (cc12c2 quotient 1 2)
(cc12c0 minus (times (sin theta) i))) (11 (cc12c1 quotient 1
2) (cc12c2 quotient 1 2) (cc12c0 minus (times (sin theta) i)
)) (12 (cc12c1 quotient 1 2) (cc12c2 quotient 1 2) (cc12c0 
minus (times (sin theta) i))) (13 (cc12c1 plus (times 2 (
expt (sin theta) 2)) (minus 1)) (cc12c2 times 2 (expt (sin 
theta) 2)) (cc12c0 quotient i (times 4 (cos theta)))) (14 (
cc12c1 plus (times 2 (expt (sin theta) 2)) (minus 1)) (
cc12c2 times 2 (expt (sin theta) 2)) (cc12c0 quotient i (
times 4 (cos theta)))) (15 (cc12c1 plus (times 2 (expt (sin
theta) 2)) (minus 1)) (cc12c2 times 2 (expt (sin theta) 2))
(cc12c0 quotient i (times 4 (cos theta)))) (22 (cc12c1 . 1)
(cc12c2 . 0) (cc12c0 quotient i (times 4 (cos theta)))) (23
(cc12c1 . 1) (cc12c2 . 0) (cc12c0 quotient i (times 4 (cos 
theta)))) (24 (cc12c1 . 1) (cc12c2 . 0) (cc12c0 quotient i (
times 4 (cos theta)))) (28 (cc12c1 quotient 1 2) (cc12c2 
quotient 1 2) (cc12c0 quotient (minus (times (sin theta) i))
3)) (29 (cc12c1 quotient 1 2) (cc12c2 quotient 1 2) (cc12c0
quotient (minus (times (sin theta) i)) 3)) (30 (cc12c1 
quotient 1 2) (cc12c2 quotient 1 2) (cc12c0 quotient (minus
(times (sin theta) i)) 3)) (31 (cc12c1 plus (times 2 (expt (
sin theta) 2)) (minus 3)) (cc12c2 times 2 (expt (sin theta)
2)) (cc12c0 quotient i (times 12 (cos theta)))) (32 (cc12c1
plus (times 2 (expt (sin theta) 2)) (minus 3)) (cc12c2 times
2 (expt (sin theta) 2)) (cc12c0 quotient i (times 12 (cos 
theta)))) (33 (cc12c1 plus (times 2 (expt (sin theta) 2)) (
minus 3)) (cc12c2 times 2 (expt (sin theta) 2)) (cc12c0 
quotient i (times 12 (cos theta)))) (43 (cc12c1 quotient 1 2
) (cc12c2 quotient 1 2) (cc12c0 quotient (times 2 (sin theta
) i) 3)) (44 (cc12c1 quotient 1 2) (cc12c2 quotient 1 2) (
cc12c0 quotient (times 2 (sin theta) i) 3)) (45 (cc12c1 
quotient 1 2) (cc12c2 quotient 1 2) (cc12c0 quotient (times
2 (sin theta) i) 3)) (46 (cc12c1 plus (times 4 (expt (cos 
theta) 2)) (minus 1)) (cc12c2 plus (times 4 (expt (cos theta
) 2)) (minus 4)) (cc12c0 quotient i (times 12 (cos theta))))
(47 (cc12c1 plus (times 4 (expt (cos theta) 2)) (minus 1)) (
cc12c2 plus (times 4 (expt (cos theta) 2)) (minus 4)) (
cc12c0 quotient i (times 12 (cos theta)))) (48 (cc12c1 plus
(times 4 (expt (cos theta) 2)) (minus 1)) (cc12c2 plus (
times 4 (expt (cos theta) 2)) (minus 4)) (cc12c0 quotient i
(times 12 (cos theta)))))) (((z pa8)) 1 ((quotient i (times
12 (cos theta))) (plus (minus (times 3 (g l a v))) (times 4
(g l v) (expt (sin theta) 2)) (minus (times 3 (g l v))))) ((
z qb 33) (z qb 33) (z qb 33))) (((z pa7)) 1 ((quotient i (
times 12 (cos theta))) (plus (times 3 (g l a v)) (times 8 (g
l v) (expt (cos theta) 2)) (minus (times 5 (g l v))))) ((z 
qt 48) (z qt 48) (z qt 48))) (((z pa1)) 2 ((quotient i (
times 12 (cos theta))) (plus (times (g l a v) cc12c2) (times
(g l v) cc12c1)) (plus (times (g l a v) cc12c1) (minus (
times (g l a v) cc12c2)) (times (g l v) cc12c1) (times (g l
v) cc12c2))) ((z qt 48) (z qt 48) (z qt 48) (z qt 48) (z qt
48) (z qt 48)) ((48 (cc12c1 plus (times 4 (expt (cos theta)
2)) (minus 1)) (cc12c2 plus (times 4 (expt (cos theta) 2)) (
minus 4))) (47 (cc12c1 plus (times 4 (expt (cos theta) 2)) (
minus 1)) (cc12c2 plus (times 4 (expt (cos theta) 2)) (minus
4))) (46 (cc12c1 plus (times 4 (expt (cos theta) 2)) (minus
1)) (cc12c2 plus (times 4 (expt (cos theta) 2)) (minus 4)))
(33 (cc12c1 plus (times 2 (expt (sin theta) 2)) (minus 3)) (
cc12c2 times 2 (expt (sin theta) 2))) (32 (cc12c1 plus (
times 2 (expt (sin theta) 2)) (minus 3)) (cc12c2 times 2 (
expt (sin theta) 2))) (31 (cc12c1 plus (times 2 (expt (sin 
theta) 2)) (minus 3)) (cc12c2 times 2 (expt (sin theta) 2)))
)) (((z pa5)) 2 (cc12c0 (plus (times (g l a v) cc12c2) (
times (g l v) cc12c1)) (plus (minus (times 2 (g l a v) (expt
(sin theta) 2))) (times (g l a v) cc12c1) (times 2 (g l v) (
expt (sin theta) 2)) (times (g l v) cc12c1))) ((z qb 33) (z
qb 33) (z qb 33) (z qb 33) (z qb 33) (z qb 33)) ((33 (cc12c1
plus (times 2 (expt (sin theta) 2)) (minus 3)) (cc12c0 
quotient i (times 12 (cos theta)))) (32 (cc12c1 plus (times
2 (expt (sin theta) 2)) (minus 3)) (cc12c0 quotient i (times
12 (cos theta)))) (31 (cc12c1 plus (times 2 (expt (sin theta
) 2)) (minus 3)) (cc12c0 quotient i (times 12 (cos theta))))
(15 (cc12c1 plus (times 2 (expt (sin theta) 2)) (minus 1)) (
cc12c0 quotient i (times 4 (cos theta)))) (14 (cc12c1 plus (
times 2 (expt (sin theta) 2)) (minus 1)) (cc12c0 quotient i
(times 4 (cos theta)))) (13 (cc12c1 plus (times 2 (expt (sin
theta) 2)) (minus 1)) (cc12c0 quotient i (times 4 (cos theta
)))))) (((p pa3) (p pa5) (p pa1) (p pa4)) 1 (cc12c0 (g l v))
((p qt 45) (p qt 45) (p qt 45) (p qt 45) (p qt 45) (p qt 45)
(p qt 45) (p qt 45) (p qt 45))) (((p pa8)) 1 ((quotient (
minus (times (sin theta) i)) 3) (g l v)) ((p qb 30) (p qb 30
) (p qb 30))) (((p pa6) (p pa7)) 1 ((quotient (times 2 (sin
theta) i) 3) (g l v)) ((p qt 45) (p qt 45) (p qt 45)))) ((13
l pa6 pa5 v) 1 ((quotient (times (sqrt 2) i) 4) (plus (g l a
v) (g l v))) ((nut tau 16) (numu mu 17) (nue ef 18) (qt qb 
34) (qc qs 35) (qu qd 36)) ((16) (17) (18) (34) (35) (36)))
((14 l pa5 pa6 v) 1 ((quotient (times (sqrt 2) i) 4) (plus (
g l a v) (g l v))) ((qd qu 42) (qs qc 41) (qb qt 40) (ef nue
21) (mu numu 20) (tau nut 19)) ((19) (20) (21) (40) (41) (42
))) ((15 pa2 v va vb r1 r2 r3) 0 (cc15 (plus (minus (times (
cons r1 v) (cons va vb))) (times (cons r1 vb) (cons v va)) (
times (cons r2 v) (cons va vb)) (minus (times (cons r2 va) (
cons v vb))) (times (cons r3 va) (cons v vb)) (minus (times
(cons r3 vb) (cons v va))))) ((p 26) (z 27)) ((26 (cc15 
times (sin theta) i)) (27 (cc15 times (cos theta) i)))) ((16
l pa3) 1 (cc16c0 1) ((qt 63) (qc 62) (qu 61) (qb 60) (qs 59)
(qd 58) (tau 55) (mu 54) (ef 53)) ((53 (cc16c0 quotient (
times fme i) (times 2 wm))) (54 (cc16c0 quotient (times fmmu
i) (times 2 wm))) (55 (cc16c0 quotient (times fmtau i) (
times 2 wm))) (58 (cc16c0 quotient (times fmd i) (times 2 wm
))) (59 (cc16c0 quotient (times fms i) (times 2 wm))) (60 (
cc16c0 quotient (times fmb i) (times 2 wm))) (61 (cc16c0 
quotient (times fmu i) (times 2 wm))) (62 (cc16c0 quotient (
times fmc i) (times 2 wm))) (63 (cc16c0 quotient (times fmt
i) (times 2 wm))))) ((17) 0 ((quotient (minus (times 3 (expt
hm 2) i)) (times 4 (expt wm 2))) 1) nil ((56))) ((18) 0 ((
quotient (minus (times 3 (expt hm 2) i)) (times 2 wm)) 1) 
nil ((57))));
 $
feynmanrulelist_new:='nil;
 $
czslist:='nil;
 $
ctp_list:='(((1 ic5 ic1 ic2 ic3 ic4) (cas1 (times (cg ic1 
ic2 ic5 3) (cg ic3 ic4 ic5 3))) (cas2 (times (cg ic3 ic2 ic5
3) (cg ic4 ic1 ic5 3))) ((plus cas1 cas2) (times (cg ic2 ic4
ic5 3) (cg ic1 ic3 ic5 3)))));
 $
end;
 $