Research Article  Open Access
MoinudDin Junjua, Saima Akram, Nusrat Yasmin, Fiza Zafar, "A New JarrattType FourthOrder Method for Solving System of Nonlinear Equations and Applications", Journal of Applied Mathematics, vol. 2015, Article ID 805278, 14 pages, 2015. https://doi.org/10.1155/2015/805278
A New JarrattType FourthOrder Method for Solving System of Nonlinear Equations and Applications
Abstract
Solving systems of nonlinear equations plays a major role in engineering problems. We present a new family of optimal fourthorder Jarratttype methods for solving nonlinear equations and extend these methods to solve system of nonlinear equations. Convergence analysis is given for both cases to show that the order of the new methods is four. Cost of computations, numerical tests, and basins of attraction are presented which illustrate the new methods as better alternates to previous methods. We also give an application of the proposed methods to wellknown Burger's equation.
1. Introduction
Let us consider the problem of constructing an optimal iterative method to find a simple zero of a nonlinear equationas well as its extension for solving nonlinear system of equations . There is no ambiguity that the quadratically convergent Newton method (NM) is one of the best root finding methods based on two evaluations of function for approximating the solution of a nonlinear equation (1) and is given asThe natural extension of Newton’s method for system of nonlinear equations is given [1] as A large number of variants of the classical Newton method to obtain higherorder multistep schemes with better efficiencies have appeared using numerous techniques [2–4]. However, the research remained focussed on obtaining optimal and computationally efficient methods. The main objective was not only to boost up the order and speed of convergence but also to reduce the computational cost. This objective was achieved by defining twostep methods. A very few of them had less computational cost. Weerakoon and Fernando [5] proposed an accelerated cubically convergent twostep variant of Newton’s method:where . Another cubically convergent without memory method was suggested by Frontini and Sormani [2]:However, the above methods were not optimal. Jarratt [6] constructed an optimal method given byand its extension for the systems of nonlinear equations is given by [7]Several researchers have also focussed their attention on obtaining Jarratttype iterative schemes for solving nonlinear equation [8–11]. Khattri and Abbasbandy [10] gave their contribution by using concept of parameter approach to develop an optimal fourthorder Jarratttype iterative scheme requiring one evaluation of the function and two evaluations of first derivative. The following is a special case of their scheme:Soleymani et al. [11] replaced the parameter approach by the weight function approach to construct a family of optimal Jarratttype twostep methods:where and . A special case of their scheme is as follows:Babajee et al. [12] extended (10) for the multivariate case as follows:where, and is the identity matrix. Chun et al. [8] used Halley’s method with an approximation of second derivative by using weight function approach in the second step of their scheme to obtain optimal variant of Jarratttype fourthorder method. The classical Jarratt family of fourthorder methods was obtained as special case. A special case of his scheme is given asIn this contribution, we present and analyze an optimal family of fourthorder convergent iterative schemes using weight function in the second step. We, then, extend it for multivariate case. The rest of the paper is organized as follows. Section 2 is comprised of the construction and convergence analysis of the new family of methods for the single nonlinear equations. Section 3 consists of the extension of our method to the system of nonlinear equations. Section 4 includes numerical tests. Sections 5 and 6 provide the cost of computations and basins of attraction, respectively, for the sake of comparison of the new methods with the existing ones in this domain. In Section 7 an application of the proposed method to wellknown Burger’s equation is presented and concluding remarks are given in Section 8.
2. New Family of JarrattType Methods and Its Convergence Analysis
In this section, we propose a new optimal fourthorder Jarratttype scheme for computing the zeros of a univariate nonlinear function. We give a twostep scheme in which the first step is similar to Jarratt’s scheme. In the second step, we use weight functions approach as follows:where , , and and represent realvalued weight functions chosen such that new scheme (14) achieves optimal fourthorder convergence as stated in Theorem 1.
Theorem 1. Let be a simple zero of sufficiently differentiable function in an open interval containing . Then, for , the new withoutmemory scheme (14) has optimal convergence of order four under the following conditions on weight functions:and it satisfies the error equation given bywhere , .
Proof. Let be the error at th computing step. By using Taylor’s expansion of about the root , we getwhere , . The first derivative in the first step of our scheme can be calculated asUsing (17) and (18) in the first step of (14), the Taylor expansions for and are given byIn the similar manner, for the second step of (14), we haveBy using (17), (18), and (20), we getAgain using Taylor’s expansion, we haveSubstituting , (21) and (23) in (14) and using (15), we have which gives the following error equation:Hence, it can be seen that the new scheme has optimal fourthorder convergence.
3. Extension of the New Family for Multivariate Case and Its Analysis
In this section, we extend our new family of optimal Jarratttype schemes to solve systems of nonlinear equations. We give a special case for our new family by defining weight functions satisfying the conditions of Theorem 1 as follows:Thus, we achieve a new Jarratttype fourthorder method:Now, let us define to be a sufficiently Fréchet differentiable in a convex set , where is an open convex neighborhood of the root , is th approximate root of the exact root , and ; we can write the th derivative of at , as a linear function such that as described in [13]. Thus, it is defined asAlso, where are the Hessian matrices of . Therefore, Taylor’s series for function of variables can be written aswhere , is continuous and nonsingular, and is closer to . Now, we extend our scheme for solving a system of nonlinear equations as follows:where , , is the identity matrix, and . By using the above Taylor expansion we can prove the following theorem.
Theorem 2. Let be an open convex set containing the root of and let be fourtime Fréchet differentiable in such that the Jacobian matrix is continuous and nonsingular in . Then the new method (31) has convergence of order four.
Proof. Let be the solution of nonlinear system , and let be an initial guess close to ; then by Taylor’s expansion of about , we havewhere , , and . We can calculateto getThe Jacobian matrix has the following Taylor expansion:to calculateUsing (33) and (36) in second step of (31), we attain with the following error term:which shows the proposed method (31) is fourthorder convergent.
4. Numerical Results
We, now, check the effectiveness of our new optimal fourthorder family of methods (27) (MSNF1) by comparing it with Newton’s method (2) (NM), Khattri’s method (8) (KM), Soleymani’s method (10) (SM), and Chun’s method (13) (CM). Algorithms have been executed in Maple software and are tested for the examples given in Table 1. The approximate solutions for solving single equations are calculated by using a precision of decimal digits with the stopping criterion . In Tables 2, 3, 4, 5, and 6, number of iterations, “,” absolute values of functions and absolute values of the difference between approximated roots and exact roots, , for each iterative method are given. We also compare our proposed method (31) (MSNF2) for solving systems of nonlinear equations with Newton’s method (3) (NM), Babajee’s method (11) (BM), and Khattri’s method (8) (KM). For solving systems of equations results are computed with a precision of decimal digits. The stopping criterion is . Table 7 consists of test functions for solving systems of nonlinear equations along with their exact zeros. Table 8 shows absolute values of the difference between two consecutive approximations of the root and absolute functional values satisfying the above stopping criterion for each of the methods.




