Discrete Temimi-Ansari method for solving a class of stochastic nonlinear differential equations

: In this paper, a numerical method to solve a class of stochastic nonlinear differential equations is introduced. The proposed method is based on the Temimi-Ansari method. The special states of the four systems are studied to show that the proposed method is efﬁcient and applicable. These systems are stochastic Langevin's equation, Ginzburg-Landau equation, Davis-Skodje, and Brusselator systems. The results clarify the accuracy and efficacy of the presented new method with no need for any restrictive assumptions for nonlinear terms.

The main contribution of the current paper is to implement a numerical technique based on Temimi and Ansari (TAM) [25] to solve the following model of stochastic nonlinear differential equations: ( , )+( ( , )+ (t)) (t), (0)= du F t u G t u f n u a dt  = (1) where (t) u is the unknown function, represents the independent variable, . The Temimi and Ansari method (TAM) has been successfully applied to solve many classes of differential equations, such as nonlinear second-order multipoint boundary value problems [25], ordinary differential equations [26], Fokker-Planck's equations [27], and korteweg-de Varies equations [28]. In this paper, we modified the Temimi-Ansari method called the discrete Temimi-Ansari method DTAM to handle the stochastic nonlinear differential Eq (1). The proposed method provides a prodigious performance by merging the traditional TAM with the finite difference numerical scheme. To test the proposed method, we solve four stochastic models. These models are stochastic Langevin's equation, Ginzburg-Landau equation, Davis-Skodje, and Brusselator systems. These problems illustrate the accuracy of the proposed method as compared with the stochastic Runge Kutta method (SRK), and semi-analytic methods.
This work is structured as follows: Section 2, introduces the discrete Temimi-Ansari method (DTAM). In the third section, the stochastic models were solved using the proposed method. Finally, conclusions are presented in Section 4.

