the fusion crosssection is now correct calculated.
added some documentation.

Worked on the surfival function f. Still does not work. This probably due to the very simple integration used.