Differential equations are mostly prominent in the applications of chemistry, physics and mathematical models related to economic and biological processes, it is classified into two types, i.e. ordinary and partial, differential equations are divided into linear and non-linear, too. and there are many ways to solve this equation, such as Bernstein’s polynomial which plays a prominent role in various branches of mathematics, and it has been used by many researchers to solve integral equations, and approximation theory    . Saad N. AL-Azawi in  used Bernstein polynomials, for solving systems of the one-dimensional volterra linear integral equations of the second kind. Mohammed A. S. and Ishak H. in  solved the nonlinear stiff system of ordinary differential equations by applying Bernstein’s polynomial. Yadollah O. and Sara D. in  used Bernstien’s polynomials application to solve nonlinear fredholm integro-differential equations. As well as Ahmed F. Q. and Ekhlass S. AL-Rawi in  presented modified Bernstein Polynomials with Adomian decomposition method for solving ordinary and Partial differential equations. Adomian decomposition is a method used to solve linear and non-linear equations, for example, partial differential equations, algebraic, differential, and many systems, integral-differential equations, delay differential equations, … etc., this method gives approximate solutions and accurate calculations of linear and nonlinear differential and integral equations, and the solution can be achieved to any degree of approximation, the benefit of this method is that it provides us with a direct way to resolve issues, in other words, without the need for disturbances, planning, enormous calculations and any other transfers . The researchers Abdul M. Wazwaz and Salah M. EL-sayed in  conducted a comparative study between new modification and the modified decomposition method by introducing a new modification of the adomian decomposition method. Wazwaz, Randolph R. and Jun-Sheng D. in  used Adomian decomposition method for solving the systems of the Volterra integral forms of the Lane-Emden equations. M. Almazmumy, F. A. Hendi, H. O. Bakodah, and H. Alzumi in  presented some modifications to the Adomian decomposition Method in Ordinary Differential Equations to solve problems of Initial Value. Lina S. and weiguo W. in  used improved Adomian decomposition method is employed to solve analytic approximate solutions of the coupled fractional Burgers equations and the single nonlinear fractional reaction-diffusion equation with nonlinear terms of any order.
In this paper, we derive a new formula for solving non-linear ordinary and partial differential equations, using bernstein’s polynomial with the adomian decomposition formula. In section 2, the basic ideas of Adomian’s polynomial and modified Bernstein’s polynomial are described. Section 3, a new formula for solving ordinary and partial nonlinear differential equations will be derived based on the combination of Bernstein’s polynomial and Adomian decomposition formulas. The results and comparisons of the numerical solutions are presented in section 4, and concluding remarks are given in section 5.
2. Adomian and Bernstein Polynomials
The non-linear function N(u) can be expressed by means of a power series whose radio of convergence is infinite, that is 
, . (1)
Assuming the above hypotheses, the series whose terms are the Adomian polynomials results to be a generalization of the Taylor’s series
Is worthy to note that (2) is a rearranged expression of the series (1), and note that, due to hypothesis, this series is convergent. Consider now, the parametrization proposed by G. Adomian in  given by
where is a parameter in R and f is a complex-valued function such that . With this choosing of f and using the hypotheses above stated, the series (3) is absolutely convergent.
Substituting (3) in (2), we obtain
is the absolute convergence. We can rearrange in order to obtain the series of the form . Using (4) we can obtain the coefficients de , and finally, we deduce the Adomian’s polynomials. That is,
Using Equation (7) making and taking derivative at both sides of the equation, we can make the following identification
Hence we have obtained equation:
The Bernstein polynomials of degree n in the interval are defined by
There are n + 1 nth-degree Bernstein polynomials. For mathematical convenience, we usually set , if or .
The Bernstein approximation to a function is the polynomial .
If the 2mth order derivative is bounded in the interval then for each we have modified Bernstein polynomial .
where is the jth central moment of random variable with a binomial distribution with parameters n and x. Clearly, , . It is well known that the sequence satisfies the following recurrence
3. Modified Bernstein-Adomian Polynomials Method (MBAPM)
In this section, a new formula for solving ordinary and partial nonlinear differential equations will be derived based on the combination of Bernstein’s polynomial and Adomian decomposition formulas.
Consider the nonlinear differential equation
where L is the highest-order derivative which is assumed to be invertible, R is a linear differential operator of less order than L, N is the nonlinear operator and g is the source term. If we apply the operator which is the inverse of the L to the Equation (15), we get
Let us suppose the solution of the Equation (15):
Now we will treat the nonlinear term in Equation (17) using the modified Bernstein polynomial developed in the period as:
Applying Taylor’s formula to the term around the point
By converting Taylor’s boundary into an iterative formula based on the Adomian polynomial, we get:
Using the Equations (20)-(21), we get:
where f is calculated from the source term and the given condition(s) which are assumed to be prescribed. We now construct the recursive relation as:
It can be easily said that the solution is
In this section, we apply our new formula to solve examples of nonlinear ordinary and partial differential equations; we adopt the following four examples.
Consider the ﬁrst order nonlinear ordinary differential equation of the form  :
subject to the initial condition
Hence, the exact solution of Equation (22) is given by
Now, using Equation (21), we have n = 3, in interval [0, 1], that is a = 0 and b = 1.
The absolute and mean square errors are presented in Table 1 and Figure 1.
Consider the first order initial value problem with nonlinearity  :
The exact solution of this problem can be expressed as
Figure 1. Illustrates the absolute error between the numerical solution and the exact solution when n = 5.
Table 1. Comparison between MBAPM and the exact solution for example 1 when m = 2.
Now, using Equation (21), we have n = 3, in interval [0, 1]
Consider the following hyperbolic nonlinear problem  :
with the initial condition
Equation (29) has the exact solution ,.
Now, using Equation (21), we have n = 3, [0, 1]
Let us consider the Problem
, , (30)
with the initial condition
Equation (19) has the exact solution  :
Now, using Equation (21), we have n = 3, [0, 1]
We note from Table 3, Table 4 and Figure 3, Figure 4 that the mean square errors are 10−11 and increases when the value n increases until it reaches the exact solution when.
Figure 2. Illustrates the absolute error between the numerical solution and the exact solution when n = 5.
Figure 3. Illustrates the absolute error between the numerical solution and the exact solution when n = 3.
Figure 4. Illustrates the absolute error between the numerical solution and the exact solution when n = 5.
Table 2. Comparison between MBAPM and the exact solution for example 2 when m = 2.
Table 3. Comparison between MBAPM and the exact solution for example 3 when m = 2.
Table 4. Comparison between MBAPM and the exact solution for example 4 when m = 2.
In this study, a new formula based on modified Bernstein’s polynomial has been derived to solve the ordinary and partial nonlinear differential equations by adding the Adomian decomposition formula to Bernstein’s terms. The main benefit of this method is to convert the solution using Bernstein’s polynomial from the method of converting a nonlinear differential equation to a nonlinear system of equations whose solution requires the use of other numerical methods such as fixed point, Newton-Raphson or other methods to iterative method that finds the solution directly as shown in the examples (1 - 4). We use Bernstein expansions with Adomian decomposition method of the source term to obtain more accurate results. Figures 1-4 enable us to see that the difference between the numerical solutions and exact solutions by graphically. Tables are also given to show the variation of the absolute errors for larger approximation, namely for larger n. Maple 15 is used for calculations and sketching graphs.
The research is supported by College of Computer Sciences and Mathematics, University of Mosul, Republic of Iraq.
 Yousefi, S.A., Barikbin, Z. and Dehghan, M. (2012) Ritz-Galerkin Method with Bernstein Polynomial Basis for Finding the Product Solution Form of Heat Equation with Non-Classic Boundary Conditions. International Journal of Numerical Methods for Heat, 22, 39-48.
 Maleknejad, K., Hashemizadeh, E. and Basirat, B. (2012) Computational Method Based on Bernstein Operational Matrices for Nonlinear Volterra-Fredholm-Ham- merstein Integral Equations. Communications in Nonlinear Science and Numerical Simulation, 17, 52-61.
 Ordokhani, Y. and Far, S.D. (2011) Application of the Bernstein Polynomials for Solving the Nonlinear Fredholm Integro-Differential Equations. Journal of Applied Mathematics and Bioinformatics, 1, 13-31.
 AL-Shbool, M. and Hashim, I. (2015) Bernstein Polynomials for Solving Nonlinear Stiff System of Ordinary Differential Equations. AIP Conference Proceedings, 1678, Article ID: 060015.
 Qasim, A.F. and AL-Rawi, E.S. (2018) Adomian Decomposition Method with Modified Bernstein Polynomials for Solving Ordinary and Partial Differential Equations. Journal of Applied Mathematics, 2018, Article ID: 1803107. https://doi.org/10.1155/2018/1803107
 Wazwaz, A.M. and El-Sayed, S.M. (2001) A New Modification of the Adomian Decomposition Method for Linear and Nonlinear Operators. Applied Mathematics and Computation, 122, 393-405. https://doi.org/10.1016/S0096-3003(00)00060-6
 Wazwaz, A.M., Rach, R. and Duan, J.S. (2014) A Study on the Systems of the Volterra Integral Forms of the Lane-Emden Equations by the Adomian Decomposition Method. Mathematical Methods in the Applied Sciences, 37, 10-19. https://doi.org/10.1002/mma.2776
 Almazmumy, M., Hendi, F.A., Bakodah, H.O. and Alzumi, H. (2012) Recent Modifications of Adomian Decomposition Method for Initial Value Problem in Ordinary Differential Equations. American Journal of Computational Mathematics, 20, 228.
 Song, L. and Wang, W. (2013) A New Improved Adomian Decomposition Method and Its Application to Fractional Differential Equations. Applied Mathematical Modelling, 37, 1590-1598. https://doi.org/10.1016/j.apm.2012.03.016
 González-Gaxiola, O. and Bernal-Jaquez, R. (2017) Applying Adomian Decomposition Method to Solve Burgess Equation with a Non-Linear Source. International Journal of Applied and Computational Mathematics, 3, 213-224. https://doi.org/10.1007/s40819-015-0100-4
 Bellman, R., Kashef, B.G. and Casti, J. (1972) Differential Quadrature: A Technique for the Rapid Solution of Nonlinear Partial Differential Equations. Journal of Computational Physics, 10, 40-52. https://doi.org/10.1016/0021-9991(72)90089-7