increase default quadgk tolerance to sqrt(eps), similar to Matlab's; 100*eps is far too agressive in the case of badly-behaved integrands