Change the pmode iteration loop from file reading to generation within the subrountines, in order to improve the speed further
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