Algorithm 691: improving quadpack automatic integration routines