Numerical algorithms are required when the solutions to mathematical problems cannot be expressed analytically. Adaptive algorithms expend the... Show moreNumerical algorithms are required when the solutions to mathematical problems cannot be expressed analytically. Adaptive algorithms expend the right amount of computational e↵ort to meet the error tolerance–easy problems require less e↵ort and hard problems require more e↵ort. Here we introduce a new globally adaptive trapezoidal rule algorithm satisfies a general error criterion, which includes both absolute and relative error tolerances. Then we derive the computational cost of the new algorithm and the complexity of this problem. M.S. in Applied Mathematics, May 2018 Show less
This thesis investigates how to solve univariate integration problems using numerical methods, including the trapezoidal rule and the Simpson... Show moreThis thesis investigates how to solve univariate integration problems using numerical methods, including the trapezoidal rule and the Simpson's rule. Most existing guaranteed algorithms are not adaptive and require too much a priori information. Most existing adaptive algorithms do not have valid justification for their results. The goal is to create adaptive algorithms utilizing the two above-mentioned methods with guarantees. The classes of integrands studied in this thesis are cones. The algorithms are analytically proved to be a success if the integrand lies in the cone. The algorithms are adaptive and automatically adjust the computational costs based on the integrand values. The lower and upper bounds on the computational costs for both algorithms are derived. The lower bounds on the complexity of the problems are derived as well. By comparing the upper bounds on the computational cost and the lower bounds on the complexity, our algorithms are shown to be asymptotically optimal. Numerical experiments are implemented. Show less