Numerical quadrature methods for integrals of singular periodic functions and their application to singular and weakly singular integral equations

High accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Such periodic equations are used in the solution of planar elliptic boundary value problems, elasticity, potential theory, conformal mapping, boundary element methods, free surface flows, etc. The use of the quadrature methods is demonstrated with numerical examples

Multi-step Richardson-Romberg Extrapolation: Remarks on Variance Control and complexity

We propose a multi-step Richardson-Romberg extrapolation method for the computation of expectations $E f(X_{_T})$ of a diffusion $(X_t)_{t\in [0,T]}$ when the weak time discretization error induced by the Euler scheme admits an expansion at an order $R\ge 2$. The complexity of the estimator grows as $R^2$ (instead of $2^R$) and its variance is asymptotically controlled by considering some consistent Brownian increments in the underlying Euler schemes. Some Monte carlo simulations carried with path-dependent options (lookback, barriers) which support the conjecture that their weak time discretization error also admits an expansion (in a different scale). Then an appropriate Richardson-Romberg extrapolation seems to outperform the Euler scheme with Brownian bridge.Comment: 28 pages, \`a para\^itre dans Monte Carlo Methods and Applications Journa

Efficient Numerical Methods for Pricing American Options under Lévy Models

Two new numerical methods for the valuation of American and Bermudan options are proposed, which admit a large class of asset price models for the underlying. In particular, the methods can be applied with Lévy models that admit jumps in the asset price. These models provide a more realistic description of market prices and lead to better calibration results than the well-known Black-Scholes model. The proposed methods are not based on the indirect approach via partial differential equations, but directly compute option prices as risk-neutral expectation values. The expectation values are approximated by numerical quadrature methods. While this approach is initially limited to European options, the proposed combination with interpolation methods also allows for pricing of Bermudan and American options. Two different interpolation methods are used. These are cubic splines on the one hand and a mesh-free interpolation by radial basis functions on the other hand. The resulting valuation methods allow for an adaptive space discretization and error control. Their numerical properties are analyzed and, finally, the methods are validated and tested against various single-asset and multi-asset options under different market models

Romberg Type Cubature over Arbitrary Triangles

We develop an extrapolation algorithm for numerical integration over arbitary non-standard triangles in IR², which parallels the well-known univariate Romberg method. This is done by a suitable generalization of the trapezoidal rule over triangles, for which we can prove the existence of an asymptotic expansion. Our approach relies mainly on two ideas: The use of barycentric coordinates and the interpretation of the trapezoidal rule as the integral over an interpolating linear spline function. Since our method works for arbitrary triangles, it yields - via triangulation - a method for cubature over arbitrary, possibly non-convex, polygon regions in IR². Moreover, also numerical integration over convex polyhedra in IR d, d > 2 , can be accomplished without difficulties. Numerical examples show the stability and efficiency of the algorithm

Concerning periodic solutions to non-linear discrete Volterra equations with finite memory

In this paper we discuss the existence of periodic solutions of discrete (and discretized) non-linear Volterra equations with finite memory. The literature contains a number of results on periodic solutions of non linear Volterra integral equations with finite memory, of a type that arises in biomathematics. The “summation” equations studied here can arise as discrete models in their own right but are (as we demonstrate) of a type that arise from the discretization of such integral equations. Our main results are in two parts: (i) results for discrete equations and (ii) consequences for quadrature methods applied to integral equations. The first set of results are obtained using a variety of fixed point theorems. The second set of results address the preservation of properties of integral equations on discretizing them. An expository style is adopted and examples are given to illustrate the discussion