Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.

Numerical Treatment of Initial Value Problems of Nonlinear Ordinary Differential Equations by Duan-Rach-Wazwaz Modified Adomian Decomposition Method

We employ the Duan-Rach-Wazwaz modified Adomian decomposition method for solving initial value problems for the systems of nonlinear ordinary differential equations numerically. In order to confirm practicality, robustness and reliability of the method, we compare the results from the modified Adomian decomposition method with those from the MATHEMATICA solutions and also from the fourth-order Runge Kutta method solutions in some cases. Furthermore, we apply Padé approximants technique to improve the solutions of the modified decomposition method whenever the exact solutions exist.

Solution of Laguerre’s Differential Equations via Modified Adomian Decomposition Method

This paper presents a technique for obtaining an exact solution for the well-known Laguerre’s differential equations that arise in the modeling of several phenomena in quantum mechanics and engineering. We utilize an efficient procedure based on the modified Adomian decomposition method to obtain closed-form solutions of the Laguerre’s and the associated Laguerre’s differential equations. The proposed technique makes sense as the attitudes of the acquired solutions towards the neighboring singular points are correctly taken care of.

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.

THE NUMERICAL STABILITY OF THE BLOCK θ-METHODS FOR DELAY DIFFERENTIAL EQUATIONS

This paper focuses on the numerical stability of the block θ methods adapted to differential equations with a delay argument. For the block θ methods, an interpolation procedure is introduced which leads to the numerical processes that satisfy an important asymptotic stability condition related to the class of test problems y′(t)=ay(t)+by(t-τ) with a,b∈C, Re(a)<-|b| and τ>0. We prove that the block θ method is GP stable if and only if the method is A stable for ordinary differential equations. Furthermore, it is proved that the P and GP stability are equivalent for the block θ method.

The θ－Methods in Numerical Solution of Systems of Differential Equations with Two Delay Terms

This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.

Numerical Study of Fractional Differential Equations of Lane-Emden Type by Method of Collocation

Lane-Emden differential equations of order fractional has been studied.Numerical solution of this type is considered by collocation method. Some of examples are illustrated. The comparison between numerical and analytic methods has been introduced.

Existence of Periodic Solutions for Odd Order Ordinary Differential Equations via the Homotopy Method

This paper deals with the problems of finding periodic solutions for the third order ordinary differential equations of the form (1) where T is a fixed positive number and f satisfies some additional conditions which will be stated later.The periodicity problem has been one of main topics in the qualitative theory of ordinary

SECOND-ORDER ACCURATE DIFFERENCE METHOD FOR THE SINGULARLY PERTURBED PROBLEM OF FOURTH-ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.

On the Application of Adomian Decomposition Method to Special Equations in Physical Sciences

The current manuscript makes use of the prominent iterative procedure, called the Adomian Decomposition Method (ADM), to tackle some important special differential equations. The equations of curiosity in this study are the singular equations that arise in many physical science applications. Thus, through the application of the ADM, a generalized recursive scheme was successfully derived and further utilized to obtain closed-form solutions for the models under consideration. The method is, indeed, fascinating as respective exact analytical solutions are accurately acquired with only a small number of iterations.

THE STABILITY OF LINEAR MULTISTEP METHODS FOR SYSTEMS OF DELAY DIFFERENTIAL EQUATIONS

This paper deals with the numerical solution of initial value problems for systems of differential equations with a delay argument. The numerical stability of a linear multistep method is investigated by analysing the solution of the lest equation y’(t)=Ay(t) + By(1-t),where A,B denote constant complex N×N-matrices,and t】0.We investigate carefully the characterization of the stability region.

An Effective Numerical Method for the Solution of a Stochastic Coronavirus(2019-nCovid)Pandemic Model

Nonlinear stochastic modeling plays a significant role in disciplines such as psychology,finance,physical sciences,engineering,econometrics,and biological sciences.Dynamical consistency,positivity,and boundedness are fundamental properties of stochastic modeling.A stochastic coronavirus model is studied with techniques of transition probabilities and parametric perturbation.Well-known explicit methods such as Euler Maruyama,stochastic Euler,and stochastic Runge–Kutta are investigated for the stochastic model.Regrettably,the above essential properties are not restored by existing methods.Hence,there is a need to construct essential properties preserving the computational method.The non-standard approach of finite difference is examined to maintain the above basic features of the stochastic model.The comparison of the results of deterministic and stochastic models is also presented.Our proposed efficient computational method well preserves the essential properties of the model.Comparison and convergence analyses of the method are presented.

