coupling_comment:='( (g3 " strong interaction ")  (g "

weak interaction ") )
 $
delta_qnumber_list:='nil
 $
pa_decay_list:='((pr3 (45 (pr1 pr1) 2)) (pr4 (21 (pr1 pr1) 2
)) (pr5 (43 (pr1 pr1) 2)) (pr6 (41 (pr1 pr1) 2)) (pr7 (39 (
pr1 pr1) 2)) (pr8 (37 (pr1 pr1) 2)) (pr9 (20 (pr1 pr1) 2)) (
pr10 (35 (pr1 pr1) 2)) (pr11 (33 (pr1 pr1) 2)) (pr12 (31 (
pr1 pr1) 2) (30 (pr2 pr2) 2)) (pr15 (49 (pr3 pr1) 2) (48 (
pr3 pr1) 1) (47 (pr5 pr1) 2) (46 (pr5 pr1) 1) (25 (pr2 pr1)
2) (24 (pr2 pr1) 1)) (pr13 (28 (pr1 pr1) 2) (27 (pr2 pr2) 2)
) (pr14 (19 (pr1 pr1) 2) (18 (pr2 pr2) 2)) (pr16 (57 (pr3 
pr3) 2) (56 (pr5 pr3) 2) (55 (pr5 pr5) 2) (54 (pr6 pr3) 2) (
53 (pr6 pr5) 2) (52 (pr7 pr3) 2) (51 (pr7 pr5) 2) (50 (pr8 
pr3) 2) (44 (pr3 pr2) 2) (42 (pr5 pr2) 2) (40 (pr6 pr2) 2) (
38 (pr7 pr2) 2) (36 (pr8 pr2) 2) (34 (pr10 pr2) 2) (32 (pr11
pr2) 2) (29 (pr12 pr2) 2) (26 (pr13 pr2) 2) (23 (pr15 pr1) 2
) (22 (pr15 pr1) 2) (17 (pr4 pr2) 2) (16 (pr9 pr2) 2) (15 (
pr14 pr2) 2)))
 $
two_pa_vertex_list:='nil
 $
stable_pal0:='(p nut numu nue ef qu qd qs mu pr1 pr2 qc pr3
pr4 pr5 pr6 pr7 pr8 pr9 tau pr10 pr11 pr12 pr15 pr13 pr14 
pr16 qb w qt)
 $
mass_order_list:='(h0 qt z w qb pr16 pr14 pr13 pr15 pr12 
pr11 pr10 tau pr9 pr8 pr7 pr6 pr5 pr4 pr3 qc pr2 pr1 mu qs 
qd qu ef nue numu nut gsg p gs)
 $
pal_order_list:='(pr16 pr1 pr2 pr3 w pr15 pr5 h0 pr7 pr6 pa1
z pr14 pr13 pr12 pr8 pa5 pa4 pr9 pr4 pr11 pr10 p pa2 pa3)
 $
