1. This graph is fitted P31+n with rebinned to 10. Theorotical cross section is around 2 and the experimental value is 1.74, which relatively match.

2. This the the distribution of l = 0 with bins = 10.

3. This the the distribution of l = 1 with bins = 10.

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

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

#!/bin/bash
rm -r P32/
mkdir P32
for((c=1;c<=1000;c++))
do
./a.out << EOF
31 0.5                                     # amass spin  
0 1000000                                #min max
0 62044.46094 0.8401e-4                            #ell D S
0.69793 0.139586                                 #<G> SigmaG
0 0                                      #Gx sigma
0                                        #bias

EOF
cp out_all1 P32/run1_0l_$c
cp out_all2 P32/run2_0l_$c
rm out*
done
for((c=1;c<=1000;c++))
do
./a.out << EOF
31 0.5                                     # amass spin  
0 1000000                                #min max
1 33466.96484 6.1765e-4                            #ell D S
1.25400 0.2508                                 #<G> SigmaG
0 0                                      #Gx sigma
0                                        #bias
EOF
cp out_all1 P32/run1_1l_$c
cp out_all2 P32/run2_1l_$c
cp out_all3 P32/run3_1l_$c
cp out_all4 P32/run4_1l_$c
rm out*
done