A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

Similar documents
Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

A new Approach for Solving Linear Ordinary Differential Equations

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

A Hybrid Variational Iteration Method for Blasius Equation

Lecture 2 Solution of Nonlinear Equations ( Root Finding Problems )

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

Grid Generation around a Cylinder by Complex Potential Functions

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

: Numerical Analysis Topic 2: Solution of Nonlinear Equations Lectures 5-11:

1 GSW Iterative Techniques for y = Ax

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

Solution for singularly perturbed problems via cubic spline in tension

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

Deriving the X-Z Identity from Auxiliary Space Method

One-sided finite-difference approximations suitable for use with Richardson extrapolation

Chapter Newton s Method

Inexact Newton Methods for Inverse Eigenvalue Problems

Fixed point method and its improvement for the system of Volterra-Fredholm integral equations of the second kind

1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions

CHAPTER 7 CONSTRAINED OPTIMIZATION 2: SQP AND GRG

4DVAR, according to the name, is a four-dimensional variational method.

Finite Element Modelling of truss/cable structures

Lecture 21: Numerical methods for pricing American type derivatives

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Septic B-Spline Collocation Method for the Numerical Solution of the Modified Equal Width Wave Equation

Counterexamples to the Connectivity Conjecture of the Mixed Cells

Chapter 4: Root Finding

Solution of Singularly Perturbed Differential Difference Equations Using Higher Order Finite Differences

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

High resolution entropy stable scheme for shallow water equations

CISE301: Numerical Methods Topic 2: Solution of Nonlinear Equations

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

EEE 241: Linear Systems

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Vector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

A MODIFIED NEWTON METHOD FOR SOLVING NON-LINEAR ALGEBRAIC EQUATIONS

Appendix B. The Finite Difference Scheme

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

FUZZY TOPOLOGICAL DIGITAL SPACE OF FLAT ELECTROENCEPHALOGRAPHY DURING EPILEPTIC SEIZURES

Research Article Green s Theorem for Sign Data

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

Convexity preserving interpolation by splines of arbitrary degree

On the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros

Newton s Method for One - Dimensional Optimization - Theory

Existence results for a fourth order multipoint boundary value problem at resonance

The Study of Teaching-learning-based Optimization Algorithm

APPENDIX A Some Linear Algebra

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

6.3.4 Modified Euler s method of integration

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

New numerical scheme for solving Troesch s Problem

Bézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0

Maejo International Journal of Science and Technology