v_pal_list:='((14 (h0 -1 3)) (13 (h0 -1 4)) (12 (pa2 1 1) (
h0 -1 1)) (11 (w 1 1) (pa1 -1 1)) (10 (pa5 1 0) (pa4 0 1) (w
0 1)) (9 (pa4 1 0) (pa5 0 1) (w 1 0)) (8 (pa3 1 1) (pa1 -1 1
)) (7 (w 2 2)) (6 (w 1 1) (p -1 1) (z -1 1)) (5 (w 1 1) (pa1
-1 2)) (4 (w 1 1) (h0 -1 1)) (3 (z -1 2) (h0 -1 1)) (2 (w 1
1) (h0 -1 2)) (1 (z -1 2) (h0 -1 2)) (15 (pr16 -1 1) (pr14 
-1 1) (pr2 -1 1)) (16 (pr16 -1 1) (pr9 -1 1) (pr2 -1 1)) (17
(pr16 -1 1) (pr4 -1 1) (pr2 -1 1)) (18 (pr14 -1 1) (pr2 -1 2
)) (19 (pr14 -1 1) (pr1 1 1)) (20 (pr9 -1 1) (pr1 1 1)) (21
(pr4 -1 1) (pr1 1 1)) (22 (pr15 1 0) (pr16 -1 1) (pr1 0 1))
(23 (pr1 1 0) (pr16 -1 1) (pr15 0 1)) (24 (pr15 1 0) (pr2 -1
1) (pr1 0 1)) (25 (pr1 1 0) (pr2 -1 1) (pr15 0 1)) (26 (pr13
-1 1) (pr16 -1 1) (pr2 -1 1)) (27 (pr13 -1 1) (pr2 -1 2)) (
28 (pr13 -1 1) (pr1 1 1)) (29 (pr12 -1 1) (pr16 -1 1) (pr2 
-1 1)) (30 (pr12 -1 1) (pr2 -1 2)) (31 (pr12 -1 1) (pr1 1 1)
) (32 (pr11 -1 1) (pr16 -1 1) (pr2 -1 1)) (33 (pr11 -1 1) (
pr1 1 1)) (34 (pr10 -1 1) (pr16 -1 1) (pr2 -1 1)) (35 (pr10
-1 1) (pr1 1 1)) (36 (pr8 -1 1) (pr16 -1 1) (pr2 -1 1)) (37
(pr8 -1 1) (pr1 1 1)) (38 (pr7 -1 1) (pr16 -1 1) (pr2 -1 1))
(39 (pr7 -1 1) (pr1 1 1)) (40 (pr6 -1 1) (pr16 -1 1) (pr2 -1
1)) (41 (pr6 -1 1) (pr1 1 1)) (42 (pr5 -1 1) (pr16 -1 1) (
pr2 -1 1)) (43 (pr5 -1 1) (pr1 1 1)) (44 (pr3 -1 1) (pr16 -1
1) (pr2 -1 1)) (45 (pr3 -1 1) (pr1 1 1)) (46 (pr15 1 0) (pr5
-1 1) (pr1 0 1)) (47 (pr1 1 0) (pr5 -1 1) (pr15 0 1)) (48 (
pr15 1 0) (pr3 -1 1) (pr1 0 1)) (49 (pr1 1 0) (pr3 -1 1) (
pr15 0 1)) (50 (pr8 -1 1) (pr3 -1 1) (pr16 -1 1)) (51 (pr7 
-1 1) (pr5 -1 1) (pr16 -1 1)) (52 (pr7 -1 1) (pr3 -1 1) (
pr16 -1 1)) (53 (pr6 -1 1) (pr5 -1 1) (pr16 -1 1)) (54 (pr6
-1 1) (pr3 -1 1) (pr16 -1 1)) (55 (pr5 -1 2) (pr16 -1 1)) (
56 (pr5 -1 1) (pr3 -1 1) (pr16 -1 1)) (57 (pr3 -1 2) (pr16 
-1 1)))
 $
pa_rv_list:='((pr3 (g (1 44 45 48 49 50 52 54 56 57))) (pr16
(g (1 15 16 17 22 23 26 29 32 34 36 38 40 42 44 50 51 52 53
54 55 56 57))) (pr5 (g (1 42 43 46 47 51 53 55 56))) (pr6 (g
(1 40 41 53 54))) (pr7 (g (1 38 39 51 52))) (pr8 (g (1 36 37
50))) (pr1 (g (1 19 20 21 22 23 24 25 28 31 33 35 37 39 41 
43 45 46 47 48 49))) (pr15 (g (1 22 23 24 25 46 47 48 49)))
(pr2 (g (1 15 16 17 18 24 25 26 27 29 30 32 34 36 38 40 42 
44))) (pr10 (g (1 34 35))) (pr11 (g (1 32 33))) (pr12 (g (1
29 30 31))) (pr13 (g (1 26 27 28))) (pr4 (g (1 17 21))) (pr9
(g (1 16 20))) (pr14 (g (1 15 18 19))) (z (g (2 6 1) (1 3)))
(h0 (g (2 13 2 1) (1 14 12 4 3))) (w (g (2 7 6 5 2) (1 11 10
9 4))) (pa3 (g (1 8))) (pa1 (g (2 5) (1 11 8))) (pa4 (g (1 
10 9))) (pa5 (g (1 10 9))) (pa2 (g (1 12))) (p (g (2 6))))
 $
cvpd_list:='((14 0) (13 0) (11 1) (7 0) (6 0) (5 0) (4 0) (3
0) (2 0) (1 0) (15 0) (16 0) (17 0) (18 0) (19 0) (20 0) (21
0) (22 2) (23 2) (24 2) (25 2) (26 2) (27 2) (28 2) (29 2) (
30 2) (31 2) (32 2) (33 2) (34 2) (35 2) (36 2) (37 2) (38 2
) (39 2) (40 2) (41 2) (42 2) (43 2) (44 2) (45 2) (46 0) (
47 0) (48 0) (49 0) (50 0) (51 0) (52 0) (53 0) (54 0) (55 0
) (56 0) (57 0))
 $
