Change the second part integration to Dawson function and error function,the performance is improved by 10%.
Add some commets for the Dawson function and Error function in secondPart.c