A Mathematica implementation of Double Adaptive Quadrature