Loading image file: /usr/lib/reduce/pslbuild/red/reduce.img 
Reduce (Free PSL version, revision 4715),  9-Aug-2018 ...

1: 
(drop model_state midmodel_add massolve equationlist feyrshow midmodel 
particle_table symbole_list_gen model_gen model_gen1 particle_table_gen 
feynman_rules model_input matterintinput symbole_list add_vertices 
physical_parameters physical_parameters_gen)

(sout alldiagramlist outadir fout amp_information option)

(showmore selfenergy justnumber ndiagfile ncolor nout nsy nsym state set_ext_sym
mkpsymboleok)

phys_cut := {}

phyinput_add := {}

ok

hpux

no

model_generator1

model_generator2

model_generator3

*** local variable n_coupling in procedure diagram_generator not used 

diagram_generator

write_diagram_information

write_next_amplitude

write_next_kinematic

generation_kinematics

process_chose23

gen_phyinput

get_x

amplitude_evaluation

model_chose

model_chose1

write_build_model_first

write_build_model

write_start_process

write_start_process1

write_model_information

write_process_information

process_chose

gmprocess

gen_wdlist

add_d0

add_d

d_alg

drop_zero_pa

mk_fpa_diag

gen_fit_parameter_list

find_n_it

gen_resonace_list

get_fasl_type

*** nonlocal use of undeclared variable new_parameter_list in procedure gmhtml 

gmhtml

ghtml

start_fdc
fdc_home:="/home/fdcpwa/fdc-pwa3.0/"$


op_system:='linux$


$

end$

model_state := model_state

"debug/massolve"

"debug/equationlist"

"debug/debug_input"

feyrshow := feynman_rule_show

drop := vertices_drop_list

add_vertices

"system/midmodel"

"system/midmodel_add"

"system/physical_parameters"

"debug/debug_input"

"diagram/"

"diagram/sout"

amp_information

option
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%---FDC---%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %  Now, you have just load the FDC package for Feynman Diagram Calculation  %
 %  The First step, you should choose a Physical model by perform:           %
 %  model_chose(aa);                                                         %
 %  Where aa refer to the file in which the model definition information is  %
 %  stored. The example of it could be found in file:                        %
 %  /home/fdcpwa/fdc-pwa3.0/example/model.def
 %  If the model is already built,  you can go directly to chose             %
 %  Physical process by perform:                                             %
 %  process_chose(aa);                                                       %
 %  Where aa refer to the file in which the process definition information   %
 %  is stored. The example of it could be found in file:                     %
 %  /home/fdcpwa/fdc-pwa3.0/example/process.def
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2:  You choose the The Phenomenology Lagrangian Model for $J/\Psi$ Physics
 in the directory ../model/, which is inputed by 
Jian-Xiong Wang and Rui-Chang Niu on April 15, 2022.
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%---FDC---%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %  The model is already built, so that you can go directly to chose         %
 %  Physical process by perform:                                             %
 %  process_chose(aa);                                                       %
 %  Where aa refer to the file in which the process definition information   %
 %  is stored. The example of it could be found in file:                     %
 %  /home/fdcpwa/fdc-pwa3.0/example/process.def
 %  If you want to rebuild the model with some other choice, You must delete %
 %  all lines in the file ../model/model_state
 %  Then run REDUCE and FDC again.                                           %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3: 
