In this paper, we present a fast numerical strategy for filtering stochastic differential equations with multiscale features. This method is designed such that it does not violate the practical linear observability condition and, more importantly, it does not require the computationally expensive cross correlation statistics between multiscale variables that are typically needed in standard fil...

In this paper, a new difference scheme based on C1-quintic splines is derived for the numerical solution of the stiff delay differential equations. Convergence results shows that the methods have a convergence of order five. Moreover, the stability analysis properties of these methods have been studied. Finally, numerical results illustrating the behavior of the methods when faced with some dif...

A class of ^-stable advanced multistep methods is derived for the numerical solution of initial value problems in ordinary differential equations. The methods, of all orders of accuracy up to ten, only require values of y' and are self starting. Results for the asymptotic behaviour of the discretization error and for estimating local truncation error are also obtained. The practical implementat...

We consider the combustion equation as one of the candidates from the class of stiff ordinary differential equations. A solution over a length of time that is inversely proportional to δ > 0 (where δ > 0 is a small disturbance of the pre-ignition state) is sought. This problem has a transient at the midpoint of the integration interval. The solution changes from being non-stiff to stiff, and af...

In this paper an analysis will be presented for the ADI (alternating direction implicit) method of Peaceman and Rachford applied to initial-boundary value problems for partial differential equations in two space dimensions. We shall use the method of lines approach. Motivated by developments in the field of stiff nonlinear ordinary differential equations, our analysis will focus on problems whe...

Abstract. Space discretization of some time-dependent partial differential equations gives rise to stiff systems of ordinary differential equations. In this case, implicit methods should be used and therefore, in general, nonlinear systems must be solved. The solutions to these systems are approximated by iterative schemes and, in order to obtain an efficient code, good initializers should be u...

In this paper, an efficient modification is introduced into the well-known homotopy perturbation method which proves very effective to control the convergence region of approximate solution. By using this scheme, explicit exact solution is calculated in the form of a convergent power series with easily computable components. The proposed algorithm is tested on some neutral stiff functional-diff...

We introduce a solver for stiff ordinary differential equations (ODEs) that is based on the deferred correction scheme for the corresponding Picard integral equation. Our solver relies on the assumption that the solution can be accurately represented by a combination of carefully selected complex exponentials. The solver’s accuracy and stability rely on the computation of highly accurate quadra...

In this paper we study the behavior of a family of implicit numerical methods applied to stochastic differential equations with multiple time scales. We show by a combination of analytical arguments and numerical examples that implicit methods in general fail to capture the effective dynamics at the slow time scale. This is due to the fact that such implicit methods cannot correctly capture non...