gbsymmetrylist:='(strange baryon len lmuon ltau)
 $
fit_para_list:='((57 f49 f48) (56 f47 f46) (55 f45 f44) (54
f43 f42) (53 f41 f40) (52 f39 f38) (51 f37 f36) (50 f35 f34)
(49 f33 f32) (48 f33 f32) (47 f31 f30) (46 f31 f30) (45 f29)
(44 f28) (43 f27) (42 f26) (41 f25) (40 f24) (39 f23) (38 
f22) (37 f21) (36 f20) (35 f19) (34 f18) (33 f17) (32 f16) (
31 f15) (30 f14) (29 f13) (28 f12) (27 f11) (26 f10) (25 f9)
(24 f9) (23 f8) (22 f8) (21 f7) (20 f6) (19 f5) (18 f4) (17
f3) (16 f2) (15 f1))
 $
palist:='(pr1 pr2 pr3 pr4 pr5 pr6 pr7 pr8 pr9 pr10 pr11 pr12
pr13 pr14 pr15 pr16 gs w z p gsg h0 qb qs qd qt qc qu tau mu
ef nut numu nue)
 $
give_up_compact:='nil
 $
any_type_list:='((any_ext pr2 pr1) (any_nf pr1 pr2 pr3 pr4 
pr5 pr6 pr7 pr8 pr9 pr10 pr11 pr12 pr13 pr14 pr15 pr16 gs w
z p h0) (any pr1 pr2 pr3 pr4 pr5 pr6 pr7 pr8 pr9 pr10 pr11 
pr12 pr13 pr14 pr15 pr16 gs w z p h0 qb qs qd qt qc qu tau 
mu ef nut numu nue) (any_f qb qs qd qt qc qu tau mu ef nut 
numu nue) (any_l nue numu nut ef mu tau) (any_q qu qc qt qd
qs qb) (any_b h0 gs w z p) (any_b1 gs w z p) (any_b0 h0))
 $
effective_numbers:='nil
 $
bound_particle_list:='nil
 $
bound_state_vertex_list:='nil
 $
charge_op_sign_list:='((12 1) (12 1))
 $
color_vertices:='nil
 $
fcolor_vertices:='nil
 $
bound_pv_list:='nil
 $
fermionprog:='(sf)
 $
bosonprog:='(s11 s1 s0)
 $
propagatorl:='(s2 s11 s1 s0 sf)
 $
splist:='((s11 gs) (s1 w z p) (s0 pr1 pr2 pr4 pr9 pr14 pr16
gsg h0) (sf qt qc qu tau mu ef qb qs qd nut numu nue) (s2 
pr3 pr5 pr6 pr7 pr8 pr10 pr11 pr12 pr13 pr15))
 $
