n_particles:='3$ 
n_total_spin:='6$ 
http:='nil$ 
sub_model:='ok$ 
direct_model:='ok$ 
just_v_type:='nil$ 
nnvertex:='173$ 
pdl_add:='((pr1 0 p1m w1d 3 1 0 0 0 0 1 0 0 0) (pr2 2 p2m 
w2d 0 0 0 0 0 0 1 0 0 0) (pr3 4 p3m w3d 0 0 0 0 0 0 1 0 0 0)
(pr4 4 p4m w4d 3 1 0 0 0 0 1 0 0 0) (pr5 2 p5m w5d 3 1 0 0 0
0 1 0 0 0) (pr6 4 p6m w6d 0 0 0 0 0 0 1 0 0 0) (pr7 2 p7m 
w7d 0 0 0 0 0 0 1 0 0 0) (pr8 4 p8m w8d 0 0 0 0 0 0 1 0 0 0)
(pr9 4 p9m w9d 0 0 0 0 0 0 1 0 0 0) (pr10 2 p10m w10d 3 1 0
0 0 0 1 0 0 0) (pr11 4 p11m w11d 0 0 0 0 0 0 1 0 0 0) (pr12
4 p12m w12d 6 1 0 0 0 0 1 0 0 0) (pr13 4 p13m w13d 0 0 0 0 0
0 1 0 0 0) (pr14 6 p14m w14d 0 0 0 0 0 0 1 0 0 0) (pr15 2 
p15m w15d 0 0 0 0 0 0 1 0 0 0) (pr16 4 p16m w16d 3 1 0 0 0 0
1 0 0 0) (pr17 4 p17m w17d 0 0 0 0 0 0 1 0 0 0) (pr18 2 p18m
w18d 0 0 0 0 0 0 1 0 0 0) (pr19 6 p19m w19d 0 0 0 0 0 0 1 0
0 0) (pr20 8 p20m w20d 3 1 0 0 0 0 1 0 0 0) (pr21 2 p21m 
w21d 0 0 0 0 0 0 1 0 0 0) (pr22 0 p22m w22d 0 0 0 0 0 0 1 0
0 0))$ 
spflist_add:='(list)$ 
v_t_list:='nil$ 
chargelist:='((pr1 1) (pr2 0) (pr3 0) (pr4 1) (pr5 1) (pr6 0
) (pr7 0) (pr8 0) (pr9 0) (pr10 1) (pr11 0) (pr12 2) (pr13 0
) (pr14 0) (pr15 0) (pr16 1) (pr17 0) (pr18 0) (pr19 0) (
pr20 1) (pr21 0) (pr22 0) (gsg 0) (gs 0) (p 0) (z 0) (w 1) (
w -1) (z 0) (p 0) (qd (quotient (minus 1) 3)) (qu (quotient
2 3)) (qu (quotient 2 3)) (qd (quotient (minus 1) 3)) (ef -1
) (nue 0) (ef -1) (h0 0) (numu 0) (mu -1) (nut 0) (tau -1) (
mu -1) (tau -1) (qc (quotient 2 3)) (qs (quotient (minus 1)
3)) (qt (quotient 2 3)) (qb (quotient (minus 1) 3)) (qc (
quotient 2 3)) (qt (quotient 2 3)) (qs (quotient (minus 1) 3
)) (qb (quotient (minus 1) 3)))$ 
phyinput_add:='(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 p17m w17d p18m w18d
p19m w19d p20m w20d p21m w21d p22m w22d)$ 
platexname_add:='( 
      (pr1  "\pi") 
      (pr2  "\rho[770]^0") 
      (pr3  "f_2[1270]") 
      (pr4  "a_2[1320]") 
      (pr5  "\pi_1[1400]") 
      (pr6  "f_2[1430]") 
      (pr7  "\rho[1450]^0") 
      (pr8  "f^{'}_2[1525]") 
      (pr9  "f_2[1565]") 
      (pr10  "\pi_1[1600]") 
      (pr11  "X[1600]^0") 
      (pr12  "X[1600]") 
      (pr13  "f_2[1640]") 
      (pr14  "\rho_3[1690]^0") 
      (pr15  "\rho[1700]^0") 
      (pr16  "a_2[1700]") 
      (pr17  "f_2[1810]") 
      (pr18  "\rho[1900]^0") 
      (pr19  "\rho_3[1990]^0") 
      (pr20  "a_4[2040] ") 
      (pr21  "\rho[2150]^0") 
      (pr22  "\eta_{c}[2980]") 
     )$ 
 
 
inputlist_add:='( 
      (p1m  0.13957) 
      (w1d  0.0001) 
      (p2m  0.769) 
      (w2d  0.1507) 
      (p3m  1.275) 
      (w3d  0.1855) 
      (p4m  1.3181) 
      (w4d  0.1041) 
      (p5m  1.376) 
      (w5d  0.3) 
      (p6m  1.43) 
      (w6d  0.081) 
      (p7m  1.465) 
      (w7d  0.31) 
      (p8m  1.525) 
      (w8d  0.076) 
      (p9m  1.542) 
      (w9d  1.131) 
      (p10m  1.6) 
      (w10d  0.312) 
      (p11m  1.6) 
      (w11d  0.4) 
      (p12m  1.6) 
      (w12d  0.4) 
      (p13m  1.638) 
      (w13d  0.099) 
      (p14m  1.691) 
      (w14d  0.16) 
      (p15m  1.7) 
      (w15d  0.24) 
      (p16m  1.726) 
      (w16d  0.256) 
      (p17m  1.815) 
      (w17d  0.197) 
      (p18m  1.91) 
      (w18d  0.145) 
      (p19m  1.99) 
      (w19d  0.18) 
      (p20m  2.04) 
      (w20d  0.387) 
      (p21m  2.149) 
      (w21d  0.363) 
      (p22m  2.9792) 
      (w22d  0.013) 
     )$ 
 
