Numerical solution of functional integral equations by the variational iteration method

In the present article, we apply the variational iteration method to obtain the numerical solution of the functional integral equations. This method does not need to be dependent on linearization, weak nonlinearity assumptions or perturbation theory. Application of this method in finding the approximate solution of some examples confirms its validity. The results seem to show that the method is very effective and convenient for solving such equations.     

Numerical solution of scattering problems using pseudo-differential impedance operators

Direct numerical solutions of scattering problems based on boundary-integral equations are computationally costly at high frequencies. A numerical method is presented that efficiently computes accurate approximations to unknown surface quantities given known surface data (an approximate Dirichlet-to-Neumann map). The method is based on a pseudo-differential impedance operator (PIO) numerically implemented using rational approximations. An example of a PIO is the so-called on-surface radiation condition (OSRC) method. For a convex obstacle, the method can be viewed as applying a parabolic equation directly on the surface of a scatterer. In contrast to past OSRCs, the use of rational approximations provides accuracy for wide scattering angles which is needed for grazing angles of incidence. Generalization to impedance operators for two-dimensional acoustic scatterers with concave parts is presented. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A degenerate kernel method for eigenvalue problems of compact integral operators

We consider the approximation of eigenfunctions of a compact integral operator with a smooth kernel by a degenerate kernel method. By interpolating the kernel of the integral operator in both the variables, we prove that the error bounds for eigenvalues and for the distance between the spectral subspaces are of the orders h ^2r and h ^r respectively. By iterating the eigenfunctions we show that the error bounds for eigenfunctions are of the orders h ^2r . We give the numerical results.      

Dominant eigenvalue problem for positive integral operators and its solution by Monte Carlo method

In this paper, a method of numerical solution to the dominant eigenvalue problem for positive integral operators is presented. This method is based on results of the theory of positive operators developed by Krein and Rutman. The problem is solved by Monte Carlo method constructing random variables in such a way that differences between results obtained and the exact ones would be arbitrarily small. Some numerical results are shown.     

Numerical solution of integral equations

Summary A new method for the solution of integral equations is presented. The method is based on direct approximation of Dirac's delta operator by linear combination of integral operators. This avoids some pitfalls which arise in more conventional numerical procedures for integral equations.Research sponsored by the U. S. Atomic Energy Commission under contract with the Union Carbide Corporation. the Union Carbide Corporation.     

Numerical solution of inverse spectral problems for Sturm-Liouville operators with discontinuous potentials

We consider Sturm-Liouville differential operators on a finite interval with discontinuous potentials having one jump. As the main result we obtain a procedure of recovering the location of the discontinuity and the height of the jump. Using our result, we apply a generalized Rundell-Sacks algorithm of Rafler and Böckmann for a more effective reconstruction of the potential and present some numerical examples.     

Numerical solution of Fredholm integral equations of the second kind by using integral mean value theorem II. High dimensional problems

In this work, we generalize the numerical method discussed in [Z. Avazzadeh, M. Heydari, G.B. Loghmani, Numerical solution of Fredholm integral equations of the second kind by using integral mean value theorem, Appl. math. modelling, 35 (2011) 2374–2383] for solving linear and nonlinear Fredholm integral and integro-differential equations of the second kind. The presented method can be used for solving integral equations in high dimensions. In this work, we describe the integral mean value method (IMVM) as the technical algorithm for solving high dimensional integral equations. The main idea in this method is applying the integral mean value theorem. However the mean value theorem is valid for multiple integrals, we apply one dimensional integral mean value theorem directly to fulfill required linearly independent equations. We solve some examples to investigate the applicability and simplicity of the method. The numerical results confirm that the method is efficient and simple.     

Numerical solution of the three-dimensional Dirichlet problem for inhomogeneous media by the method of integral equations

We consider the three-dimensional Dirichlet problem for equations of elliptic type in inhomogeneous media. The problem can be reduced to a system of loaded Fredholm integral equations of the second kind over the volume. We prove the uniqueness of a classical solution of the problem. We suggest a numerical solution algorithm of iterative type. An example of the numerical solution of the problem is considered, and the convergence of the iterative procedure is demonstrated numerically.     

Numerical solution of the abel integral equation

A numerical method for the solution of the Abel integral equation is presented. The known function is approximated by a sum of Chebyshev polynomials. The solution can then be expressed as a sum of generalized hypergeometric functions, which can easily be evaluated, using a simple recurrence relation.     

Numerical solution of some three-dimensional problems of elasticity theory by Vekua's method

We consider the solution of the problem of elastic equilibrium of a three-dimensional orthotropic plate in the absence of displacements on the end surfaces under the action of forces applied to the lateral surfaces. The solution of the original problem by Vekua's method is reduced to the solution of a recursive sequence of two-dimensional problems. A numerical solution of these problems is obtained by computer using the finite-difference method. The effect of the number of Legendre polynomials on the accuracy with which the boundary conditions are satisfied is investigated.Translated from Vychislitel'naya i Prikladnaya Matematika, No. 59, pp. 77–84, 1986.     

Numerical solution by the quasi-reversibility method of unstable problems for the first-order evolution equation

Prior bounds are derived on the solution of the perturbed problem in different versions of the quasi-reversibility method used for approximate solution of unstable problems for first-order evolution equations. An example of such a problem is provided by the problem backward in time for the equation of heat conduction. Approximate solution of perturbed problems by difference methods is considered. The investigation of the difference schemes of the quasi-reversibility method relies on the general theory of p-stability of difference schemes. Specific features of solution of problems with non-self-adjoint operators are considered. Efficient difference schemes are constructed for multidimensional problems.Translated from Matematicheskoe Modelirovanie i Reshenie Obratnykh Zadach. Matematicheskoi Fiziki, pp. 93–124, 1993.     

Numerical solution of eigentuple-eigenvector problems in Hilbert spaces by a gradient method

Summary A gradient technique previously developed for computing the eigenvalues and eigenvectors of the general eigenproblemAx=Bx is generalized to the eigentuple-eigenvector problem . Among the applications of the latter are (1) the determination of complex (,x) forAx=Bx using only real arithmetic, (2) a 2-parameter Sturm-Liouville equation and (3) -matrices. The use of complex arithmetic in the gradient method is also discussed. Computational results are presented.This research was partially supported by NSF Grants MPS74-13332 and MCS76-09172     

Numerical solution of Volterra integral equations by spline functions

A method for obtaining spline function approximations to solutions of non-singular Volterra equations of the second kind is presented. Convergence results as well as numerical examples comparing the method with known techniques are given.     

Numerical solution of stationary viscous incompressible fluid problems by the method of fictitious domains

We consider some issues of numerical implementation of the fictitious domain method for viscous incompressible fluid problems. Plane stationary problems are solved by successive approximations in the nonlinearity. Plane heat convection problems in the Boussinesq approximation are also considered. Solution examples of some specimen problems are presented.Translated from Vychislitel'naya Matematika i Matematicheskoe Obespechenie EVM, pp. 255–262, 1985.     

Numerical solution of Volterra integral equation

Summary This paper discusses the use of Gregory's formula for numerical integration of Volterra linear integral equations of the second type. The order of magnitude of the truncation error and the asymptotic behavior of this error are obtained by means of recursive inequalities.