... dftintwhich uses dftcor to compute equation (13.9.1) for general a, b, ω,andh(t). Several pointswithin this program bear mentioning: The parameters M and NDFT correspond to M and Nin the above discussion. ... inputvaluesof ω,∆,a,b ,and an array with the eight values h0, ,h3,hM−3, ,hM,it returns therealand imaginary parts of the endpointcorrections in equation (13.9.13), and thefactor W(θ). The code ... cancellationsto high powers of θ. It istherefore necessarytocompute the right-hand sides in double precision, even whenthe corrections are desired only to single precision.It is also necessary to use...