pa_decay_list:='((pr2 (23 (pr1 pr1) 2)) (pr3 (59 (pr1 pr1) 2
)) (pr4 (97 (pr2 pr1) 2) (96 (pr2 pr1) 1)) (pr5 (45 (pr2 pr1
) 2) (44 (pr2 pr1) 1)) (pr6 (56 (pr1 pr1) 2)) (pr7 (93 (pr4
pr1) 2) (92 (pr4 pr1) 2) (20 (pr1 pr1) 2)) (pr8 (130 (pr4 
pr1) 2) (129 (pr4 pr1) 2) (81 (pr5 pr1) 2) (80 (pr5 pr1) 2)
(55 (pr1 pr1) 2)) (pr9 (126 (pr4 pr1) 2) (125 (pr4 pr1) 2) (
105 (pr2 pr2) 2) (79 (pr5 pr1) 2) (78 (pr5 pr1) 2) (54 (pr1
pr1) 2)) (pr10 (99 (pr3 pr1) 2) (98 (pr3 pr1) 1) (83 (pr6 
pr1) 2) (82 (pr6 pr1) 1) (40 (pr2 pr1) 2) (39 (pr2 pr1) 1))
(pr11 (123 (pr4 pr1) 2) (122 (pr4 pr1) 2) (104 (pr2 pr2) 2)
(77 (pr5 pr1) 2) (76 (pr5 pr1) 2) (53 (pr1 pr1) 2)) (pr12 (
121 (pr4 pr1) 2) (120 (pr4 pr1) 1) (75 (pr5 pr1) 2) (74 (pr5
pr1) 1) (52 (pr1 pr1) 2) (51 (pr1 pr1) 1)) (pr13 (118 (pr4 
pr1) 2) (117 (pr4 pr1) 2) (103 (pr2 pr2) 2) (73 (pr5 pr1) 2)
(72 (pr5 pr1) 2) (50 (pr1 pr1) 2)) (pr14 (151 (pr4 pr1) 2) (
150 (pr4 pr1) 2) (142 (pr5 pr1) 2) (141 (pr5 pr1) 2) (101 (
pr1 pr1) 2)) (pr15 (89 (pr4 pr1) 2) (88 (pr4 pr1) 2) (35 (
pr5 pr1) 2) (34 (pr5 pr1) 2) (17 (pr1 pr1) 2)) (pr16 (116 (
pr3 pr1) 2) (115 (pr3 pr1) 1) (114 (pr6 pr1) 2) (113 (pr6 
pr1) 1) (112 (pr8 pr1) 2) (111 (pr8 pr1) 1) (110 (pr9 pr1) 2
) (109 (pr9 pr1) 1) (71 (pr2 pr1) 2) (70 (pr2 pr1) 1) (69 (
pr7 pr1) 2) (68 (pr7 pr1) 1)) (pr17 (108 (pr4 pr1) 2) (107 (
pr4 pr1) 2) (102 (pr2 pr2) 2) (63 (pr5 pr1) 2) (62 (pr5 pr1)
2) (61 (pr10 pr1) 2) (60 (pr10 pr1) 2) (47 (pr1 pr1) 2)) (
pr18 (87 (pr4 pr1) 2) (86 (pr4 pr1) 2) (67 (pr16 pr1) 2) (66
(pr16 pr1) 2) (32 (pr5 pr1) 2) (31 (pr5 pr1) 2) (30 (pr10 
pr1) 2) (29 (pr10 pr1) 2) (16 (pr1 pr1) 2)) (pr19 (149 (pr4
pr1) 2) (148 (pr4 pr1) 2) (147 (pr16 pr1) 2) (146 (pr16 pr1)
2) (139 (pr5 pr1) 2) (138 (pr5 pr1) 2) (137 (pr10 pr1) 2) (
136 (pr10 pr1) 2) (100 (pr1 pr1) 2)) (pr20 (173 (pr3 pr1) 2)
(172 (pr3 pr1) 1) (171 (pr6 pr1) 2) (170 (pr6 pr1) 1) (169 (
pr8 pr1) 2) (168 (pr8 pr1) 1) (167 (pr9 pr1) 2) (166 (pr9 
pr1) 1) (165 (pr11 pr1) 2) (164 (pr11 pr1) 1) (163 (pr12 pr1
) 2) (162 (pr12 pr1) 1) (161 (pr13 pr1) 2) (160 (pr13 pr1) 1
) (159 (pr17 pr1) 2) (158 (pr17 pr1) 1) (157 (pr2 pr1) 2) (
156 (pr2 pr1) 1) (155 (pr7 pr1) 2) (154 (pr7 pr1) 1) (153 (
pr15 pr1) 2) (152 (pr15 pr1) 1)) (pr21 (106 (pr3 pr2) 2) (85
(pr4 pr1) 2) (84 (pr4 pr1) 2) (65 (pr16 pr1) 2) (64 (pr16 
pr1) 2) (27 (pr5 pr1) 2) (26 (pr5 pr1) 2) (25 (pr10 pr1) 2)
(24 (pr10 pr1) 2) (15 (pr1 pr1) 2)) (pr22 (145 (pr20 pr1) 2)
(144 (pr20 pr1) 2) (143 (pr14 pr2) 2) (140 (pr19 pr2) 2) (
135 (pr3 pr3) 2) (134 (pr4 pr4) 2) (133 (pr6 pr3) 2) (132 (
pr6 pr6) 2) (131 (pr8 pr3) 2) (128 (pr8 pr6) 2) (127 (pr9 
pr3) 2) (124 (pr9 pr6) 2) (119 (pr13 pr3) 2) (95 (pr4 pr5) 2
) (94 (pr4 pr5) 2) (91 (pr4 pr10) 2) (90 (pr4 pr10) 2) (58 (
pr4 pr1) 2) (57 (pr4 pr1) 2) (49 (pr16 pr1) 2) (48 (pr16 pr1
) 2) (46 (pr2 pr2) 2) (43 (pr5 pr5) 2) (42 (pr7 pr2) 2) (41
(pr7 pr7) 2) (38 (pr10 pr5) 2) (37 (pr10 pr5) 2) (36 (pr15 
pr2) 2) (33 (pr18 pr2) 2) (28 (pr21 pr2) 2) (22 (pr5 pr1) 2)
(21 (pr5 pr1) 2) (19 (pr10 pr1) 2) (18 (pr10 pr1) 2)))$ 
htionvvlist_add:='((173 ((pr1 -1 r3) (pr3 0 r2 vd vf) (pr20
1 r1 v va vb vc)) (plus (minus (times (cons r1 vf) (cons r2
va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f148 g)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f147 g)))) (172 ((pr20 -1 r1 v va vb vc) (pr3 0 r2
vd vf) (pr1 1 r3)) (plus (minus (times (cons r1 vf) (cons r2
va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f148 g)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f147 g)))) (171 ((pr1 -1 r3) (pr6 0 r2 vd vf) (pr20
1 r1 v va vb vc)) (plus (minus (times (cons r1 vf) (cons r2
va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f146 g)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f145 g)))) (170 ((pr20 -1 r1 v va vb vc) (pr6 0 r2
vd vf) (pr1 1 r3)) (plus (minus (times (cons r1 vf) (cons r2
va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f146 g)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f145 g)))) (169 ((pr1 -1 r3) (pr8 0 r2 vd vf) (pr20
1 r1 v va vb vc)) (plus (minus (times (cons r1 vf) (cons r2
va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f144 g)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f143 g)))) (168 ((pr20 -1 r1 v va vb vc) (pr8 0 r2
vd vf) (pr1 1 r3)) (plus (minus (times (cons r1 vf) (cons r2
va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f144 g)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f143 g)))) (167 ((pr1 -1 r3) (pr9 0 r2 vd vf) (pr20
1 r1 v va vb vc)) (plus (minus (times (cons r1 vf) (cons r2
va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f142 g)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f141 g)))) (166 ((pr20 -1 r1 v va vb vc) (pr9 0 r2
vd vf) (pr1 1 r3)) (plus (minus (times (cons r1 vf) (cons r2
va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f142 g)) (
minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (eps r1
r2 v vd) f141 g)))) (165 ((pr1 -1 r3) (pr11 0 r2 vd vf) (
pr20 1 r1 v va vb vc)) (plus (minus (times (cons r1 vf) (
cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f140
g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf) (
eps r1 r2 v vd) f139 g)))) (164 ((pr20 -1 r1 v va vb vc) (
pr11 0 r2 vd vf) (pr1 1 r3)) (plus (minus (times (cons r1 vf
) (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) 
f140 g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf
) (eps r1 r2 v vd) f139 g)))) (163 ((pr1 1 r3) (pr12 -1 r2 
vd vf) (pr20 1 r1 v va vb vc)) (plus (minus (times (cons r1
vf) (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd)
f138 g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf
) (eps r1 r2 v vd) f137 g)))) (162 ((pr20 -1 r1 v va vb vc)
(pr12 1 r2 vd vf) (pr1 -1 r3)) (plus (minus (times (cons r1
vf) (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd)
f138 g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf
) (eps r1 r2 v vd) f137 g)))) (161 ((pr1 -1 r3) (pr13 0 r2 
vd vf) (pr20 1 r1 v va vb vc)) (plus (minus (times (cons r1
vf) (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd)
f136 g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf
) (eps r1 r2 v vd) f135 g)))) (160 ((pr20 -1 r1 v va vb vc)
(pr13 0 r2 vd vf) (pr1 1 r3)) (plus (minus (times (cons r1 
vf) (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd)
f136 g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf
) (eps r1 r2 v vd) f135 g)))) (159 ((pr1 -1 r3) (pr17 0 r2 
vd vf) (pr20 1 r1 v va vb vc)) (plus (minus (times (cons r1
vf) (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd)
f134 g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf
) (eps r1 r2 v vd) f133 g)))) (158 ((pr20 -1 r1 v va vb vc)
(pr17 0 r2 vd vf) (pr1 1 r3)) (plus (minus (times (cons r1 
vf) (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd)
f134 g)) (minus (times (cons r2 va) (cons r2 vb) (cons vc vf
) (eps r1 r2 v vd) f133 g)))) (157 ((pr1 -1 r3) (pr2 0 r2 vd
) (pr20 1 r1 v va vb vc)) (times (cons r2 va) (cons r2 vb) (
cons r2 vc) (eps r1 r2 v vd) f132 g)) (156 ((pr20 -1 r1 v va
vb vc) (pr2 0 r2 vd) (pr1 1 r3)) (minus (times (cons r2 va)
(cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f132 g))) (155 ((
pr1 -1 r3) (pr7 0 r2 vd) (pr20 1 r1 v va vb vc)) (times (
cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1 r2 v vd) f131
g)) (154 ((pr20 -1 r1 v va vb vc) (pr7 0 r2 vd) (pr1 1 r3))
(minus (times (cons r2 va) (cons r2 vb) (cons r2 vc) (eps r1
r2 v vd) f131 g))) (153 ((pr1 -1 r3) (pr15 0 r2 vd) (pr20 1
r1 v va vb vc)) (times (cons r2 va) (cons r2 vb) (cons r2 vc
) (eps r1 r2 v vd) f130 g)) (152 ((pr20 -1 r1 v va vb vc) (
pr15 0 r2 vd) (pr1 1 r3)) (minus (times (cons r2 va) (cons 
r2 vb) (cons r2 vc) (eps r1 r2 v vd) f130 g))) (151 ((pr1 -1
r3) (pr4 1 r2 vc vd) (pr14 0 r1 v va vb)) (plus (times (cons
r1 vd) (cons r2 va) (cons r2 vb) (eps r1 r2 v vc) f129 g) (
times (cons r2 va) (cons vb vd) (eps r1 r2 v vc) f128 g))) (
150 ((pr14 0 r1 v va vb) (pr4 -1 r2 vc vd) (pr1 1 r3)) (plus
(minus (times (cons r1 vd) (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f129 g)) (minus (times (cons r2 va) (cons vb vd) (
eps r1 r2 v vc) f128 g)))) (149 ((pr1 -1 r3) (pr4 1 r2 vc vd
) (pr19 0 r1 v va vb)) (plus (times (cons r1 vd) (cons r2 va
) (cons r2 vb) (eps r1 r2 v vc) f127 g) (times (cons r2 va)
(cons vb vd) (eps r1 r2 v vc) f126 g))) (148 ((pr19 0 r1 v 
va vb) (pr4 -1 r2 vc vd) (pr1 1 r3)) (plus (minus (times (
cons r1 vd) (cons r2 va) (cons r2 vb) (eps r1 r2 v vc) f127
g)) (minus (times (cons r2 va) (cons vb vd) (eps r1 r2 v vc)
f126 g)))) (147 ((pr1 -1 r3) (pr16 1 r2 vc vd) (pr19 0 r1 v
va vb)) (plus (times (cons r1 vd) (cons r2 va) (cons r2 vb)
(eps r1 r2 v vc) f125 g) (times (cons r2 va) (cons vb vd) (
eps r1 r2 v vc) f124 g))) (146 ((pr19 0 r1 v va vb) (pr16 -1
r2 vc vd) (pr1 1 r3)) (plus (minus (times (cons r1 vd) (cons
r2 va) (cons r2 vb) (eps r1 r2 v vc) f125 g)) (minus (times
(cons r2 va) (cons vb vd) (eps r1 r2 v vc) f124 g)))) (145 (
(pr1 -1 r3) (pr22 0 r2) (pr20 1 r1 v va vb vc)) (times (cons
r2 v) (cons r2 va) (cons r2 vb) (cons r2 vc) f123 g)) (144 (
(pr20 -1 r1 v va vb vc) (pr22 0 r2) (pr1 1 r3)) (times (cons
r2 v) (cons r2 va) (cons r2 vb) (cons r2 vc) f123 g)) (143 (
(pr14 0 r1 v va vb) (pr2 0 r2 vc) (pr22 0 r3)) (minus (times
(cons r2 va) (cons r2 vb) (eps r1 r2 v vc) f122 g))) (142 ((
pr1 -1 r3) (pr5 1 r2 vc) (pr14 0 r1 v va vb)) (times (cons 
r2 va) (cons r2 vb) (eps r1 r2 v vc) f121 g i)) (141 ((pr14
0 r1 v va vb) (pr5 -1 r2 vc) (pr1 1 r3)) (minus (times (cons
r2 va) (cons r2 vb) (eps r1 r2 v vc) f121 g i))) (140 ((pr19
0 r1 v va vb) (pr2 0 r2 vc) (pr22 0 r3)) (minus (times (cons
r2 va) (cons r2 vb) (eps r1 r2 v vc) f120 g))) (139 ((pr1 -1
r3) (pr5 1 r2 vc) (pr19 0 r1 v va vb)) (times (cons r2 va) (
cons r2 vb) (eps r1 r2 v vc) f119 g i)) (138 ((pr19 0 r1 v 
va vb) (pr5 -1 r2 vc) (pr1 1 r3)) (minus (times (cons r2 va)
(cons r2 vb) (eps r1 r2 v vc) f119 g i))) (137 ((pr1 -1 r3)
(pr10 1 r2 vc) (pr19 0 r1 v va vb)) (times (cons r2 va) (
cons r2 vb) (eps r1 r2 v vc) f118 g i)) (136 ((pr19 0 r1 v 
va vb) (pr10 -1 r2 vc) (pr1 1 r3)) (minus (times (cons r2 va
) (cons r2 vb) (eps r1 r2 v vc) f118 g i))) (135 ((pr3 0 r1
v va) (pr3 0 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f117 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f116 g)))) (134 ((pr4 -1
r1 v va) (pr4 1 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f115 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f114 g)))) (133 ((pr6 0
r1 v va) (pr3 0 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f113 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f112 g)))) (132 ((pr6 0
r1 v va) (pr6 0 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f111 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f110 g)))) (131 ((pr8 0
r1 v va) (pr3 0 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f109 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f108 g)))) (130 ((pr1 -1
r3) (pr4 1 r2 vb vc) (pr8 0 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f107 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f106 g)))) (129 ((pr8 0
r1 v va) (pr4 -1 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f107 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f106 g)))) (128 ((pr8 0
r1 v va) (pr6 0 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f105 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f104 g)))) (127 ((pr9 0
r1 v va) (pr3 0 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f103 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f102 g)))) (126 ((pr1 -1
r3) (pr4 1 r2 vb vc) (pr9 0 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f101 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f100 g)))) (125 ((pr9 0
r1 v va) (pr4 -1 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f101 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f100 g)))) (124 ((pr9 0
r1 v va) (pr6 0 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f99 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f98 g)))) (123 ((pr1 -1
r3) (pr4 1 r2 vb vc) (pr11 0 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f97 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f96 g)))) (122 ((pr11 0
r1 v va) (pr4 -1 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f97 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f96 g)))) (121 ((pr1 -1
r3) (pr4 -1 r2 vb vc) (pr12 1 r1 v va)) (plus (minus (times
(cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f95 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f94 g)))) (120 ((pr12 -1
r1 v va) (pr4 1 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f95 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f94 g)))) (119 ((pr13 0
r1 v va) (pr3 0 r2 vb vc) (pr22 0 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f93 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f92 g)))) (118 ((pr1 -1
r3) (pr4 1 r2 vb vc) (pr13 0 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f91 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f90 g)))) (117 ((pr13 0
r1 v va) (pr4 -1 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f91 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f90 g)))) (116 ((pr1 -1
r3) (pr3 0 r2 vb vc) (pr16 1 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f89 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f88 g)))) (115 ((pr16 -1
r1 v va) (pr3 0 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f89 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f88 g)))) (114 ((pr1 -1
r3) (pr6 0 r2 vb vc) (pr16 1 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f87 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f86 g)))) (113 ((pr16 -1
r1 v va) (pr6 0 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f87 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f86 g)))) (112 ((pr1 -1
r3) (pr8 0 r2 vb vc) (pr16 1 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f85 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f84 g)))) (111 ((pr16 -1
r1 v va) (pr8 0 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f85 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f84 g)))) (110 ((pr1 -1
r3) (pr9 0 r2 vb vc) (pr16 1 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f83 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f82 g)))) (109 ((pr16 -1
r1 v va) (pr9 0 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f83 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f82 g)))) (108 ((pr1 -1
r3) (pr4 1 r2 vb vc) (pr17 0 r1 v va)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f81 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f80 g)))) (107 ((pr17 0
r1 v va) (pr4 -1 r2 vb vc) (pr1 1 r3)) (plus (minus (times (
cons r1 vc) (cons r2 va) (eps r1 r2 v vb) f81 g)) (minus (
times (cons va vc) (eps r1 r2 v vb) f80 g)))) (106 ((pr3 0 
r1 v va) (pr21 0 r2 vb) (pr2 0 r3 vc)) (plus (times (cons r1
vb) (cons r1 vc) (cons r2 v) (cons r2 va) f79 g) (times (
cons r1 vb) (cons r2 v) (cons va vc) f78 g) (times (cons r1
vc) (cons r2 v) (cons va vb) f76 g) (times (cons r2 v) (cons
r2 va) (cons vb vc) f77 g) (times (cons v vb) (cons va vc) 
f75 g))) (105 ((pr9 0 r1 v va) (pr2 0 r2 vb) (pr2 0 r3 vc))
(plus (times (cons r1 vb) (cons r1 vc) (cons r2 v) (cons r2
va) f74 g) (times (cons r1 vb) (cons r2 v) (cons va vc) f73
g) (times (cons r1 vc) (cons r2 v) (cons va vb) f71 g) (
times (cons r2 v) (cons r2 va) (cons vb vc) f72 g) (times (
cons v vb) (cons va vc) f70 g))) (104 ((pr11 0 r1 v va) (pr2
0 r2 vb) (pr2 0 r3 vc)) (plus (times (cons r1 vb) (cons r1 
vc) (cons r2 v) (cons r2 va) f69 g) (times (cons r1 vb) (
cons r2 v) (cons va vc) f68 g) (times (cons r1 vc) (cons r2
v) (cons va vb) f66 g) (times (cons r2 v) (cons r2 va) (cons
vb vc) f67 g) (times (cons v vb) (cons va vc) f65 g))) (103
((pr13 0 r1 v va) (pr2 0 r2 vb) (pr2 0 r3 vc)) (plus (times
(cons r1 vb) (cons r1 vc) (cons r2 v) (cons r2 va) f64 g) (
times (cons r1 vb) (cons r2 v) (cons va vc) f63 g) (times (
cons r1 vc) (cons r2 v) (cons va vb) f61 g) (times (cons r2
v) (cons r2 va) (cons vb vc) f62 g) (times (cons v vb) (cons
va vc) f60 g))) (102 ((pr17 0 r1 v va) (pr2 0 r2 vb) (pr2 0
r3 vc)) (plus (times (cons r1 vb) (cons r1 vc) (cons r2 v) (
cons r2 va) f59 g) (times (cons r1 vb) (cons r2 v) (cons va
vc) f58 g) (times (cons r1 vc) (cons r2 v) (cons va vb) f56
g) (times (cons r2 v) (cons r2 va) (cons vb vc) f57 g) (
times (cons v vb) (cons va vc) f55 g))) (101 ((pr14 0 r1 v 
va vb) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 
va) (cons r2 vb) f54 g)) (100 ((pr19 0 r1 v va vb) (pr1 -1 
r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 va) (cons r2 vb)
f53 g)) (99 ((pr1 -1 r3) (pr10 1 r2 vb) (pr3 0 r1 v va)) (
minus (times (cons r2 va) (eps r1 r2 v vb) f52 g i))) (98 ((
pr3 0 r1 v va) (pr10 -1 r2 vb) (pr1 1 r3)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f52 g i))) (97 ((pr1 -1 r3) (
pr2 0 r2 vb) (pr4 1 r1 v va)) (times (cons r2 va) (eps r1 r2
v vb) f51 g)) (96 ((pr4 -1 r1 v va) (pr2 0 r2 vb) (pr1 1 r3)
) (minus (times (cons r2 va) (eps r1 r2 v vb) f51 g))) (95 (
(pr22 0 r3) (pr5 -1 r2 vb) (pr4 1 r1 v va)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f50 g i))) (94 ((pr4 -1 r1 v va
) (pr5 1 r2 vb) (pr22 0 r3)) (minus (times (cons r2 va) (eps
r1 r2 v vb) f50 g i))) (93 ((pr1 -1 r3) (pr7 0 r2 vb) (pr4 1
r1 v va)) (times (cons r2 va) (eps r1 r2 v vb) f49 g)) (92 (
(pr4 -1 r1 v va) (pr7 0 r2 vb) (pr1 1 r3)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f49 g))) (91 ((pr22 0 r3) (pr10
-1 r2 vb) (pr4 1 r1 v va)) (minus (times (cons r2 va) (eps 
r1 r2 v vb) f48 g i))) (90 ((pr4 -1 r1 v va) (pr10 1 r2 vb)
(pr22 0 r3)) (minus (times (cons r2 va) (eps r1 r2 v vb) f48
g i))) (89 ((pr1 -1 r3) (pr15 0 r2 vb) (pr4 1 r1 v va)) (
times (cons r2 va) (eps r1 r2 v vb) f47 g)) (88 ((pr4 -1 r1
v va) (pr15 0 r2 vb) (pr1 1 r3)) (minus (times (cons r2 va)
(eps r1 r2 v vb) f47 g))) (87 ((pr1 -1 r3) (pr18 0 r2 vb) (
pr4 1 r1 v va)) (times (cons r2 va) (eps r1 r2 v vb) f46 g))
(86 ((pr4 -1 r1 v va) (pr18 0 r2 vb) (pr1 1 r3)) (minus (
times (cons r2 va) (eps r1 r2 v vb) f46 g))) (85 ((pr1 -1 r3
) (pr21 0 r2 vb) (pr4 1 r1 v va)) (times (cons r2 va) (eps 
r1 r2 v vb) f45 g)) (84 ((pr4 -1 r1 v va) (pr21 0 r2 vb) (
pr1 1 r3)) (minus (times (cons r2 va) (eps r1 r2 v vb) f45 g
))) (83 ((pr1 -1 r3) (pr10 1 r2 vb) (pr6 0 r1 v va)) (minus
(times (cons r2 va) (eps r1 r2 v vb) f44 g i))) (82 ((pr6 0
r1 v va) (pr10 -1 r2 vb) (pr1 1 r3)) (minus (times (cons r2
va) (eps r1 r2 v vb) f44 g i))) (81 ((pr1 -1 r3) (pr5 1 r2 
vb) (pr8 0 r1 v va)) (minus (times (cons r2 va) (eps r1 r2 v
vb) f43 g i))) (80 ((pr8 0 r1 v va) (pr5 -1 r2 vb) (pr1 1 r3
)) (minus (times (cons r2 va) (eps r1 r2 v vb) f43 g i))) (
79 ((pr1 -1 r3) (pr5 1 r2 vb) (pr9 0 r1 v va)) (minus (times
(cons r2 va) (eps r1 r2 v vb) f42 g i))) (78 ((pr9 0 r1 v va
) (pr5 -1 r2 vb) (pr1 1 r3)) (minus (times (cons r2 va) (eps
r1 r2 v vb) f42 g i))) (77 ((pr1 -1 r3) (pr5 1 r2 vb) (pr11
0 r1 v va)) (minus (times (cons r2 va) (eps r1 r2 v vb) f41
g i))) (76 ((pr11 0 r1 v va) (pr5 -1 r2 vb) (pr1 1 r3)) (
minus (times (cons r2 va) (eps r1 r2 v vb) f41 g i))) (75 ((
pr1 -1 r3) (pr5 -1 r2 vb) (pr12 1 r1 v va)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f40 g i))) (74 ((pr12 -1 r1 v 
va) (pr5 1 r2 vb) (pr1 1 r3)) (minus (times (cons r2 va) (
eps r1 r2 v vb) f40 g i))) (73 ((pr1 -1 r3) (pr5 1 r2 vb) (
pr13 0 r1 v va)) (minus (times (cons r2 va) (eps r1 r2 v vb)
f39 g i))) (72 ((pr13 0 r1 v va) (pr5 -1 r2 vb) (pr1 1 r3))
(minus (times (cons r2 va) (eps r1 r2 v vb) f39 g i))) (71 (
(pr1 -1 r3) (pr2 0 r2 vb) (pr16 1 r1 v va)) (times (cons r2
va) (eps r1 r2 v vb) f38 g)) (70 ((pr16 -1 r1 v va) (pr2 0 
r2 vb) (pr1 1 r3)) (minus (times (cons r2 va) (eps r1 r2 v 
vb) f38 g))) (69 ((pr1 -1 r3) (pr7 0 r2 vb) (pr16 1 r1 v va)
) (times (cons r2 va) (eps r1 r2 v vb) f37 g)) (68 ((pr16 -1
r1 v va) (pr7 0 r2 vb) (pr1 1 r3)) (minus (times (cons r2 va
) (eps r1 r2 v vb) f37 g))) (67 ((pr1 -1 r3) (pr18 0 r2 vb)
(pr16 1 r1 v va)) (times (cons r2 va) (eps r1 r2 v vb) f36 g
)) (66 ((pr16 -1 r1 v va) (pr18 0 r2 vb) (pr1 1 r3)) (minus
(times (cons r2 va) (eps r1 r2 v vb) f36 g))) (65 ((pr1 -1 
r3) (pr21 0 r2 vb) (pr16 1 r1 v va)) (times (cons r2 va) (
eps r1 r2 v vb) f35 g)) (64 ((pr16 -1 r1 v va) (pr21 0 r2 vb
) (pr1 1 r3)) (minus (times (cons r2 va) (eps r1 r2 v vb) 
f35 g))) (63 ((pr1 -1 r3) (pr5 1 r2 vb) (pr17 0 r1 v va)) (
minus (times (cons r2 va) (eps r1 r2 v vb) f34 g i))) (62 ((
pr17 0 r1 v va) (pr5 -1 r2 vb) (pr1 1 r3)) (minus (times (
cons r2 va) (eps r1 r2 v vb) f34 g i))) (61 ((pr1 -1 r3) (
pr10 1 r2 vb) (pr17 0 r1 v va)) (minus (times (cons r2 va) (
eps r1 r2 v vb) f33 g i))) (60 ((pr17 0 r1 v va) (pr10 -1 r2
vb) (pr1 1 r3)) (minus (times (cons r2 va) (eps r1 r2 v vb)
f33 g i))) (59 ((pr3 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (
times (cons r2 v) (cons r2 va) f32 g)) (58 ((pr1 -1 r3) (
pr22 0 r2) (pr4 1 r1 v va)) (times (cons r2 v) (cons r2 va)
f31 g)) (57 ((pr4 -1 r1 v va) (pr22 0 r2) (pr1 1 r3)) (times
(cons r2 v) (cons r2 va) f31 g)) (56 ((pr6 0 r1 v va) (pr1 
-1 r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 va) f30 g)) (
55 ((pr8 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2
v) (cons r2 va) f29 g)) (54 ((pr9 0 r1 v va) (pr1 -1 r2) (
pr1 1 r3)) (times (cons r2 v) (cons r2 va) f28 g)) (53 ((
pr11 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) (
cons r2 va) f27 g)) (52 ((pr1 -1 r3) (pr1 -1 r2) (pr12 1 r1
v va)) (times (cons r2 v) (cons r2 va) f26 g)) (51 ((pr12 -1
r1 v va) (pr1 1 r2) (pr1 1 r3)) (times (cons r2 v) (cons r2
va) f26 g)) (50 ((pr13 0 r1 v va) (pr1 -1 r2) (pr1 1 r3)) (
times (cons r2 v) (cons r2 va) f25 g)) (49 ((pr1 -1 r3) (
pr22 0 r2) (pr16 1 r1 v va)) (times (cons r2 v) (cons r2 va)
f24 g)) (48 ((pr16 -1 r1 v va) (pr22 0 r2) (pr1 1 r3)) (
times (cons r2 v) (cons r2 va) f24 g)) (47 ((pr17 0 r1 v va)
(pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) (cons r2 va) f23
g)) (46 ((pr2 0 r1 v) (pr2 0 r2 va) (pr22 0 r3)) (minus (
times (eps r1 r2 v va) f22 g))) (45 ((pr1 -1 r3) (pr2 0 r2 
va) (pr5 1 r1 v)) (times (eps r1 r2 v va) f21 g i)) (44 ((
pr5 -1 r1 v) (pr2 0 r2 va) (pr1 1 r3)) (minus (times (eps r1
r2 v va) f21 g i))) (43 ((pr5 -1 r1 v) (pr5 1 r2 va) (pr22 0
r3)) (minus (times (eps r1 r2 v va) f20 g))) (42 ((pr7 0 r1
v) (pr2 0 r2 va) (pr22 0 r3)) (minus (times (eps r1 r2 v va)
f19 g))) (41 ((pr7 0 r1 v) (pr7 0 r2 va) (pr22 0 r3)) (minus
(times (eps r1 r2 v va) f18 g))) (40 ((pr1 -1 r3) (pr2 0 r2
va) (pr10 1 r1 v)) (times (eps r1 r2 v va) f17 g i)) (39 ((
pr10 -1 r1 v) (pr2 0 r2 va) (pr1 1 r3)) (minus (times (eps 
r1 r2 v va) f17 g i))) (38 ((pr22 0 r3) (pr5 -1 r2 va) (pr10
1 r1 v)) (minus (times (eps r1 r2 v va) f16 g))) (37 ((pr10
-1 r1 v) (pr5 1 r2 va) (pr22 0 r3)) (minus (times (eps r1 r2
v va) f16 g))) (36 ((pr15 0 r1 v) (pr2 0 r2 va) (pr22 0 r3))
(minus (times (eps r1 r2 v va) f15 g))) (35 ((pr1 -1 r3) (
pr5 1 r2 va) (pr15 0 r1 v)) (times (eps r1 r2 v va) f14 g i)
) (34 ((pr15 0 r1 v) (pr5 -1 r2 va) (pr1 1 r3)) (minus (
times (eps r1 r2 v va) f14 g i))) (33 ((pr18 0 r1 v) (pr2 0
r2 va) (pr22 0 r3)) (minus (times (eps r1 r2 v va) f13 g)))
(32 ((pr1 -1 r3) (pr5 1 r2 va) (pr18 0 r1 v)) (times (eps r1
r2 v va) f12 g i)) (31 ((pr18 0 r1 v) (pr5 -1 r2 va) (pr1 1
r3)) (minus (times (eps r1 r2 v va) f12 g i))) (30 ((pr1 -1
r3) (pr10 1 r2 va) (pr18 0 r1 v)) (times (eps r1 r2 v va) 
f11 g i)) (29 ((pr18 0 r1 v) (pr10 -1 r2 va) (pr1 1 r3)) (
minus (times (eps r1 r2 v va) f11 g i))) (28 ((pr21 0 r1 v)
(pr2 0 r2 va) (pr22 0 r3)) (minus (times (eps r1 r2 v va) 
f10 g))) (27 ((pr1 -1 r3) (pr5 1 r2 va) (pr21 0 r1 v)) (
times (eps r1 r2 v va) f9 g i)) (26 ((pr21 0 r1 v) (pr5 -1 
r2 va) (pr1 1 r3)) (minus (times (eps r1 r2 v va) f9 g i)))
(25 ((pr1 -1 r3) (pr10 1 r2 va) (pr21 0 r1 v)) (times (eps 
r1 r2 v va) f8 g i)) (24 ((pr21 0 r1 v) (pr10 -1 r2 va) (pr1
1 r3)) (minus (times (eps r1 r2 v va) f8 g i))) (23 ((pr2 0
r1 v) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) f7 g)) (22
((pr1 -1 r3) (pr22 0 r2) (pr5 1 r1 v)) (times (cons r2 v) f6
g i)) (21 ((pr5 -1 r1 v) (pr22 0 r2) (pr1 1 r3)) (times (
cons r2 v) f6 g i)) (20 ((pr7 0 r1 v) (pr1 -1 r2) (pr1 1 r3)
) (times (cons r2 v) f5 g)) (19 ((pr1 -1 r3) (pr22 0 r2) (
pr10 1 r1 v)) (times (cons r2 v) f4 g i)) (18 ((pr10 -1 r1 v
) (pr22 0 r2) (pr1 1 r3)) (times (cons r2 v) f4 g i)) (17 ((
pr15 0 r1 v) (pr1 -1 r2) (pr1 1 r3)) (times (cons r2 v) f3 g
)) (16 ((pr18 0 r1 v) (pr1 -1 r2) (pr1 1 r3)) (times (cons 
r2 v) f2 g)) (15 ((pr21 0 r1 v) (pr1 -1 r2) (pr1 1 r3)) (
times (cons r2 v) f1 g)))$ 
vertices_add_comment:='"
In this model, We have added some effective 
vertices, to do some calculation related to some resonance states
"$ 
end; 