hvvvl:='((57 ((vvv x 57 r2 r1 vb v va vc) (b2 pr3 x 0 r1 v 
va) (b2 pr3 x 0 r2 vb vc) (b0 pr16 x 0 r3))) (56 ((vvv x 56
r2 r1 vb v va vc) (b2 pr5 x 0 r1 v va) (b2 pr3 x 0 r2 vb vc)
(b0 pr16 x 0 r3))) (55 ((vvv x 55 r2 r1 vb v va vc) (b2 pr5
x 0 r1 v va) (b2 pr5 x 0 r2 vb vc) (b0 pr16 x 0 r3))) (54 ((
vvv x 54 r2 r1 vb v va vc) (b2 pr6 x 0 r1 v va) (b2 pr3 x 0
r2 vb vc) (b0 pr16 x 0 r3))) (53 ((vvv x 53 r2 r1 vb v va vc
) (b2 pr6 x 0 r1 v va) (b2 pr5 x 0 r2 vb vc) (b0 pr16 x 0 r3
))) (52 ((vvv x 52 r2 r1 vb v va vc) (b2 pr7 x 0 r1 v va) (
b2 pr3 x 0 r2 vb vc) (b0 pr16 x 0 r3))) (51 ((vvv x 51 r2 r1
vb v va vc) (b2 pr7 x 0 r1 v va) (b2 pr5 x 0 r2 vb vc) (b0 
pr16 x 0 r3))) (50 ((vvv x 50 r2 r1 vb v va vc) (b2 pr8 x 0
r1 v va) (b2 pr3 x 0 r2 vb vc) (b0 pr16 x 0 r3))) (49 ((vvv
x 49 r2 r1 vb v va vc) (b0 pr1 x -1 r3) (b2 pr3 x 0 r2 vb vc
) (b2 pr15 x 1 r1 v va))) (48 ((vvv x 48 r2 r1 vb v va vc) (
b2 pr15 x -1 r1 v va) (b2 pr3 x 0 r2 vb vc) (b0 pr1 x 1 r3))
) (47 ((vvv x 47 r2 r1 vb v va vc) (b0 pr1 x -1 r3) (b2 pr5
x 0 r2 vb vc) (b2 pr15 x 1 r1 v va))) (46 ((vvv x 46 r2 r1 
vb v va vc) (b2 pr15 x -1 r1 v va) (b2 pr5 x 0 r2 vb vc) (b0
pr1 x 1 r3))) (45 ((vvv x 45 r2 va v) (b2 pr3 x 0 r1 v va) (
b0 pr1 x -1 r2) (b0 pr1 x 1 r3))) (44 ((vvv x 44 r2 va v) (
b2 pr3 x 0 r1 v va) (b0 pr16 x 0 r2) (b0 pr2 x 0 r3))) (43 (
(vvv x 43 r2 va v) (b2 pr5 x 0 r1 v va) (b0 pr1 x -1 r2) (b0
pr1 x 1 r3))) (42 ((vvv x 42 r2 va v) (b2 pr5 x 0 r1 v va) (
b0 pr16 x 0 r2) (b0 pr2 x 0 r3))) (41 ((vvv x 41 r2 va v) (
b2 pr6 x 0 r1 v va) (b0 pr1 x -1 r2) (b0 pr1 x 1 r3))) (40 (
(vvv x 40 r2 va v) (b2 pr6 x 0 r1 v va) (b0 pr16 x 0 r2) (b0
pr2 x 0 r3))) (39 ((vvv x 39 r2 va v) (b2 pr7 x 0 r1 v va) (
b0 pr1 x -1 r2) (b0 pr1 x 1 r3))) (38 ((vvv x 38 r2 va v) (
b2 pr7 x 0 r1 v va) (b0 pr16 x 0 r2) (b0 pr2 x 0 r3))) (37 (
(vvv x 37 r2 va v) (b2 pr8 x 0 r1 v va) (b0 pr1 x -1 r2) (b0
pr1 x 1 r3))) (36 ((vvv x 36 r2 va v) (b2 pr8 x 0 r1 v va) (
b0 pr16 x 0 r2) (b0 pr2 x 0 r3))) (35 ((vvv x 35 r2 va v) (
b2 pr10 x 0 r1 v va) (b0 pr1 x -1 r2) (b0 pr1 x 1 r3))) (34
((vvv x 34 r2 va v) (b2 pr10 x 0 r1 v va) (b0 pr16 x 0 r2) (
b0 pr2 x 0 r3))) (33 ((vvv x 33 r2 va v) (b2 pr11 x 0 r1 v 
va) (b0 pr1 x -1 r2) (b0 pr1 x 1 r3))) (32 ((vvv x 32 r2 va
v) (b2 pr11 x 0 r1 v va) (b0 pr16 x 0 r2) (b0 pr2 x 0 r3)))
(31 ((vvv x 31 r2 va v) (b2 pr12 x 0 r1 v va) (b0 pr1 x -1 
r2) (b0 pr1 x 1 r3))) (30 ((vvv x 30 r2 va v) (b2 pr12 x 0 
r1 v va) (b0 pr2 x 0 r2) (b0 pr2 x 0 r3))) (29 ((vvv x 29 r2
va v) (b2 pr12 x 0 r1 v va) (b0 pr16 x 0 r2) (b0 pr2 x 0 r3)
)) (28 ((vvv x 28 r2 va v) (b2 pr13 x 0 r1 v va) (b0 pr1 x 
-1 r2) (b0 pr1 x 1 r3))) (27 ((vvv x 27 r2 va v) (b2 pr13 x
0 r1 v va) (b0 pr2 x 0 r2) (b0 pr2 x 0 r3))) (26 ((vvv x 26
r2 va v) (b2 pr13 x 0 r1 v va) (b0 pr16 x 0 r2) (b0 pr2 x 0
r3))) (25 ((vvv x 25 r2 va v) (b0 pr1 x -1 r3) (b0 pr2 x 0 
r2) (b2 pr15 x 1 r1 v va))) (24 ((vvv x 24 r2 va v) (b2 pr15
x -1 r1 v va) (b0 pr2 x 0 r2) (b0 pr1 x 1 r3))) (23 ((vvv x
23 r2 va v) (b0 pr1 x -1 r3) (b0 pr16 x 0 r2) (b2 pr15 x 1 
r1 v va))) (22 ((vvv x 22 r2 va v) (b2 pr15 x -1 r1 v va) (
b0 pr16 x 0 r2) (b0 pr1 x 1 r3))) (21 ((vvv x 21) (b0 pr4 x
0 r1) (b0 pr1 x -1 r2) (b0 pr1 x 1 r3))) (20 ((vvv x 20) (b0
pr9 x 0 r1) (b0 pr1 x -1 r2) (b0 pr1 x 1 r3))) (19 ((vvv x 
19) (b0 pr14 x 0 r1) (b0 pr1 x -1 r2) (b0 pr1 x 1 r3))) (18
((vvv x 18) (b0 pr14 x 0 r1) (b0 pr2 x 0 r2) (b0 pr2 x 0 r3)
)) (17 ((vvv x 17) (b0 pr16 x 0 r1) (b0 pr4 x 0 r2) (b0 pr2
x 0 r3))) (16 ((vvv x 16) (b0 pr16 x 0 r1) (b0 pr9 x 0 r2) (
b0 pr2 x 0 r3))) (15 ((vvv x 15) (b0 pr16 x 0 r1) (b0 pr14 x
0 r2) (b0 pr2 x 0 r3))) (1 ((vvv x 1 v va) (b1 z x 0 r4 va)
(b1 z x 0 r3 v) (b0 h0 x 0 r1) (b0 h0 x 0 r2))) (2 ((vvv x 2
v va) (b1 w x -1 r4 va) (b1 w x 1 r3 v) (b0 h0 x 0 r2) (b0 
h0 x 0 r1))) (3 ((vvv x 3 v va) (b1 z x 0 r3 va) (b1 z x 0 
r2 v) (b0 h0 x 0 r1))) (4 ((vvv x 4 v va) (b1 w x -1 r3 va)
(b1 w x 1 r2 v) (b0 h0 x 0 r1))) (5 ((vvv x 5 pa1 v va vb vc
) (b1 w x 1 r1 v) (b1 w x -1 r2 va) (b1 pa1 x 0 r4 vc) (b1 
pa1 x 0 r3 vb))) (6 ((vvv x 6 v va vb vc) (b1 w x 1 r1 v) (
b1 w x -1 r2 va) (b1 p x 0 r4 vc) (b1 z x 0 r3 vb))) (7 ((
vvv x 7 v va vb vc) (b1 w x -1 r4 vc) (b1 w x -1 r3 vb) (b1
w x 1 r2 va) (b1 w x 1 r1 v))) (8 ((vvv x 8 l pa1 pa3 v) (
psi pa3 x -1 r1 id1) (psi pa3 x 1 r2 id1) (b1 pa1 x 0 r3 v))
) (9 ((vvv x 9 l pa5 pa4 v) (psi pa4 x -1 r1 id1) (psi pa5 x
1 r2 id1) (b1 w x -1 r3 v))) (10 ((vvv x 10 l pa4 pa5 v) (
psi pa5 x -1 r1 id1) (psi pa4 x 1 r2 id1) (b1 w x 1 r3 v)))
(11 ((vvv x 11 pa1 v va vb r1 r2 r3) (b1 w x 1 r1 va) (b1 w
x -1 r2 vb) (b1 pa1 x 0 r3 v))) (12 ((vvv x 12 l pa2) (psi 
pa2 x -1 r1 id1) (psi pa2 x 1 r2 id1) (b0 h0 x 0 r3))) (13 (
(vvv x 13) (b0 h0 x 0 r4) (b0 h0 x 0 r3) (b0 h0 x 0 r2) (b0
h0 x 0 r1))) (14 ((vvv x 14) (b0 h0 x 0 r3) (b0 h0 x 0 r2) (
b0 h0 x 0 r1))))
 $
