1. This graph is the exponential decay Si28+n total. Theorotical cross section is between 0-1.5 while the experimental should be towards 1.15

2. This the the distribution of l = 0.

3. This the the distribution of l = 1.

4. This is the TALYS output for Si28+n.

5. mygoe.f input code which consists of TALYS parameters (Gg, D, S).

#!/bin/bash
rm -r Si29/
mkdir Si29
for((c=1;c<=1000;c++))
do
./a.out << EOF
28 0                                     # amass spin  
0 1000000                                #min max
0 211419.28125 1.0356e-4                            #ell D S
0.56467 0.112934                                 #<G> SigmaG
0 0                                      #Gx sigma
0                                        #bias

EOF
cp out_all1 Si29/run1_0l_$c
cp out_all2 Si29/run2_0l_$c
rm out*
done
for((c=1;c<=1000;c++))
do
./a.out << EOF
28 0                                     # amass spin  
0 1000000                                #min max
1 84723.37500 5.9809e-4                            #ell D S
0.98687 0.197374                                 #<G> SigmaG
0 0                                      #Gx sigma
0                                        #bias
EOF
cp out_all1 Si29/run1_1l_$c
cp out_all2 Si29/run2_1l_$c
cp out_all3 Si29/run3_1l_$c
cp out_all4 Si29/run4_1l_$c
rm out*
done