The integrant of f(r) is now working and giving the same results as the mathematica file.

An error was found in the values for the cross section. Some values where missing some digits
The surfival function is now the same as the mathematica code.
The function A is now working