International Journal of Mathematics and Mathematical Sciences

Volume 2019, Article ID 2872867, 8 pages

https://doi.org/10.1155/2019/2872867

## Treating Transcendental Functions in Partial Differential Equations Using the Variational Iteration Method with Bernstein Polynomials

College of Computer Sciences and Mathematics, University of Mosul, Iraq

Correspondence should be addressed to Ahmed Farooq Qasim; moc.oohay@laciremundemha

Received 25 November 2018; Revised 3 January 2019; Accepted 11 February 2019; Published 3 March 2019

Academic Editor: Shyam L. Kalla

Copyright © 2019 Ahmed Farooq Qasim and Almutasim Abdulmuhsin Hamed. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We aim through this paper to present an improved variational iteration method (VIM) based on Bernstein polynomials (BP) approximations to be used with transcendental functions. The key benefits gained from this modification are to reach stable and fairly accurate results and, at the same time, to expand the unknown function’s domain in partial differential equations (PDEs). The proposed approach introduces the Bernstein polynomials in the transcendental functions of nonlinear PDEs. A number of examples were included in order to expound the method’s capacity and reliability. From the results, we conclude that the VIM with BP is a powerful mathematical tool that can be applied to solve nonlinear PDEs.

#### 1. Introduction

It is common knowledge that a large number of phenomena are in essence nonlinear and therefore can be modeled using nonlinear differential equations [1]. However, it is not possible to find exact solutions to nonlinear differential equations. Therefore, the approach to tackle this shortcoming was to apply different numerical methods suffice to reach approximate solutions including the methods of Adomian decomposition [2, 3], homotopy perturbation [4–6], differential transform [7–9], and variational iteration [10, 11].

He [12] introduced the VIM in 1998, and since then it was applied by many scientists as an approach to solve mathematical and physical problems for its capability to simplify a difficult problem into a readily solvable problem. The VIM is considered a capable tool that can be used to solve functional equations. Commonly used numerical approaches such as finite difference or characteristics methods require tedious computations. These methods are often less accurate in their results because of their round-off error. Solutions for Schrodinger equations based on the widely used analytical methods are very constrained and can be only applied in certain scenarios while they are inapplicable to obtain solutions for a large number of equations that model real life scenarios. The VIM is a modified Lagrange multiplier method with promising results in solving numerous classes of nonlinear problems with approximations that converge to the correct solution within a short number of iterations.

Bernstein polynomials’ role is prominent in countless mathematics applications as they are utilized in solving differential equations in addition to their contribution to the approximation theory. Bézier curves, which are Bernstein polynomials restricted to the interval , played a more prominent role with the development of computer graphics [13].

A novel iterative approach based on the VIM is applied to differential equations by [14]. Olayiwola Mo [15] applied it to solve the Convection-Diffusion equation and for the class of fractional Convection-Diffusion equation [16]. The modified variational iteration method for sine-Gordon equation is suggested in [17] using Chebyshev polynomials and He’s polynomials [18].

In this paper, we adjusted the VIM to utilize the Bernstein polynomials so as to approximate the correction functions’ nonlinear terms. The solutions obtained using our modified method are comparable to VIM. In our approach, we applied the Bernstein polynomial approximations to nonlinear functions in the correction functions before iterating the numerical solution of PDEs. The main objective is to achieve reliable accurate PDEs solutions if using VIM leads to unstable solutions. The feasibility of our approach is expounded through a number of examples.

The rest of the paper is divided into five sections. Section 2 introduces brief ideas of VIM. Section 3 presents one- and two-dimensional Bernstein polynomials. In Section 4, an analysis of the convergence is presented. Section 5 elucidates three examples that expound the feasibility of using VIM with the Bernstein polynomials. Section 6 is dedicated for the conclusions.

#### 2. VIM Basics

The concept of VIM is based on the general Lagrange’s multiplier method. The key advantage of this method is reaching a mathematical problem solution by linearization assumption to be used as an initial approximation that promptly converges to an exact solution [19].

The following nonlinear differential equation elucidates the basic concept of the VIM, as follows:where L and N are the linear and nonlinear operators, respectively, while g(x) is an inhomogeneous term. VIM is used to formalize the following correction function:where is a general Lagrangian multiplier which can be determined through the variational theory and integration by parts, represents a restricted variation , and the subscript symbolizes the nth-order approximation.

The Lagrange multiplier is optimally determined through integration by parts. The consecutive approximations , of the solution are determined using the Lagrange multiplier from the first step using any selective function . Thus, the solution will be as follows.

#### 3. Bernstein Polynomials

Polynomials are considered easily defined, calculated, differentiated, and integrated mathematical tools. The Bernstein based polynomials are tested to approximate the functions as they result in improved approximations to a function with a few terms. This feature has led to their widespread use in applied mathematics, physics, and computer based geometric designs, combined with other methods like Galerkin and collocation to solve both differential and integral equations.

The nth degree Bernstein Polynomials form a complete basis over , formalized asin which represent the binomial coefficients.

Based on [13], we use the generalized Bernstein polynomials of nth degree which form a complete basis over as follows.

If the function of two real variables over the interval square are

then the two-variable Bernstein polynomial of degree, equivalent to the function , is determined in the following formula [20]:where

We use (7) to convert the transcendental functions (trigonometric, exponential, and Logarithm functions) that appear in homogenous and nonhomogenous PDEs, such as the sine-Gordon equation, into the series approximation by means of Bernstein polynomial.

The following nonlinear partial differential equation contains the transcendental functions term as follows:such that is the transcendental functions.

Using (2), we have the following.Therefore, we approach the transcendent functions using Bernstein’s approximation, because it can only be integrated by numerical methods.

#### 4. Analysis of the Convergence

The variational iteration method converts the partial differential equations to a recurrence sequence of functions. The limit of that sequence leads to the solution of the PDE.

Now, letbe the second- order partial differential equation and be the transcendent functions. Using VIM with Bernstein polynomial in (12), we havesuch that is He’s monographs, i.e., [21]. To find the value , (13) becomes as follows.

In fact the solution of (14) is considered as the fixed point with the initial condition [21].

Theorem 1 (see [22]). *Let a sequence satisfy and as . Then for any function ,In other words, is Bernstein polynomials and uniformly convergent to .*

Theorem 2 (see [21] (Banach’s fixed point theorem)). *Assume that is a Banach space andis a nonlinear mapping, and suppose thatfor some constant . Then A has a unique fixed point. Furthermore, the sequencewith an arbitrary choice of , converges to the fixed point of A andand, according to Theorem 2, for the nonlinear mappinga sufficient condition for convergence of the variational iteration method is strictly contraction of A.*

Corollary 3. *Using Theorems 1 and 2 we conclude thatsuch that is Bernstein approximation; then the sufficient condition for convergence of the variational iteration method is strictly contraction of A. Furthermore, the sequence (13) converges to the fixed point of A which is also the solution of the partial differential equation (9).*

#### 5. Numerical Examples

In this section, the following three examples of nonlinear PDEs are solved using the VIM with Bernstein polynomials. The results are generated using Maple 13. The accuracy of the results is determined when compared to the exact solutions.

*Example 1. *Let us take the sine-Gordon equation [17].the exact solution of (9) is as follows.Applying the Bernstein polynomials approximation for in (7) with n=m=2 results in the following.Applying the VIM on (22), we getwhere the Lagrange multiplier is as follows.Consider the first iterate numerical solution.The absolute error of VIM-Bernstein and the exact solution are presented in Tables 1-2 and Figures 1-2.