Calculation of expansion coefficients of series in Chebyshev polynomials for a solution to a Cauchy problem

Two approaches are proposed to determine an initial approximation for the coefficients of an expansion of the solution to a Cauchy problem for ordinary differential equations in the form of series in shifted Chebyshev polynomials of the first kind. This approximation is used in an analytical method to solve ordinary differential equations using orthogonal expansions.     

Iterative solutions of nonlinear initial value differential equations in Chebyshev series using Lie series

Summary In this paper Lie series are presented in Chebyshev form and applied to the iterative solution of initial value problems in differential equations. The resulting method, though algebraically complicated, is of theoretical interest as a generalisation of Taylor series methods and iterative Chebyshev methods. The theory of the method is discussed and the solutions of some simple scalar equations are analysed to illustrate the behaviour of the process.     

A Chebyshev interval method for nonlinear dynamic systems under uncertainty

This paper proposes a new interval analysis method for the dynamic response of nonlinear systems with uncertain-but-bounded parameters using Chebyshev polynomial series. Interval model can be used to describe nonlinear dynamic systems under uncertainty with low-order Taylor series expansions. However, the Taylor series-based interval method can only suit problems with small uncertain levels. To account for larger uncertain levels, this study introduces Chebyshev series expansions into interval model to develop a new uncertain method for dynamic nonlinear systems. In contrast to the Taylor series, the Chebyshev series can offer a higher numerical accuracy in the approximation of solutions. The Chebyshev inclusion function is developed to control the overestimation in interval computations, based on the truncated Chevbyshev series expansion. The Mehler integral is used to calculate the coefficients of Chebyshev polynomials. With the proposed Chebyshev approximation, the set of ordinary differential equations (ODEs) with interval parameters can be transformed to a new set of ODEs with deterministic parameters, to which many numerical solvers for ODEs can be directly applied. Two numerical examples are applied to demonstrate the effectiveness of the proposed method, in particular its ability to effectively control the overestimation as a non-intrusive method.     

The Sinc-collocation method for solving the Thomas–Fermi equation

A numerical technique for solving nonlinear ordinary differential equations on a semi-infinite interval is presented. We solve the Thomas–Fermi equation by the Sinc-Collocation method that converges to the solution at an exponential rate. This method is utilized to reduce the nonlinear ordinary differential equation to some algebraic equations. This method is easy to implement and yields very accurate results.     

Weak second order S-ROCK methods for Stratonovich stochastic differential equations

It is well known that the numerical solution of stiff stochastic ordinary differential equations leads to a step size reduction when explicit methods are used. This has led to a plethora of implicit or semi-implicit methods with a wide variety of stability properties. However, for stiff stochastic problems in which the eigenvalues of a drift term lie near the negative real axis, such as those arising from stochastic partial differential equations, explicit methods with extended stability regions can be very effective. In the present paper our aim is to derive explicit Runge–Kutta schemes for non-commutative Stratonovich stochastic differential equations, which are of weak order two and which have large stability regions. This will be achieved by the use of a technique in Chebyshev methods for ordinary differential equations.     

A Note on the Numerical Solution of Quasi-linear Parabolic Partial Differential Equations with Error Estimates

A method is considered for the numerical solution of quasi-linearpartial differential equations. The partial differential equationis reduced to a set of ordinary differential equations usinga Chebyshev series expansion. The exact solution of this setof ordinary differential equations is shown to be the solutionof a perturbed form of the original equation. This enables errorestimates to be found for linear and mildly non-linear problems.     

Numerical approach to differential equations of fractional order

In this paper, the variational iteration method and the Adomian decomposition method are implemented to give approximate solutions for linear and nonlinear systems of differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper presents a numerical comparison between the two methods for solving systems of fractional differential equations. Numerical results show that the two approaches are easy to implement and accurate when applied to differential equations of fractional order.     

Solving systems of fractional differential equations using differential transform method

This paper presents approximate analytical solutions for systems of fractional differential equations using the differential transform method. The fractional derivatives are described in the Caputo sense. The application of differential transform method, developed for differential equations of integer order, is extended to derive approximate analytical solutions of systems of fractional differential equations. The solutions of our model equations are calculated in the form of convergent series with easily computable components. Some examples are solved as illustrations, using symbolic computation. The numerical results show that the approach is easy to implement and accurate when applied to systems of fractional differential equations. The method introduces a promising tool for solving many linear and nonlinear fractional differential equations.     

Extrapolation methods in Lie groups

Summary. The numerical solution of differential equations on Lie groups by extrapolation methods is investigated. The main principles of extrapolation for ordinary differential equations are extended on the general case of differential equations in noncommutative Lie groups. An asymptotic expansion of the global error is given. A symmetric method is given and quadratic asymptotic expansion of the global error is proved. The theoretical results are verified by numerical experiments. Received September 27, 1999 / Revised version received February 14, 2000 / Published online April 5, 2001     

Numerical study of modified Adomian's method applied to Burgers equation

The application of Adomian's decomposition method to partial differential equations, when the exact solution is not reached, demands the use of truncated series. But the solution's series may have small convergence radius and the truncated series may be inaccurate in many regions. In order to enlarge the convergence domain of the truncated series, Padé approximants (PAs) to the Adomian's series solution have been tested and applied to partial and ordinary differential equations, with good results. In this paper, PAs, both in x$x$ and t$t$ directions, applied to the truncated series solution given by Adomian's decomposition technique for Burgers equation, are tested. Numerical and graphical illustrations show that this technique can improve the accuracy and enlarge the domain of convergence of the solution. It is also shown in this paper, that the application of Adomian's method to the ordinary differential equations set arising from the discretization of the spatial derivatives by finite differences, the so-called method of lines, may reduce the convergence domain of the solution's series.     