Discrete Temimi-Ansari method and convergence analysis
Consider the following differential equation in the form is with initial conditions: where and represent the linear and nonlinear operators, respectively, and ( ) represents the inhomogeneous term. The Temimi-Ansari method was used to solve differential Eq (3) as follows: First consider the initial approximate function 0 (t), which is the solution of the following initial value problem 0 00 To find the next approximate function 1 (t), the following problem must be solved The n th approximate functions ) (t u n can be evaluated in the same way. In this method, it is very serious to note that each of () n ut is separately a solution to model (3a). We document that this scheme is easy to implement and has the distinction that each solution is a refinement of the previous iteration. Sequential solutions must be verified versus the previous iteration to confirm the convergence of solutions. The analytical solution converges to the exact solution as the number of iterations increases. Depending on this, an analytical solution can be obtained with good agreement with the exact solution as An error analysis of a single nonlinear ordinary differential equation, generalization of these theorems applied to systems of differential equations, and a complete study of convergence criteria for the TAM iterative scheme are presented in [26,29,30].
To supply the convergence analysis for the proposed method in this paper, we begin by submitting the following steps for the proposed approximate analytical scheme: is defined as where the term () n ut is the analytical solution of the TAM. According to these criteria, appropriate conditions for the convergence of this technique are studied with the following theorems. Proof. Fruition [29,30].
This theory is a special case of fixed-point theory and it is sufficient to prove the convergence of TAM.
The solution obtained by TAM converges to the exact solution provided that: 0 < < 1   such that ut The TAM approach can be used to solve the differential Eq (1) with a random function excitation. However, only a few iterations can be calculated due to the difficulty of integrating random functions. We adopted discrete Temimi-Ansari (DTAM) to solve stochastic nonlinear differential Eq (1) as follows: Suppose: where to discretize the initial approximate function ) and by using The next discrete approximate function and the n th discrete approximate functions ) ( 1 + i n t u can be computed as follows The Wiener process ) (t w is a random function, so the solution must be based on k runs with different samples of the Wiener process ) (t w . Then the iterative Eqs (13) and (15) can be written in the following form: The time step h should be used to warrant convergence of systems (16a)-(16c). Applying the convergence standard of the fixed-point iteration by differentiating the right-hand side of Eq (16b) concerning 0,k u we obtain then we have: The condition is a sufficient condition for the time step used in FDM for convergence. Due to similarity in the other two Eqs (16a) and (16c), the same condition can be also used in the numerical FDM approximations. Finally, by taking the mean and the variance of the solution sequences { ,1 , ,2 , ,3 , … . , , } we obtain the mean and the variance of the solution. This elegant combination of the traditional TAM and finite difference numerical scheme produces a fast and effective tool in handling stochastic nonlinear differential equations and this is one of its uncountable merits. Thus, we simply calculate many iterations without time consumption in comparison with traditional methods such as HAM, WHEP, and even TAM which are unable to perform many iterations due to the complexity of calculations of the mean and the variance.

Numerical simulation
In this section, four examples are given to illustrate the applicability of the proposed method, and all of them are performed on a computer by using the Mathematica program.

Example 1.
Consider the stochastic Langevin's equation [31] of the form: The constructed iterative scheme (16) for stochastic Langevin's Eq (20) where 0, , 0     . Eq (22) from the theory of superconductivity is derived by Ginzburg and Landau to describe a phase transition in a deterministic sense [32,33]    Example 3. Consider Davis-Skodje system is [33] The chemical reaction system (24) The stable equilibrium point for the nonlinear system (24) is (0, 0) proven in [34]. Figure This nonlinear system is unforced periodic oscillations in certain chemical reactions. We simulate this system when 1.9  = and 0.1  = , and the proposed iterative scheme of it is given by

Discussion
The TAM approach can solve the differential Eq (1) with a random function excitation. However, only a few iterations can be calculated because of the stiffness of integrating random functions. We suggested discrete Temimi-Ansari (DTAM) to solve the stochastic nonlinear differential equation. We establish an unusual stochastic method to approximate a wide class of stochastic differential equations by combining the classical TAM with the finite difference numerical scheme.
The ability of the proposed scheme has been illustrated by several nonlinear stochastic differential equations. In particular, Figures 1-5 show that the proposed method is able to solve the nonlinear stochastic Langevin's equation (20) and . Also, calculating the amount of the convergence parameter for Eqs (20) and (22) as in (11) is shown in Tables 1 and 2. These tables refer to the values of less than one and verify that the proposed solutions are convergent. In Figures 6 and 7, we examine the behavior of the stochastic Davis-Skodje system (24) and the Brusselator system (26). It can be seen that the proposed method solutions tend to the solution (0, 0). The results reveal that the proposed method is reliable and can be applied to stiff stochastic problems in applied sciences. This wonderful scheme produces a fast and functional approach in treating the stochastic nonlinear differential equations and this is one of its uncountable advantages. The classical semi-analytic method needs a huge time and effort to reach an acceptable solution because of the complexity of mean and variance solutions. Our proposed method beats this shortage simply. We can expand the above analysis to many random parameters besides the noise enjoined in the previous examples.

Conclusions
In some cases, the traditional TAM is not suitable to solve the stochastic nonlinear differential Eq (1), and we devolved this method to handle this class. The DTAM had been successfully applied to find the solutions of the stochastic nonlinear differential equations. The efficiency and accuracy of the proposed method were demonstrated by solving stochastic Langevin's equation, Ginzburg-Landau equation, Davis-Skodje, and Brusselator systems. Through the figures, it can be seen clearly that the discrete TAM solutions converge when the number of iterations is increased. The motivation of our work is achieved by comparing the stochastic Runge Kutta method with DTAM. Future directions might include solving various complicated stochastic linear and nonlinear differential equations, partial stochastic differential equations, fractional differential equations, fractional stochastic differential equations, delayed fractional stochastic differential equations, and fuzzy differential equations.