( ) [ ( k) ( k) ( x) ( ) ( ) ( ) [ ] ξ [ ] [ ] [ ] ( )( ) i ( ) ( )( ) 2! ( ) = ( ) 3 Interpolation. Polynomial Approximation.

ME 501A Seminar in Engineering Analysis Page 1

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

On a direct solver for linear least squares problems

MMA and GCMMA two methods for nonlinear optimization

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Neumann Asymptotic Eigenvalues of Sturm-liouville Problem with Three Turning Points

Relaxation Methods for Iterative Solution to Linear Systems of Equations

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

1. Statement of the problem

ELE B7 Power Systems Engineering. Power Flow- Introduction

Numerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method

Quadratic Bezier Homotopy Function for Solving System of Polynomial Equations

PART 8. Partial Differential Equations PDEs

Modified Mass Matrices and Positivity Preservation for Hyperbolic and Parabolic PDEs

System in Weibull Distribution

Binomial transforms of the modified k-fibonacci-like sequence

Least squares cubic splines without B-splines S.K. Lucas

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

A Fast Computer Aided Design Method for Filters

Heuristic Algorithm for Finding Sensitivity Analysis in Interval Solid Transportation Problems

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

SHIFTED JACOBI COLLOCATION METHOD BASED ON OPERATIONAL MATRIX FOR SOLVING THE SYSTEMS OF FREDHOLM AND VOLTERRA INTEGRAL EQUATIONS

IV. Performance Optimization

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

Expectation Maximization Mixture Models HMMs

CHAPTER 4d. ROOTS OF EQUATIONS

Nonlinear Overlapping Domain Decomposition Methods

Sampling Theory MODULE V LECTURE - 17 RATIO AND PRODUCT METHODS OF ESTIMATION

SIO 224. m(r) =(ρ(r),k s (r),µ(r))

RELIABILITY ASSESSMENT

Summary with Examples for Root finding Methods -Bisection -Newton Raphson -Secant

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Transcription:

Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS Kanttha Chompuvsed Department of Mathematcs and Appled Statstcs, Faculty of Scence and Technology, Nahon Ratchasma Rajabhat Unversty, Nahon Ratchasma, Thaland Receved 1-1-8, Revsed 1--16; Accepted 1--15 ABSTRACT Solvng systems of nonlnear equaton s a great mportant whch arses n varous branches of scence and engneerng. In the last decades, several numercal technques were proposed to solve these problems. In ths study, we propose a modfed of teratve method whch s based on the dea of and Fed pont teraton method. The proposed method has been llustrated wth several eamples from the reference. The numercal results ndcate that ths proposed method provde the good performance of teratons. Keywords: System of Nonlnear Equaton, Newton Method, Fed Pont Iteraton Scence Publcatons 1. INTRODUCTION Solvng systems of nonlnear equatons s a great mportance, because these systems frequently arse n varous branches of pure and appled scences. The general form of a system of nonlnear equatons s Equaton 1: f 1( 1,,, n) =,f ( 1,,, n) =,f n( 1,,, n) = (1) where, each functon f can be thought of as mappng a vector = ( 1,,, n ) of the n dmensonal space R n, nto the real lne R. The system can alternatvely be represented by defnng a functonal F, mappng R n nto R n by. F( 1,,, n ) = (f 1 ( 1,,, n ),,f n ( 1,,, n )) T Usng vector notaton to represent the varables 1,,, n, a system (1) can be wrtten as the form: F() = The functons f 1, f,,f n are called the coordnate functons of F (Burden and Farres, 1). Recently, several teratve methods have been used to solve nonlnear equatons and the system of nonlnear equatons (Awawdeh, 9; Noor, 1; Cordero et al., 11; Sharma and Sharma, 11; Vahd et al., 1). 4 Wang (11) usng a thrd order famly of Newton-Le teraton method for solvng nonlnear equatons; Ozel (1) has consdered a new decomposton method for solvng the system of nonlnear equatons. Saha (1) has presented a modfed method to solvng nonlnear equatons by hybrdsng the results of and fed pont teraton method. Km et al. (1) developed a new scheme for the constructon of teratve methods for the soluton of nonlnear equatons and gvng a new class of methods from any teratve method. Furthermore, several teratve methods have been developed for solvng the system of nonlnear equatons by usng varous technques such as Newton s method, Revsed Adoman decomposton method, homotopy perturbaton method, Householder teratve method (Darvsh, 9; Noor and Waseem, 9; Hossen and Kafash, 1; Darvsh and Shn, 11; Hafz and Bahgat, 1a; 1b; Noor et al., 1). It s the purpose of ths study to ntroduce a new mprovement of by fed pont teraton method. We etend the Saha (1) method to solve systems of nonlnear equatons. Some eamples are tested and the obtaned results suggest that ths newly mprovement technque ntroduces a promsng tool and powerful mprovement for solvng a System of Nonlnear Equatons.

Kanttha Chompuvsed /Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 1.1. Descrpton of an Iteratve Method Consder a nonlnear equaton: f() = () We assume that the Equaton () admts a unque soluton *. In the all nown teratve formula used to fnd the real root s Equaton : +1 = - f ( ) Scence Publcatons f( ) () In fed pont teraton method () wll be rewrtten n the form: = g() (4) Equaton 4 whch s equvalent to () wll converge to a real root n the nterval D f g`() <1 for all n D provded the ntal appromaton s chosen n D. Choose the ntal appromaton then (, g ( )) s a pont on the curve Equaton 5: y = g() (5) The equaton of the tangent to the curve gven by (5) at the pont (, g ( )) s Equaton 6: y - g( ) = g ( )( - ) (6) Now we consder the lne Equaton 7: y = (7) Substtutng y = n (6) we have: - g( ) = g ( )( - ) g( ) - g ( ) [1- g ( )] = g( ) - g ( ) = 1- g ( ) whch produces the followng teraton scheme Equatons 8: g( ) - g ( ) +1 = 1- g ( ) 1.. The N-Dmensonal Case (8) The (Gautsch, 11; Sauer, 11) s commonly used for solvng such systems Equaton 9: 5 F() = (9) where, F: Ω R n R n s defned Equaton 1: +1 = - F ( ) F( ) (1) where, F ( ) s the Jacoban matr n pont. In fed pont teraton method (9) wll be rewrtten n the form = g() We rewrte Equaton 8 to solve the nonlnear system F() =, ths produces the followng teraton scheme Equatons 11: = [I - g ( )] [g( ) - g ( ) ] (11) -1 +1 where, I s an dentty matr. 1.. Numercal Eamples We present some eamples to llustrate the effcency of our proposed methods, we solve four systems of nonlnear equatons and one of a nonlnear boundary value problem. The followng tables show the Number of Iteratons (NI) to receve the requred soluton. For all test problems the stop crtera s F() <1 9. Eample 1 Consder the followng system of nonlnear equatons: + - = - = 1 1 * * * T T The eact solutons are = ( 1, ) = (1,1). To solve ths system, we set = (.1,) T as an ntal value. The results are presented n Table 1. Eample Consder the followng system of nonlnear equatons (Hossen and Kafash, 1): + - 6 + = - - 6 - = 1 1 1 The eact solutons are * * * T T 1 = (, ) = (1,1) = (.57779,.51574475988) T To solve ths system, we set = (.5,.5) as an ntal value. The results are presented n Table.

Kanttha Chompuvsed /Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 Table 1. Numercal results for Eample 1 ------------------------------------------------------ ----------------------------------------------------- NI X 1 X X 1 X 1 1.5. 1.414156.818471 5.84997 1. 1.767476 1.8941756 5.1465777 1.1176476 1.176444.89185 4 1.571718 1.45777.87495814.814856 5 6.49996 1.1 1.9689551 1.4786 6.17718 1..999759947.999744 7 1.7788 1. 1. 1.5 8 1.16758157 1. 1. 1. 9 1.16468 1. 1 1.71459 1. 11 1. 1. 1 1. 1. Table. Numercal results for Eample ------------------------------------------------------- ---------------------------------------------------------- NI X 1 X X 1 X 1.48717845 -.8914615.48717845 -.8914615.5187686 -.568678.5187686 -.568678.5184855 -.55754.5184855 -.55754 4.518485 -.557478.518485 -.557478 Table. Numercal results for Eample -------------------------------------------------------------------- ------------------------------------------------------------------------- NI X 1 X X X 1 X X 1 4.784155 1.9877697 -.47598777 5.6516177.55774445 -.47598777.541.89418971 -.476619.5119999 -.1717 -.488991715.4975994.4114 -.515965.5588.1168 -.55754 4.511778.1617795 -.519818.5 -. -.5598776 5.54198.497788 -.561 6.57456.86841 -.58611 7.588.1699 -.5595 8.55.5 -.559876 9.5. -.5598776 Table 4. Numercal results for Eample 4 Number of teratons ------------------------------------------------ Method m = 5 m = 75 m = 1 6 6 6 4 4 4 Table 5. Numercal results for Eample 5 Number of teratons ------------------------------------------------ Method M = 5 M = 75 M = 1 6 7 7 6 7 7 Scence Publcatons 6 Eample Consder the followng system of nonlnear equatons (Awawdeh, 9): 1 - cos( ) -.5 = 1-81( +.1) + sn +1.6 = -1 1π - e + + = The eact solutons = (,, ) = (.5,,-.5598775598). To solve are * * * * T T 1

Kanttha Chompuvsed /Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ths system, we set = (5, 4, ) T as an ntal value. The results are presented n Table. Eample 4 Consder the followng system of nonlnear equatons (Darvsh and Shn, 11): Scence Publcatons - cos( -1) =, = 1,,,m The eact solutons = (,,, ) = (1,1,,1). To solve ths system, are * * * * T T 1 m we set = (.5,.5,,.5) T as an ntal value. The results are presented n Table 4. Eample 5 Consder the nonlnear boundary value problem (Noor and Waseem, 9): y = -(y) - y + ln, 1, y(1) =, y() = ln Whose eact solutons s y = In. We consder the followng partton of the nterval: 1 = 1, n =, j = + jh, h =, j = 1,,,m -1 m Let us defne now: y = y( ) =, y m = ln, y = f( ), = 1,,,m -1 If we dscretze the problem by usng the second order fnte dfferences method defned by the numercal formulas: y - y y = h y+1 - y + y-1 =, = 1,,,m -1 h +1-1, = 1,,,m -1, y Then, we obtan a (m-1) (m-1) system of nonlnear equatons: 4y + y + 4y (h - ) - 4h ln =,4(y + y ) 1 1 +1-1 +(y - y ) + 4y (h - ) - 4h ln =, =,,m -, +1-1 ( ) ( ) ( ) 4 In + y + In y + 4y h 4h In = m m m 1 m 1 we tae X wth () y = ln, = 1,,,m -1, as a 1 startng pont. The results are presented n Table 5. 7.CONCLUSION In ths study, we have demonstrated the applcablty of the modfed method for the system of nonlnear equatons wth the help of some concrete eamples. The results show that: the proposed problem can be solved by the proposed method.. REFERENCES Awawdeh, F., 9. On new teratve method for solvng systems of nonlnear equatons. Numercal Algorthms, 54: 95-49. DOI: 1.17/s1175-9-94-8 Burden, R.L. and J.D. Farres, 1. Numercal Analyss. 9th Edn., Cengage Learnng, Boston, MA., ISBN-1: 587519, pp: 87. Cordero, A., J.L. Hueso, E. Martnez, J.R. Torregrosa, 11. Effcent hgh-order methods based on golden rato for nonlnear systems. Appled Math. Comput., 17: 4548-4556. DOI: 1.116/j.amc.1.11.6 Darvsh, M.T. and B.C. Shn, 11. Hgh-order newtonrylov methods to solve systems of nonlnear equatons. J. KSIAM., 15: 19-. Darvsh, M.T., 9. A two-step hgh order newton-le method for solvng systems of nonlnear equatons. Int. J. Pure Appled Math., 57: 54-555. Gautsch, W., 11. Numercal Analyss. nd Edn., Sprnger, Boston, ISBN-1: 81768597, pp: 588. Hafz, M.A. and M.S.M. Bahgat, 1a. An effcent two-step teratve method for solvng system of nonlnear equatons. J. Math. Res., 4: 8-4. DOI: 1.559/jmr.v4n4p8 Hafz, M.A. and M.S.M. Bahgat, 1b. Modfed of householder teratve method for solvng nonlnear systems. J. Math. Comput. Sc. Hossen, M.M. and B. Kafash, 1. An effcent algorthm for solvng system of nonlnear equatons. Appled Math. Sc., 4: 119-11. Km, Y.L., C. Chun and W. Km., 1. Some thrdorder curvature based methods for solvng nonlnear equatons. Stud. Nonlnear Sc., 1: 7-76. Noor, M.A. and M. Waseem, 9. Some teratve methods for solvng a system of nonlnear equatons. Comput. Math. Appl., 57: 11-16. DOI: 1.116/j.camwa.8.1.67 Noor, M.A., 1. Iteratve methods for nonlnear equatons usng homotopy perturbaton technque. Appled Math. Inform. Sc., 4: 7-5.

Kanttha Chompuvsed /Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 Noor, M.A., M. Waseem, K.I. Noor and E. Al-Sad, 1. Varatonal teraton technque for solvng a system of nonlnear equatons. Optm Lett., DOI: 1.17/s1159-1-479- Ozel, M., 1. A new decomposton method for solvng system of nonlnear equatons. J. Appled Math. Comput., 15: 89-95. Saha, S., 1. A modfed method for solvng nonlnear equatons. Int. J. Comput. Sc. Intell. Comput., : 6-11. Sauer, T., 11. Numercal Analyss. nd Edn., Prentce Hall, US., ISBN-1: 1818768, pp: 646. Sharma, J.R. and R. Sharma, 11. Some thrd order methods for solvng systems of nonlnear equatons. World Acad. Sc. Eng. Technol., 6: 194-11. Vahd, A.R., S.H. Javad and S.M. Khorasan, 1. Solvng system of nonlnear equatons by restarted adoman s method. Appled Math. Comput., 6: 59-516. Wang, P., 11. A thrd-order famly of newton-le teraton methods for solvng nonlnear equatons. J. Num. Math. Stochast., : 1-19. Scence Publcatons 8