/* rate from n_TOF  */
/* Be7 + n -> p + Li7 */
a0=   6.8046e+09; a1= -1.9706e+00; a2=  2.0419e+00; a3= -1.0687e+00; a4=  2.7172e-01; a5= -2.6698e-02; 
a6=  1.9627e+08; a7=  2.8894e+07; b0=  2.8108e-01;
f[17] = a0*(1. + a1*pow(T9,1./2.) + a2*T9 + a3*pow(T9,3./2.) + a4*T9*T9 + a5*pow(T9,5./2.)) + 
+ a6*pow(T9/(1.+13.076*T9),3./2.)/pow(T9,3./2.) + a7/pow(T9,3./2.)*exp(-b0/T9);

-- AlbertoMengoni - 2018-02-19