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’...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.展开更多
An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions ...An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions with trigonometric and exponential terms satisfying different conditions are employed to generate a number of formulations. Performances of the new schemes are tested against well-known numerical integrators for selected test cases with quite satisfactory results. Convergence and stability issues of the new formulations are not addressed as the treatment of these aspects requires a separate work. The general approach introduced herein opens a wide vista for producing virtually unlimited number of formulations.展开更多
In this paper, a new one-step explicit method of fourth order is derived. The new method is proved to be A-stable and L-stable, and it gives exact results when applied to the test equation y’=λy with Re(λ)【0, Also...In this paper, a new one-step explicit method of fourth order is derived. The new method is proved to be A-stable and L-stable, and it gives exact results when applied to the test equation y’=λy with Re(λ)【0, Also several numerical examples are included.展开更多
In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an alg...In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an algorithm whose accuracy is good for arbitrary e>0 .展开更多
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 solutio...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.展开更多
An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the ex...An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the existence and uniqueness of solutions for the system were derived. Using a fractional predictorcorrector method, a numerical method was presented for the specified system. An example was given to illustrate the obtained results.展开更多
This paper investigates some known difference schemes for the numerical solution to parabolic differential equation with derivative boundary conditions by the fictitious domain method.The stability and convergence in...This paper investigates some known difference schemes for the numerical solution to parabolic differential equation with derivative boundary conditions by the fictitious domain method.The stability and convergence in L ∞ are proven.展开更多
In this paper, the principle techinique of the differentiator method, and some examples using the method to obtain the general solution and special solution of the differential equation are introduced. The essential d...In this paper, the principle techinique of the differentiator method, and some examples using the method to obtain the general solution and special solution of the differential equation are introduced. The essential difference between this method and the others is that by this method special and general solutions can be obtained directly with the operations of the differentor in the differential equation and without the enlightenment of other scientific knowledge.展开更多
Boundary value problem; for third-order ordinary differential equations with turning points are studied as follows : epsilon gamma ' ' + f(x ; epsilon) gamma ' + g(x ; epsilon) gamma ' +h(x ; epsilon) ...Boundary value problem; for third-order ordinary differential equations with turning points are studied as follows : epsilon gamma ' ' + f(x ; epsilon) gamma ' + g(x ; epsilon) gamma ' +h(x ; epsilon) gamma = 0 (- a < x < b, 0 epsilon 1), where f(x ; 0) has several multiple zero points in ( - n, b). the necessary conditions for exhibiting resonance is given, and the uniformly valid asymptotic solutions and the estimations of remainder terms are obtained.展开更多
We emphasized explicitly on the derivation and implementation of a new numerical algorithm scheme which gave stable results that show the applicability of the method. In this paper, we aimed to solve some second order...We emphasized explicitly on the derivation and implementation of a new numerical algorithm scheme which gave stable results that show the applicability of the method. In this paper, we aimed to solve some second order initial value problems of ordinary differential equations and compare the results with the theoretical solution. Using this method to solve some initial value problems of second order ordinary differential equations, we discovered that the results compared favorably with the theoretical solution which led to the conclusion that the new numerical algorithm scheme derived in the research is approximately correct and can be prescribed for any related ordinary differential equations.展开更多
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, robus...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.展开更多
In this paper, we used an interpolation function to derive a Numerical Integrator that can be used for solving first order Initial Value Problems in Ordinary Differential Equation. The numerical quality of the Integra...In this paper, we used an interpolation function to derive a Numerical Integrator that can be used for solving first order Initial Value Problems in Ordinary Differential Equation. The numerical quality of the Integrator has been analyzed to authenticate the reliability of the new method. The numerical test showed that the finite difference methods developed possess the same monotonic properties with the analytic solution of the sampled Initial Value Problems.展开更多
A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obta...A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obtained from other algorithms. It is shown that the suggested algorithm competes strongly with other existing algorithms, both in accuracy and ease of application, while demanding a shorter computation time.展开更多
In this paper, we prove existence and multiplicities of solutions for asymptotically linear ordinary differential equations satisfying Sturm-Liouville boundary value conditions with resonance. Adding assumption H3 tha...In this paper, we prove existence and multiplicities of solutions for asymptotically linear ordinary differential equations satisfying Sturm-Liouville boundary value conditions with resonance. Adding assumption H3 that is similar to (LL) in Theorem 1.1, by index theory and Morse theory, we obtain more nontrivial solutions.展开更多
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 ...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展开更多
By making use of the differential inequalities, in this paper we study the uniqueness of solutions of the two kinds of the singularly perturbed boundary value problems for the nonlinear third order ordinary differenti...By making use of the differential inequalities, in this paper we study the uniqueness of solutions of the two kinds of the singularly perturbed boundary value problems for the nonlinear third order ordinary differential equation with a small parameter ε>0: where i=1, 2; a(?)(ε), β(ε) and γ(ε) are functions defined on (0, ε_o], while ε_o>0 is a constant.This paper is the continuation of our works [4, 6].展开更多
Nonlinear partial differetial equation(NLPDE) is converted into ordinary differential equation(ODE) via a new ansatz.Using undetermined function method,the ODE obtained above is replaced by a set of algebraic equation...Nonlinear partial differetial equation(NLPDE) is converted into ordinary differential equation(ODE) via a new ansatz.Using undetermined function method,the ODE obtained above is replaced by a set of algebraic equations which are solved out with the aid of Mathematica.The exact solutions and solitary solutions of NLPDE are obtained.展开更多
In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarant...In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.展开更多
In this paper, we used an interpolation function with strong trigonometric components to derive a numerical integrator that can be used for solving first order initial value problems in ordinary differential equation....In this paper, we used an interpolation function with strong trigonometric components to derive a numerical integrator that can be used for solving first order initial value problems in ordinary differential equation. This numerical integrator has been tested for desirable qualities like stability, convergence and consistency. The discrete models have been used for a numerical experiment which makes us conclude that the schemes are suitable for the solution of first order ordinary differential equation.展开更多
A framework to obtain numerical solution of the fractional partial differential equation using Bernstein polynomials is presented. The main characteristic behind this approach is that a fractional order operational ma...A framework to obtain numerical solution of the fractional partial differential equation using Bernstein polynomials is presented. The main characteristic behind this approach is that a fractional order operational matrix of Bernstein polynomials is derived. With the operational matrix, the equation is transformed into the products of several dependent matrixes which can also be regarded as the system of linear equations after dispersing the variable. By solving the linear equations, the numerical solutions are acquired. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples are provided to show that the method is computationally efficient.展开更多
文摘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.
文摘An entirely new framework is established for developing various single- and multi-step formulations for the numerical integration of ordinary differential equations. Besides polynomials, unconventional base-functions with trigonometric and exponential terms satisfying different conditions are employed to generate a number of formulations. Performances of the new schemes are tested against well-known numerical integrators for selected test cases with quite satisfactory results. Convergence and stability issues of the new formulations are not addressed as the treatment of these aspects requires a separate work. The general approach introduced herein opens a wide vista for producing virtually unlimited number of formulations.
文摘In this paper, a new one-step explicit method of fourth order is derived. The new method is proved to be A-stable and L-stable, and it gives exact results when applied to the test equation y’=λy with Re(λ)【0, Also several numerical examples are included.
文摘In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an algorithm whose accuracy is good for arbitrary e>0 .
文摘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.
基金National Natural Science Foundation of China(No.11371087)
文摘An initial value problem was considered for a coupled differential system with multi-term Caputo type fractional derivatives. By means of nonlinear alternative of Leray-Schauder and Banach contraction principle,the existence and uniqueness of solutions for the system were derived. Using a fractional predictorcorrector method, a numerical method was presented for the specified system. An example was given to illustrate the obtained results.
文摘This paper investigates some known difference schemes for the numerical solution to parabolic differential equation with derivative boundary conditions by the fictitious domain method.The stability and convergence in L ∞ are proven.
文摘In this paper, the principle techinique of the differentiator method, and some examples using the method to obtain the general solution and special solution of the differential equation are introduced. The essential difference between this method and the others is that by this method special and general solutions can be obtained directly with the operations of the differentor in the differential equation and without the enlightenment of other scientific knowledge.
文摘Boundary value problem; for third-order ordinary differential equations with turning points are studied as follows : epsilon gamma ' ' + f(x ; epsilon) gamma ' + g(x ; epsilon) gamma ' +h(x ; epsilon) gamma = 0 (- a < x < b, 0 epsilon 1), where f(x ; 0) has several multiple zero points in ( - n, b). the necessary conditions for exhibiting resonance is given, and the uniformly valid asymptotic solutions and the estimations of remainder terms are obtained.
文摘We emphasized explicitly on the derivation and implementation of a new numerical algorithm scheme which gave stable results that show the applicability of the method. In this paper, we aimed to solve some second order initial value problems of ordinary differential equations and compare the results with the theoretical solution. Using this method to solve some initial value problems of second order ordinary differential equations, we discovered that the results compared favorably with the theoretical solution which led to the conclusion that the new numerical algorithm scheme derived in the research is approximately correct and can be prescribed for any related ordinary differential equations.
文摘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.
文摘In this paper, we used an interpolation function to derive a Numerical Integrator that can be used for solving first order Initial Value Problems in Ordinary Differential Equation. The numerical quality of the Integrator has been analyzed to authenticate the reliability of the new method. The numerical test showed that the finite difference methods developed possess the same monotonic properties with the analytic solution of the sampled Initial Value Problems.
文摘A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obtained from other algorithms. It is shown that the suggested algorithm competes strongly with other existing algorithms, both in accuracy and ease of application, while demanding a shorter computation time.
文摘In this paper, we prove existence and multiplicities of solutions for asymptotically linear ordinary differential equations satisfying Sturm-Liouville boundary value conditions with resonance. Adding assumption H3 that is similar to (LL) in Theorem 1.1, by index theory and Morse theory, we obtain more nontrivial solutions.
文摘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
基金Project supported by the National Natural Science Foundation of China.
文摘By making use of the differential inequalities, in this paper we study the uniqueness of solutions of the two kinds of the singularly perturbed boundary value problems for the nonlinear third order ordinary differential equation with a small parameter ε>0: where i=1, 2; a(?)(ε), β(ε) and γ(ε) are functions defined on (0, ε_o], while ε_o>0 is a constant.This paper is the continuation of our works [4, 6].
基金Supported by the Natural Science Foundation of Zhejiang Province(1 0 2 0 3 7)
文摘Nonlinear partial differetial equation(NLPDE) is converted into ordinary differential equation(ODE) via a new ansatz.Using undetermined function method,the ODE obtained above is replaced by a set of algebraic equations which are solved out with the aid of Mathematica.The exact solutions and solitary solutions of NLPDE are obtained.
基金Supported by Russian Fund of Fund amental Investigations(Pr.990101064)and Russian Minister of Educatin
文摘In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.
文摘In this paper, we used an interpolation function with strong trigonometric components to derive a numerical integrator that can be used for solving first order initial value problems in ordinary differential equation. This numerical integrator has been tested for desirable qualities like stability, convergence and consistency. The discrete models have been used for a numerical experiment which makes us conclude that the schemes are suitable for the solution of first order ordinary differential equation.
基金supported by the Natural Science Foundation of Hebei Province under Grant No.A2012203407
文摘A framework to obtain numerical solution of the fractional partial differential equation using Bernstein polynomials is presented. The main characteristic behind this approach is that a fractional order operational matrix of Bernstein polynomials is derived. With the operational matrix, the equation is transformed into the products of several dependent matrixes which can also be regarded as the system of linear equations after dispersing the variable. By solving the linear equations, the numerical solutions are acquired. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples are provided to show that the method is computationally efficient.