change the type of integration used for the Class I integration
changed the integration methode to a beter one. The difference between this code and mathematica remains.

Next step is to see if the cross section for charge exchange is maybe off (what I don't expect).