Related papers: Elementary numerical methods for double integrals
We study two modifications of the trapezoidal product cubature formulae, approximating double integrals over the square domain $[a,b]^2=[a,b]\times [a,b]$. Our modified cubature formulae use mixed type data: except evaluations of the…
We give a new method for numerically solving Abel integral equations of first kind. An estimation for the error is obtained. The method is based on approximations of fractional integrals and Caputo derivatives. Using trapezoidal rule and…
In this paper we present a simple method of deriving a posteriori error equalities and estimates for linear elliptic and parabolic partial differential equations. The error is measured in a combined norm taking into account both the primal…
We present a numerical method for rigorous over-approximation of a reachable set of differential inclusions. The method gives high-order error bounds for single step approximations and a uniform bound on the error over the finite time…
In this paper, we establish several new inequalities for some twice differantiable mappings. Then, we apply these inequalities to obtain new midpoint, trapezoid and perturbed trapezoid rules. Finally, some applications for special means of…
Given a real function $f$ on an interval $[a,b]$ satisfying mild regularity conditions, we determine the number of zeros of $f$ by evaluating a certain integral. The integrand depends on $f, f'$ and $f''$. In particular, by approximating…
In this paper, we propose a numerical method for computing Hadamard finite-part integrals with an integral-power singularity at an endpoint, the part of the divergent integral which is finite as a limiting procedure. In the proposed method,…
In this paper we present a new family of rules for numerical integration. This family has up to half the error of the widely used Newton-Cotes rules when a sufficient number of points is evaluated and also much better numerical stability…
This paper presents two methods for approximating a proper subset of the entries of a Hessian using only function evaluations. These approximations are obtained using the techniques called \emph{generalized simplex Hessian} and…
In our preceding paper, we have proposed an algorithm for obtaining finite-norm solutions of higher-order linear ordinary differential equations of the Fuchsian type [\sum_m p_m (x) (d/dx)^m] f(x) = 0 (where p_m is a polynomial with…
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 single-exponential and a…
In this work we consider the general functional-integral equation: \begin{equation*} y(t) = f\left(t, \int_{a}^{b} k(t,s)g(s,y(s))ds\right), \qquad t\in [a,b], \end{equation*} and give conditions that guarantee existence and uniqueness of…
The aim of this paper is to present an efficient numerical procedure to approximate the generalized Abel's integral equations of the first and second kinds. For this reason, the Taylor polynomials and the collocation method are applied.…
An explicit formula to approximate the diagonal entries of the Hessian is introduced. When the derivative-free technique called \emph{generalized centered simplex gradient} is used to approximate the gradient, then the formula can be…
In this paper we propose methods for computing Fresnel integrals based on truncated trapezium rule approximations to integrals on the real line, these trapezium rules modified to take into account poles of the integrand near the real axis.…
We derive a numerical method, based on operator splitting, to abstract parabolic semilinear boundary coupled systems. The method decouples the linear components which describe the coupling and the dynamics in the bulk and on the surface,…
We consider an initial/boundary value problem for one-dimensional fractional-order parabolic equations with a space fractional derivative of Riemann-Liouville type and order $\alpha\in (1,2)$. We study a spatial semidiscrete scheme with the…
The main purpose of this article is to facilitate the implementation of space-time finite element methods in four-dimensional space. In order to develop a finite element method in this setting, it is necessary to create a numerical…
In this paper, we introduce and analyze arbitrarily high-order quadrature rules for evaluating the two-dimensional singular integrals of the forms \begin{align} I_{i,j} = \int_{\mathbb{R}^2}\phi(x)\frac{x_ix_j}{|x|^{2+\alpha}} \d x, \quad…
A class of numerical quadrature rules is derived, with equally-spaced nodes, and unit weights except at a few points at each end of the series, for which "corrections" (not using any further information about the integrand) are added to the…