htionlist:='((g (1 57 56 55 54 53 52 51 50 49 48 47 46 45 44
43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24
23 22 21 20 19 18 17 16 15 3 4 8 9 10 11 12 14) (2 1 2 5 6 7
13)))
 $
couplist:='(g3 g)
 $
maslist2:='((gsg mgsg 0) (z zm (quotient wm (cos theta))) (
gs mgs 0))
 $
paramlist:='(((c 0 1) (quotient (expt hm 2) 2)) ((c 0 5) (
quotient (minus (times (expt g 2) (expt hm 2))) (times 8 (
expt wm 2)))) (wangm3 (quotient (times (sqrt 2) fmb g) (
times 2 wm))) (wangm2 (quotient (times (sqrt 2) fms g) (
times 2 wm))) (wangm1 (quotient (times (sqrt 2) fmd g) (
times 2 wm))) (wangm6 (quotient (times (sqrt 2) fmtau g) (
times 2 wm))) (wangm5 (quotient (times (sqrt 2) fmmu g) (
times 2 wm))) (wangm4 (quotient (times (sqrt 2) fme g) (
times 2 wm))) (wangm9 (quotient (times (sqrt 2) fmt g) (
times 2 wm))) (wangm8 (quotient (times (sqrt 2) fmc g) (
times 2 wm))) (wangm7 (quotient (times (sqrt 2) fmu g) (
times 2 wm))) (v0 (quotient (times 2 wm) g)))
 $