4: 
5: 
6: 6: 6: n_particles:='3$

 
n_total_spin:='6$

 
http:='nil$

 
sub_model:='ok$

 
direct_model:='ok$

 
just_v_type:='nil$

 
nnvertex:='81$

 
pdl_add:='((pr1 0 p1m w1d 0 0 0 0 0 0 1 0 0 0) (pr2 2 0 0 0
0 0 0 0 0 1 0 0 0) (pr3 2 p3m w3d 0 0 0 0 0 0 1 0 0 0) (pr4
2 p4m w4d 0 0 0 0 0 0 1 0 0 0) (pr5 2 p5m w5d 0 0 0 0 0 0 1
0 0 0) (pr6 2 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 2 p8m w8d 0 0 0 0 0 0 1 0 0 0) (pr9
2 p9m w9d 0 0 0 0 0 0 1 0 0 0) (pr10 6 p10m w10d 0 0 0 0 0 0
1 0 0 0) (pr11 2 p11m w11d 0 0 0 0 0 0 1 0 0 0) (pr12 6 p12m
w12d 0 0 0 0 0 0 1 0 0 0) (pr13 0 p13m w13d 0 0 0 0 0 0 1 0
0 0))$

 
spflist_add:='(list)$

 
v_t_list:='nil$

 
chargelist:='((pr1 0) (pr2 0) (pr3 0) (pr4 0) (pr5 0) (pr6 0
) (pr7 0) (pr8 0) (pr9 0) (pr10 0) (pr11 0) (pr12 0) (pr13 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 p3m w3d p4m w4d p5m w5d p6m w6d p7m
w7d p8m w8d p9m w9d p10m w10d p11m w11d p12m w12d p13m w13d)
$

 
platexname_add:='( 
      (pr1  "\eta") 
      (pr2  "\gamma") 
      (pr3  "\omega[782]") 
      (pr4  "\phi[1020]") 
      (pr5  "h_1[1170]") 
      (pr6  "h_1[1380]") 
      (pr7  "\omega[1420]") 
      (pr8  "h_1[1595]") 
      (pr9  "\omega_[1650]") 
      (pr10  "\omega_3[1670]") 
      (pr11  "\phi[1680]") 
      (pr12  "\phi_3[1850]") 
      (pr13  "\eta_{c}[2980]") 
     )$

 
 
 
inputlist_add:='( 
      (p1m  0.5473) 
      (w1d  0.0001) 
      (p3m  0.78194) 
      (w3d  0.0084) 
      (p4m  1.01941) 
      (w4d  0.0043) 
      (p5m  1.17) 
      (w5d  0.36) 
      (p6m  1.386) 
      (w6d  0.17) 
      (p7m  1.419) 
      (w7d  0.174) 
      (p8m  1.594) 
      (w8d  0.384) 
      (p9m  1.652) 
      (w9d  0.05) 
      (p10m  1.667) 
      (w10d  0.168) 
      (p11m  1.68) 
      (w11d  0.15) 
      (p12m  1.854) 
      (w12d  0.087) 
      (p13m  2.9792) 
      (w13d  0.013) 
     )$

 
 