Asymptotic theory of the global error and some techniques of error estimation

Summary The error of the approximate solution obtained by discretising a functional equation can be shown under certain conditions to possess an asymptotic expansion in terms of some parameter which is usually a representative step-length. We consider the case of two-parameter expansions, which is particularly relevant to parabolic equations. We derive results for the existence of the expansion and for the application of the classical difference correction and of defect correction. The theory is illustrated by the discussion of a simple parabolic problem     

Bilinear biorthogonal expansions and the Dunkl kernel on the real line

We study an extension of the classical Paley–Wiener space structure, which is based on bilinear expansions of integral kernels into biorthogonal sequences of functions. The structure includes both sampling expansions and Fourier–Neumann type series as special cases, and it also provides a bilinear expansion for the Dunkl kernel (in the rank 1 case) which is a Dunkl analogue of Gegenbauer’s expansion of the plane wave and the corresponding sampling expansions. In fact, we show how to derive sampling and Fourier–Neumann type expansions from the results related to the bilinear expansion for the Dunkl kernel.     

Application of Markov’s quadrature in orthogonal expansions

A method of using Markov’s quadrature with a fixed node is proposed to calculate the coefficients of the expansion of a function in a shifted Chebyshev series. Approximation properties of a partial sum of a series with approximate coefficients are considered. This approach can be used to construct some numerical-analytic methods for solving ordinary differential equations.     

An algebra of differential operators and generating functions on the set of univalent functions

With a method close to that of Kirillov [4], we define sequences of vector fields on the set of univalent functions and we construct systems of partial differential equations which have the sequence of the Faber polynomials (F_n) as a solution. Through the Faber polynomials and Grunsky coefficients, we obtain the generating functions for some of the sequences of vector fields.     

Asymptotic summation of power series

Summary. We give an asymptotic expansion in powers of of the remainder , when the sequence has a similar expansion. Contrary to previous results, explicit formulas for the computation of the coefficients are presented. In the case of numerical series (), rigorous error estimates for the asymptotic approximations are also provided. We apply our results to the evaluation of , which generalizes various summation problems appeared in the recent literature on convergence acceleration of numerical and power series. Received April 22, 1997     

A-priori error estimates of Galerkin backward differentiation methods in time-inhomogeneous parabolic problems

Summary Backward differentiation methods up to orderk=5 are applied to solve linear ordinary and partial (parabolic) differential equations where in the second case the space variables are discretized by Galerkin procedures. Using a mean square norm over all considered time levels a-priori error estimates are derived. The emphasis of the results lies on the fact that the obtained error bounds do not depend on a Lipschitz constant and the dimension of the basic system of ordinary differential equations even though this system is allowed to have time-varying coefficients. It is therefore possible to use the bounds to estimate the error of systems with arbitrary varying dimension as they arise in the finite element regression of parabolic problems.     

Orthogonal polynomials and higher order singular Sturm-Liouville systems

In addition to the classic orthogonal polynomials which satisfy second order differential equations, there are a number of orthogonal polynomials which satisfy differential equations of orders four or six. Like the classic sets, they have distributional weight functions, are the eigenfunctions for certain self-adjoint boundary-value problems, and sometimes are involved with indefinite boundary-value problems.The purpose of this survey is to summarize the work of the last decade and to exhibit the state of the art as it now stands. Of particular interest is the development of the theory of singular Sturm-Liouville systems, which is so necessary in order to describe the boundary-value problems associated with these polynomials.     

A series method for solving nonlinear two-point boundary value problems

Summary A new method for solving nonlinear boundary value problems based on Taylor-type expansions generated by the use of Lie series is derived and applied to a set of test examples. A detailed discussion is given of the comparative performance of this method under various conditions. The method is of theoretical interest but is not applicable, in its present form, to real life problems; in particular, because of the algebraic complexity of the expressions involved, only scalar second order equations have been discussed, though in principle systems of equations could be similarly treated. A continuation procedure based on this method is suggested for future investigation.     

Rooted tree analysis of the order conditions of row-type scheme for stochastic differential equations

Numerical schemes for initial value problems of stochastic differential equations (SDEs) are considered so as to derive the order conditions of ROW-type schemes in the weak sense. Rooted tree analysis, the well-known useful technique for the counterpart of the ordinary differential equation case, is extended to be applicable to the SDE case. In our analysis, the roots are bi-colored corresponding to the ordinary and stochastic differential terms, whereas the vertices have four kinds of label corresponding to the terms derived from the ROW-schemes. The analysis brings a transparent way for the weak order conditions of the scheme. An example is given for illustration.     

Numerical solutions for convection–diffusion equation using El-Gendi method

The method of El-Gendi [El-Gendi SE. Chebyshev solution of differential integral and integro-differential equations. J Comput 1969;12:282–7; Mihaila B, Mihaila I. Numerical approximation using Chebyshev polynomial expansions: El-gendi’s method revisited. J Phys A Math Gen 2002;35:731–46] is presented with interface points to deal with linear and non-linear convection–diffusion equations.The linear problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using three-level time scheme.The non-linear problem is reduced to three systems of ordinary differential. Each one of these systems is, then, solved using three-level time scheme. Numerical results for Burgers’ equation and modified Burgers’ equation are shown and compared with other methods. The numerical results are found to be in good agreement with the exact solutions.