phyinput:='(p1m w1d p2m w2d p3m w3d p4m w4d p5m w5d p6m w6d
p7m w7d p8m w8d p9m w9d p10m w10d p11m w11d p12m w12d p13m 
w13d p14m w14d p15m w15d p16m w16d g3 g theta wm hm fme fmmu
fmtau fmd fmu fms fmc fmb fmt)
 $
colorlist:='(pa6 pa4 pa5 pa2 pa3 qd qs qb qu qc qt gs)
 $
hnsuml:='(pr16 pr15 pr14 pr13 pr12 pr11 pr10 pr9 pr8 pr7 pr6
pr5 pr4 pr3 pr2 pr1 h0 gsg w gs)
 $
hsuml:='(pa3 pa2 pa5 pa4 pa6 pa1)
 $
hparameter:='((spacetime x) (sum pa3 pa2 pa5 pa4 pa6 pa1) (
nsum pr16 pr15 pr14 pr13 pr12 pr11 pr10 pr9 pr8 pr7 pr6 pr5
pr4 pr3 pr2 pr1 h0 gsg w gs) (single p z qt qc qu qb qs qd 
nut numu nue tau mu ef) (index3 ic1 ic2 ic3 ic4 ic5) (index2
id1 id2) (index1 v va vb vc vd vf va1 va2))
 $
pnametype:='(sum nsum single)
 $
hway:='(number dotnumber wayorder nameorder dotnameorder)
 $
hparaway:='((index3 wayorder) (index2 wayorder) (index1 
wayorder) (spacetime number) (sum dotnumber) (nsum 
dotwayorder) (single mdotnumber))
 $
hpsumrange:='((pa3 ef mu tau nue numu nut qd qs qb qu qc qt)
(pa2 qt qc qu qb qs qd tau mu ef) (pa5 nut numu nue qt qc qu
) (pa4 tau mu ef qb qs qd) (pa6 qu qc qt) (pa1 z p))
 $
hpindex:='(spacetime indices sum)
 $
singlist:='(pr1 pr2 pr3 pr4 pr5 pr6 pr7 pr8 pr9 pr10 pr11 
pr12 pr13 pr14 pr15 pr16 gs w z p gsg h0 qb qs qd qt qc qu 
tau mu ef nut numu nue)
 $
rangetable:='((pa5 pa6 1) (pa2 pa6 1) (pa2 pa4 1) (pa5 pa2 
pa6) (pa3 pa6 1) (pa3 pa4 1) (pa3 pa5 1) (pa3 pa2 1))
 $
hconstr:='((((0 pa3 pa1) (1 z nut numu nue) (2 pa3 z) (2 pa2
p) (1 pa1 pa2 qt qc qu qb qs qd tau mu ef pa4 pa6)) 8) (((0
pa5 pa4) (0 tau nut) (0 mu numu) (0 ef nue) (0 qb qt) (0 qs
qc) (0 qd qu)) 9 10))
 $
hconstr1:='(((pa3 pa1) 8) ((pa5 pa4) 9 10))
 $
patable1:='((b1 s11) (b1 s1) (b0 s0) (psi sf) (b2 s2))
 $
end;
 $
