Appendix A

Section A.2, "The functions eta_circ, eta_+, h and h_H"

./integMh a b tol r k

computes the integral \int_a^{b} |Mh(i t)| dt

a,b are positive integers
tol = error tolerance (per step; approximate)
r = initial number of subdivisions of an interval [n,n+1]
k = initial number of terms of Taylor series

Examples:

./integMh 1 5000 1e-11 10000 30

The final lines of the output are:
--------
Integral of |Mh(t)| from 4999 to 5000: [1.3773001e-12,1.3773423e-12]
Total integral from 1 to 5000: [4.0938731863,4.0938731877]
--------