pa_decay_list:='((pr3 (67 (pr2 pr1) 2) (15 (ef ef) 2)) (pr4
(64 (pr2 pr1) 2) (16 (ef ef) 2)) (pr5 (60 (pr2 pr1) 2)) (pr6
(55 (pr2 pr1) 2) (53 (pr3 pr1) 2)) (pr7 (48 (pr2 pr1) 2) (46
(pr3 pr1) 2) (17 (ef ef) 2)) (pr8 (40 (pr2 pr1) 2) (38 (pr3
pr1) 2) (36 (pr4 pr1) 2)) (pr9 (33 (pr2 pr1) 2) (31 (pr3 pr1
) 2) (29 (pr4 pr1) 2) (18 (ef ef) 2)) (pr10 (81 (pr2 pr1) 2)
(79 (pr3 pr1) 2) (77 (pr4 pr1) 2)) (pr11 (26 (pr2 pr1) 2) (
24 (pr3 pr1) 2) (22 (pr4 pr1) 2) (19 (ef ef) 2)) (pr12 (74 (
pr2 pr1) 2) (72 (pr3 pr1) 2) (70 (pr4 pr1) 2) (68 (pr5 pr1)
2)) (pr13 (80 (pr10 pr2) 2) (78 (pr10 pr3) 2) (76 (pr10 pr4)
2) (75 (pr10 pr5) 2) (73 (pr12 pr2) 2) (71 (pr12 pr3) 2) (69
(pr12 pr4) 2) (66 (pr3 pr2) 2) (65 (pr3 pr3) 2) (63 (pr4 pr2
) 2) (62 (pr4 pr3) 2) (61 (pr4 pr4) 2) (59 (pr5 pr2) 2) (58
(pr5 pr3) 2) (57 (pr5 pr4) 2) (56 (pr5 pr5) 2) (54 (pr6 pr2)
2) (52 (pr6 pr3) 2) (51 (pr6 pr4) 2) (50 (pr6 pr5) 2) (49 (
pr6 pr6) 2) (47 (pr7 pr2) 2) (45 (pr7 pr3) 2) (44 (pr7 pr4)
2) (43 (pr7 pr5) 2) (42 (pr7 pr6) 2) (41 (pr7 pr7) 2) (39 (
pr8 pr2) 2) (37 (pr8 pr3) 2) (35 (pr8 pr4) 2) (34 (pr8 pr5)
2) (32 (pr9 pr2) 2) (30 (pr9 pr3) 2) (28 (pr9 pr4) 2) (27 (
pr9 pr5) 2) (25 (pr11 pr2) 2) (23 (pr11 pr3) 2) (21 (pr11 
pr4) 2) (20 (pr11 pr5) 2)))$

 
htionvvlist_add:='((81 ((pr10 0 r1 v va vb) (pr2 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f82 g))) (80 ((pr10 0 r1 v va vb) (pr2 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f81 g))) (79 ((pr10 0 r1 v va vb) (pr3 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f80 g))) (78 ((pr10 0 r1 v va vb) (pr3 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f79 g))) (77 ((pr10 0 r1 v va vb) (pr4 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f78 g))) (76 ((pr10 0 r1 v va vb) (pr4 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f77 g))) (75 ((pr10 0 r1 v va vb) (pr5 0 r2 vc) (
pr13 0 r3)) (plus (times (cons r1 vc) (cons r2 v) (cons r2 
va) (cons r2 vb) f76 g) (times (cons r2 v) (cons r2 va) (
cons vb vc) f75 g))) (74 ((pr12 0 r1 v va vb) (pr2 0 r2 vc)
(pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f74 g))) (73 ((pr12 0 r1 v va vb) (pr2 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f73 g))) (72 ((pr12 0 r1 v va vb) (pr3 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f72 g))) (71 ((pr12 0 r1 v va vb) (pr3 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f71 g))) (70 ((pr12 0 r1 v va vb) (pr4 0 r2 vc) (
pr1 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1 
r2 v vc) f70 g))) (69 ((pr12 0 r1 v va vb) (pr4 0 r2 vc) (
pr13 0 r3)) (minus (times (cons r2 va) (cons r2 vb) (eps r1
r2 v vc) f69 g))) (68 ((pr12 0 r1 v va vb) (pr5 0 r2 vc) (
pr1 0 r3)) (plus (times (cons r1 vc) (cons r2 v) (cons r2 va
) (cons r2 vb) f68 g) (times (cons r2 v) (cons r2 va) (cons
vb vc) f67 g))) (67 ((pr3 0 r1 v) (pr2 0 r2 va) (pr1 0 r3))
(minus (times (eps r1 r2 v va) f66 g))) (66 ((pr3 0 r1 v) (
pr2 0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) f65
g))) (65 ((pr3 0 r1 v) (pr3 0 r2 va) (pr13 0 r3)) (minus (
times (eps r1 r2 v va) f64 g))) (64 ((pr4 0 r1 v) (pr2 0 r2
va) (pr1 0 r3)) (minus (times (eps r1 r2 v va) f63 g))) (63
((pr4 0 r1 v) (pr2 0 r2 va) (pr13 0 r3)) (minus (times (eps
r1 r2 v va) f62 g))) (62 ((pr4 0 r1 v) (pr3 0 r2 va) (pr13 0
r3)) (minus (times (eps r1 r2 v va) f61 g))) (61 ((pr4 0 r1
v) (pr4 0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va)
f60 g))) (60 ((pr5 0 r1 v) (pr2 0 r2 va) (pr1 0 r3)) (
quotient (plus (times (cons r1 r1) (cons v va) f59 g) (times
2 (cons r1 va) (cons r2 v) f59 g) (minus (times (cons r3 r3)
(cons v va) f59 g))) 2)) (59 ((pr5 0 r1 v) (pr2 0 r2 va) (
pr13 0 r3)) (quotient (plus (times (cons r1 r1) (cons v va)
f58 g) (times 2 (cons r1 va) (cons r2 v) f58 g) (minus (
times (cons r3 r3) (cons v va) f58 g))) 2)) (58 ((pr5 0 r1 v
) (pr3 0 r2 va) (pr13 0 r3)) (plus (times (cons r1 va) (cons
r2 v) f57 g) (times (cons v va) f56 g))) (57 ((pr5 0 r1 v) (
pr4 0 r2 va) (pr13 0 r3)) (plus (times (cons r1 va) (cons r2
v) f55 g) (times (cons v va) f54 g))) (56 ((pr5 0 r1 v) (pr5
0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) f53 g))
) (55 ((pr6 0 r1 v) (pr2 0 r2 va) (pr1 0 r3)) (quotient (
plus (times (cons r1 r1) (cons v va) f52 g) (times 2 (cons 
r1 va) (cons r2 v) f52 g) (minus (times (cons r3 r3) (cons v
va) f52 g))) 2)) (54 ((pr6 0 r1 v) (pr2 0 r2 va) (pr13 0 r3)
) (quotient (plus (times (cons r1 r1) (cons v va) f51 g) (
times 2 (cons r1 va) (cons r2 v) f51 g) (minus (times (cons
r3 r3) (cons v va) f51 g))) 2)) (53 ((pr6 0 r1 v) (pr3 0 r2
va) (pr1 0 r3)) (plus (times (cons r1 va) (cons r2 v) f50 g)
(times (cons v va) f49 g))) (52 ((pr6 0 r1 v) (pr3 0 r2 va)
(pr13 0 r3)) (plus (times (cons r1 va) (cons r2 v) f48 g) (
times (cons v va) f47 g))) (51 ((pr6 0 r1 v) (pr4 0 r2 va) (
pr13 0 r3)) (plus (times (cons r1 va) (cons r2 v) f46 g) (
times (cons v va) f45 g))) (50 ((pr6 0 r1 v) (pr5 0 r2 va) (
pr13 0 r3)) (minus (times (eps r1 r2 v va) f44 g))) (49 ((
pr6 0 r1 v) (pr6 0 r2 va) (pr13 0 r3)) (minus (times (eps r1
r2 v va) f43 g))) (48 ((pr7 0 r1 v) (pr2 0 r2 va) (pr1 0 r3)
) (minus (times (eps r1 r2 v va) f42 g))) (47 ((pr7 0 r1 v)
(pr2 0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) 
f41 g))) (46 ((pr7 0 r1 v) (pr3 0 r2 va) (pr1 0 r3)) (minus
(times (eps r1 r2 v va) f40 g))) (45 ((pr7 0 r1 v) (pr3 0 r2
va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) f39 g))) (44
((pr7 0 r1 v) (pr4 0 r2 va) (pr13 0 r3)) (minus (times (eps
r1 r2 v va) f38 g))) (43 ((pr7 0 r1 v) (pr5 0 r2 va) (pr13 0
r3)) (plus (times (cons r1 va) (cons r2 v) f37 g) (times (
cons v va) f36 g))) (42 ((pr7 0 r1 v) (pr6 0 r2 va) (pr13 0
r3)) (plus (times (cons r1 va) (cons r2 v) f35 g) (times (
cons v va) f34 g))) (41 ((pr7 0 r1 v) (pr7 0 r2 va) (pr13 0
r3)) (minus (times (eps r1 r2 v va) f33 g))) (40 ((pr8 0 r1
v) (pr2 0 r2 va) (pr1 0 r3)) (quotient (plus (times (cons r1
r1) (cons v va) f32 g) (times 2 (cons r1 va) (cons r2 v) f32
g) (minus (times (cons r3 r3) (cons v va) f32 g))) 2)) (39 (
(pr8 0 r1 v) (pr2 0 r2 va) (pr13 0 r3)) (quotient (plus (
times (cons r1 r1) (cons v va) f31 g) (times 2 (cons r1 va)
(cons r2 v) f31 g) (minus (times (cons r3 r3) (cons v va) 
f31 g))) 2)) (38 ((pr8 0 r1 v) (pr3 0 r2 va) (pr1 0 r3)) (
plus (times (cons r1 va) (cons r2 v) f30 g) (times (cons v 
va) f29 g))) (37 ((pr8 0 r1 v) (pr3 0 r2 va) (pr13 0 r3)) (
plus (times (cons r1 va) (cons r2 v) f28 g) (times (cons v 
va) f27 g))) (36 ((pr8 0 r1 v) (pr4 0 r2 va) (pr1 0 r3)) (
plus (times (cons r1 va) (cons r2 v) f26 g) (times (cons v 
va) f25 g))) (35 ((pr8 0 r1 v) (pr4 0 r2 va) (pr13 0 r3)) (
plus (times (cons r1 va) (cons r2 v) f24 g) (times (cons v 
va) f23 g))) (34 ((pr8 0 r1 v) (pr5 0 r2 va) (pr13 0 r3)) (
minus (times (eps r1 r2 v va) f22 g))) (33 ((pr9 0 r1 v) (
pr2 0 r2 va) (pr1 0 r3)) (minus (times (eps r1 r2 v va) f21
g))) (32 ((pr9 0 r1 v) (pr2 0 r2 va) (pr13 0 r3)) (minus (
times (eps r1 r2 v va) f20 g))) (31 ((pr9 0 r1 v) (pr3 0 r2
va) (pr1 0 r3)) (minus (times (eps r1 r2 v va) f19 g))) (30
((pr9 0 r1 v) (pr3 0 r2 va) (pr13 0 r3)) (minus (times (eps
r1 r2 v va) f18 g))) (29 ((pr9 0 r1 v) (pr4 0 r2 va) (pr1 0
r3)) (minus (times (eps r1 r2 v va) f17 g))) (28 ((pr9 0 r1
v) (pr4 0 r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va)
f16 g))) (27 ((pr9 0 r1 v) (pr5 0 r2 va) (pr13 0 r3)) (plus
(times (cons r1 va) (cons r2 v) f15 g) (times (cons v va) 
f14 g))) (26 ((pr11 0 r1 v) (pr2 0 r2 va) (pr1 0 r3)) (minus
(times (eps r1 r2 v va) f13 g))) (25 ((pr11 0 r1 v) (pr2 0 
r2 va) (pr13 0 r3)) (minus (times (eps r1 r2 v va) f12 g)))
(24 ((pr11 0 r1 v) (pr3 0 r2 va) (pr1 0 r3)) (minus (times (
eps r1 r2 v va) f11 g))) (23 ((pr11 0 r1 v) (pr3 0 r2 va) (
pr13 0 r3)) (minus (times (eps r1 r2 v va) f10 g))) (22 ((
pr11 0 r1 v) (pr4 0 r2 va) (pr1 0 r3)) (minus (times (eps r1
r2 v va) f9 g))) (21 ((pr11 0 r1 v) (pr4 0 r2 va) (pr13 0 r3
)) (minus (times (eps r1 r2 v va) f8 g))) (20 ((pr11 0 r1 v)
(pr5 0 r2 va) (pr13 0 r3)) (plus (times (cons r1 va) (cons 
r2 v) f7 g) (times (cons v va) f6 g))) (19 ((pr11 0 r1 v) (
ef -1 r2) (ef 1 r3)) (times (g l v) f5 g)) (18 ((pr9 0 r1 v)
(ef -1 r2) (ef 1 r3)) (times (g l v) f4 g)) (17 ((pr7 0 r1 v
) (ef -1 r2) (ef 1 r3)) (times (g l v) f3 g)) (16 ((pr4 0 r1
v) (ef -1 r2) (ef 1 r3)) (times (g l v) f2 g)) (15 ((pr3 0 
r1 v) (ef -1 r2) (ef 1 r3)) (times (g l 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;

7: 
matterinteraction:=el(-1)*er(1)*hig(1)*c(1,4)
+hig(-1)*er(-1)*el(1)*c(-1,4)
+hig(-1)*hig(1)*c(0,1)
+hig(-1)*hig(1)*hig(-1)*hig(1)*c(0,5)
+hig(-1)*q1dr(-1)*q1l(1)*c(-1,3)
+higtid(-1)*q1ur(-1)*q1l(1)*c(-1,2)
+q1l(-1)*q1dr(1)*hig(1)*c(1,3)
+q1l(-1)*q1ur(1)*higtid(1)*c(1,2)$


*** el declared operator 

*** er declared operator 

*** hig declared operator 

*** q1dr declared operator 

*** q1l declared operator 

*** higtid declared operator 

*** q1ur declared operator 


 $


gl_s_list:='(    )$


 $



end$

8: ;

end;

9: platexname:='(
        (p "\gamma") 
        (pa1 "\gamma~Z^0") 
        (z "Z") 
        (w w) 
        (gs g) 
        (nue "{\nu}_e") 
        (numu "{\nu}_{\mu}") 
        (nut "{\nu}_{\tau}") 
        (ef e)
        (mu "\mu")
        (tau "\tau")
        (qu u) 
        (qc c) 
        (qt t) 
        (qd d) 
        (qs s) 
        (qb b) 
        (h0 "h") 
        (gsg "{g_g}")
      );


((p "\gamma") (pa1 "\gamma~Z^0") (z "Z") (w w) (gs g) (nue "{\nu}_e") (numu 
"{\nu}_{\mu}") (nut "{\nu}_{\tau}") (ef e) (mu "\mu") (tau "\tau") (qu u) (qc c)
(qt t) (qd d) (qs s) (qb b) (h0 "h") (gsg "{g_g}"))

 
ptypelist:='( 
        (psi 1) 
        (psi3 1) 
        (psi5 1) 
        (b1 (gs 43) 42) 
        (b0 (gsg 121) (h0 121) 41)
      );


((psi 1) (psi3 1) (psi5 1) (b1 (gs 43) 42) (b0 (gsg 121) (h0 121) 41))
 

latex_s_list:='(
        (theta "{\theta_w}")  
        (xi "{\xi}")  
        (g3 "{g_s}") 
        (s0 "S_0")
        (s1 "V_1")
        (s2 "V_2")
        (s3 "V_3")
        (s4 "V_4")
        (s5 "V_5")
        (s6 "V_6")
        (s7 "V_7")
        (s8 "V_8")
        (sf "F")
        (sf1 "F_{1/2}")
        (sf3 "F_{3/2}")
        (sf5 "F_{5/2}")
        (sf7 "F_{7/2}")
        (sf9 "F_{9/2}")
        (sf11 "F_{11/2}")
        (sf13 "F_{13/2}")
        (sf15 "F_{15/2}")
     );


((theta "{\theta_w}") (xi "{\xi}") (g3 "{g_s}") (s0 "S_0") (s1 "V_1") (s2 "V_2")
(s3 "V_3") (s4 "V_4") (s5 "V_5") (s6 "V_6") (s7 "V_7") (s8 "V_8") (sf "F") (sf1
"F_{1/2}") (sf3 "F_{3/2}") (sf5 "F_{5/2}") (sf7 "F_{7/2}") (sf9 "F_{9/2}") (sf11
"F_{11/2}") (sf13 "F_{13/2}") (sf15 "F_{15/2}"))
 
end;

10: 
     *  ......The latex file model.tex was generated 
     *        Please latex it 
11: 
Time: 0 ms

12: 
Quitting
