Exponentially Convergent Trapezoidal Rules to Approximate Fractional Powers of Operators

In this paper we are interested in the approximation of fractional powers of self-adjoint positive operators. Starting from the integral representation of the operators, we apply the trapezoidal rule combined with a double-exponential transform of the integrand function. In this work we show how to improve the existing error estimates for the scalar case and also extend the analysis to operators. We report some numerical experiments to show the reliability of the estimates obtained.