Solving Large Scale Nonlinear Equations by a New ODE Numerical Integration Method

In this paper a new ODE numerical integration method was successfully applied to solving nonlinear equations. The method is of same simplicity as fixed point iteration, but the efficiency has been significantly improved, so it is especially suitable for large scale systems. For Brown’s equations, an existing article reported that when the dimension of the equation N = 40, the subroutines they used could not give a solution, as compared with our method, we can easily solve this equation even when N = 100. Other two large equations have the dimension of N = 1000, all the existing available methods have great difficulties to handle them, however, our method proposed in this paper can deal with those tough equations without any difficulties. The sigularity and choosing initial values problems were also mentioned in this paper.

PROJECTION METHODS AND APPROXIMATIONS FOR ORDINARY DIFFERENTIAL EQUATIONS

A formulation of a differential equation as projection and fixed point pi-Mem alloivs approximations using general piecnvise functions. We prone existence and uniqueness of the up proximate solution* convergence in the L2 norm and nodal supercnnvergence. These results generalize those obtained earlier by Hulme for continuous piecevjise polynomials and by Delfour-Dubeau for discontinuous pieceuiise polynomials. A duality relationship for the two types of approximations is also given.

The Stability of Runge－Kutta Methods for Systems of Delay Differential Equations

ＴｈｅＳｔａｂｉｌｉｔｙｏｆＲｕｎｇｅ－ＫｕｔｔａＭｅｔｈｏｄｓｆｏｒＳｙｓｔｅｍｓｏｆＤｅｌａｙＤｉｆｆｅｒｅｎｔｉａｌＥｑｕａｔｉｏｎｓ￥ＷＡＮＧＸｉａｏｂｉａｏ；ＬＩＵＭｉｎｇｚｈｕ(王晓彪)(刘明珠)（Ｄｅｐｔ．ｏｆＭａｔｈｅｍａｔｉｃｓ，Ｈａｒ...

A Conceptual Numerical Model of the Wave Equation Using the Complex Variable Boundary Element Method

In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.

Production of the Reduction Formula of Seventh Order Runge-Kutta Method with Step Size Control of an Ordinary Differential Equation

The purpose of the present work is to construct a nonlinear equation system (85 × 53) using Butcher’s Table and then by solving this system to find the values of all set parameters and finally the reduction formula of the Runge-Kutta (7,9) method (7th order and 9 stages) for the solution of an Ordinary Differential Equation (ODE). Since the system of high order conditions required to be solved is too complicated, we introduce a subsystem from the original system where all coefficients are found with respect to 9 free parameters. These free parameters, as well as some others in addition, are adjusted in such a way to furnish more efficient R-K methods. We use the MATLAB software to solve several of the created subsystems for the comparison of our results which have been solved analytically.

Derivation of the Reduction Formula of Sixth Order and Seven Stages Runge-Kutta Method for the Solution of an Ordinary Differential Equation

This paper is describing in detail the way we define the equations which give the formulas in the methods Runge-Kutta 6th order 7 stages with the incorporated control step size in the numerical solution of Ordinary Differential Equations (ODE). The purpose of the present work is to construct a system of nonlinear equations and then by solving this system to find the values of all set parameters and finally the reduction formula of the Runge-Kutta (6,7) method (6th order and 7 stages) for the solution of an Ordinary Differential Equation (ODE). Since the system of high order conditions required to be solved is complicated, all coefficients are found with respect to 7 free parameters. These free parameters, as well as some others in addition, are adjusted in such a way to furnish more efficient R-K methods. We use the MATLAB software to solve several of the created subsystems for the comparison of our results which have been solved analytically. Some examples for five different choices of the arbitrary values of the systems are presented in this paper.

Computing Differentially Rotating Neutron Stars Obeying Realistic Equations of State by Using Hartle’s Perturbation Method

In this paper, we use the well-known Hartle’s perturbation method in order to compute models of differentially rotating neutron stars obeying realistic equations of state. In our numerical treatment, we keep terms up to third order in the angular velocity. We present indicative numerical results for models satisfying a particular differential rotation law. We emphasize on computing the change in mass owing to this differential rotation law.

THE STABILITY OF θ-METHODS FOR PANTOGRAPH DELAY DIFFERENTIAL EQUATIONS

This paper deals with the numerical solution of initial value problems for pantograph differential equations with variable delays. We investigate the stability of one leg θ-methods in the numerical solution of these problems. Sufficient conditions for the asymptotic stability of θ-methods are given by Fourier analysis and Ergodic theory.