Cubic Bezier Homotopy Function for Solving Exponential Equations

Similar documents
Approximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy

Fuzzy Set Theory in Modeling Uncertainty Data. via Interpolation Rational Bezier Surface Function

On One Analytic Method of. Constructing Program Controls

M. Y. Adamu Mathematical Sciences Programme, AbubakarTafawaBalewa University, Bauchi, Nigeria

Chapter Lagrangian Interpolation

On computing differential transform of nonlinear non-autonomous functions and its applications

Variants of Pegasos. December 11, 2009

[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5

A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION

HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD

GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim

Performance Analysis for a Network having Standby Redundant Unit with Waiting in Repair

Solution in semi infinite diffusion couples (error function analysis)

Existence and Uniqueness Results for Random Impulsive Integro-Differential Equation

( ) () we define the interaction representation by the unitary transformation () = ()

Normal Type-2 Fuzzy Rational B-spline Curve

CH.3. COMPATIBILITY EQUATIONS. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4

Numerical solution of second-order hyperbolic telegraph equation via new cubic trigonometric B-splines approach

Robust and Accurate Cancer Classification with Gene Expression Profiling

COMPUTER SCIENCE 349A SAMPLE EXAM QUESTIONS WITH SOLUTIONS PARTS 1, 2

V.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS

Lecture VI Regression

Method of upper lower solutions for nonlinear system of fractional differential equations and applications

Comparison of Differences between Power Means 1

Ordinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s

New M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)

Volatility Interpolation

Lecture 18: The Laplace Transform (See Sections and 14.7 in Boas)

FTCS Solution to the Heat Equation

J i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.

Lecture 6: Learning for Control (Generalised Linear Regression)

Including the ordinary differential of distance with time as velocity makes a system of ordinary differential equations.

Research Article Numerical Approximation of Higher-Order Solutions of the Quadratic Nonlinear Stochastic Oscillatory Equation Using WHEP Technique

Notes on the stability of dynamic systems and the use of Eigen Values.

Comb Filters. Comb Filters

Dual Approximate Dynamic Programming for Large Scale Hydro Valleys

Testing a new idea to solve the P = NP problem with mathematical induction

by Lauren DeDieu Advisor: George Chen

Advanced Machine Learning & Perception

Mohammad H. Al-Towaiq a & Hasan K. Al-Bzoor a a Department of Mathematics and Statistics, Jordan University of

THE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS

ISSN MIT Publications

This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.

Mechanics Physics 151

TSS = SST + SSE An orthogonal partition of the total SS

Generalized double sinh-gordon equation: Symmetry reductions, exact solutions and conservation laws

Numerical Solution of Quenching Problems Using Mesh-Dependent Variable Temporal Steps

John Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany

CHAPTER 10: LINEAR DISCRIMINATION

Image Morphing Based on Morphological Interpolation Combined with Linear Filtering

Appendix H: Rarefaction and extrapolation of Hill numbers for incidence data

Relative controllability of nonlinear systems with delays in control

Li An-Ping. Beijing , P.R.China

In the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!

Implementation of Quantized State Systems in MATLAB/Simulink

Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system Y X (2) and write EOM (1) as two first-order Eqs.

Observer Design for Nonlinear Systems using Linear Approximations

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Institute for Mathematical Research, Universiti Putra Malaysia, UPM Serdang, Selangor, Malaysia

T q (heat generation) Figure 0-1: Slab heated with constant source 2 = q k

Pendulum Dynamics. = Ft tangential direction (2) radial direction (1)

Tight results for Next Fit and Worst Fit with resource augmentation

Online Supplement for Dynamic Multi-Technology. Production-Inventory Problem with Emissions Trading

ABSTRACT KEYWORDS. Bonus-malus systems, frequency component, severity component. 1. INTRODUCTION

APPROXIMATE ANALYTIC SOLUTIONS OF A NONLINEAR ELASTIC WAVE EQUATIONS WITH THE ANHARMONIC CORRECTION

MANY real-world applications (e.g. production

Quadratic Bezier Homotopy Function for Solving System of Polynomial Equations

. The geometric multiplicity is dim[ker( λi. number of linearly independent eigenvectors associated with this eigenvalue.

FI 3103 Quantum Physics

Lecture 11 SVM cont

Using Aggregation to Construct Periodic Policies for Routing Jobs to Parallel Servers with Deterministic Service Times

Structural Optimization Using Metamodels

A Novel Efficient Stopping Criterion for BICM-ID System

Linear Response Theory: The connection between QFT and experiments

P R = P 0. The system is shown on the next figure:

FRACTIONAL OPTICAL SOLITARY WAVE SOLUTIONS OF THE HIGHER-ORDER NONLINEAR SCHRÖDINGER EQUATION

SOME NOISELESS CODING THEOREMS OF INACCURACY MEASURE OF ORDER α AND TYPE β

Machine Learning Linear Regression

CS286.2 Lecture 14: Quantum de Finetti Theorems II

Attribute Reduction Algorithm Based on Discernibility Matrix with Algebraic Method GAO Jing1,a, Ma Hui1, Han Zhidong2,b

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two

The Analysis of the Thickness-predictive Model Based on the SVM Xiu-ming Zhao1,a,Yan Wang2,band Zhimin Bi3,c

Outline. Probabilistic Model Learning. Probabilistic Model Learning. Probabilistic Model for Time-series Data: Hidden Markov Model

The Dynamic Programming Models for Inventory Control System with Time-varying Demand

Solving the multi-period fixed cost transportation problem using LINGO solver

Bayesian Inference of the GARCH model with Rational Errors

Neural Networks-Based Time Series Prediction Using Long and Short Term Dependence in the Learning Process

PHYS 1443 Section 001 Lecture #4

How about the more general "linear" scalar functions of scalars (i.e., a 1st degree polynomial of the following form with a constant term )?

Unconstrained Gibbs free energy minimization for phase equilibrium calculations in non-reactive systems using an improved Cuckoo Search algorithm

Sampling Procedure of the Sum of two Binary Markov Process Realizations

EEL 6266 Power System Operation and Control. Chapter 5 Unit Commitment

The Performance of Optimum Response Surface Methodology Based on MM-Estimator

A Cell Decomposition Approach to Online Evasive Path Planning and the Video Game Ms. Pac-Man

Reconstruction of Variational Iterative Method for Solving Fifth Order Caudrey-Dodd-Gibbon (CDG) Equation

Survival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System

Econ107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)

Review of Numerical Schemes for Two Point Second Order Non-Linear Boundary Value Problems

Response of MDOF systems

Coupled Method for Solving Time-Fractional Navier-Stokes Equation

Transcription:

Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M. A. Salm Nasr,c Faculy of Compuer and Mahemacal Scences, Unvers Teknolog MARA Malaysa, 445 Shah Alam, Selangor. School of Mahemacal Scences, Unvers Sans Malaysa, 8 Mnden, Pulau Pnang, Malaysa. * s.syuhadahraml@ymal.com, a hafz97my@yahoo.com, b syazanasaharzan@yahoo.com, c masn@msk.um.edu.my Absrac In hs sudy, he new homoopy funcon has been eended from Quadrac Bezer homoopy funcon for mplemenaon wh sngle eponenal equaons. The eended funcon s a combnaon of he Bezer curve and homoopy funcon ha was based on De Caseljau algorhm. The funcon wll nclude he arge funcon and he aulary funcon. The arge funcon s he eponenal equaon whle he aulary funcon seleced mus be conrollable and easy o be solved. Then, he eended funcon has solved by usng Newon omoopy Connuaon Mehod o compue he appromae soluons, accuracy of he appromae soluons and mamum absolue error of he soluons. The analysed resuls proved ha hese omoopy funcon show beer appromae soluons, accuracy of he appromae soluons and mamum absolue error of he soluons compared o Quadrac Bezer homoopy funcon. Copyrgh 6 Penerb Akadema Baru - All rghs reserved. Keywords: Bezer curve, homoopy funcon, newon homoopy connuaon mehod, sngle of eponenal equaon. INTRODUCTION Nowadays, a research abou concep of omoopy funcon has been an neresng and wdely used. The majory of hem ulzed and solve equaon comprse of a few classfcaons, for eample, algebrac, polynomals, rgonomerc, eponenal, and logarhmc equaon []. In 9, Lahaye sar o nroduce he homoopy mehod [], bu he concep of omoopy Connuaon mehod for solvng equaon sars o be denfed and suded n 97. Prevous sudy, Wu have made a lo of sudy abou he omoopy Connuaon mehod. Wu begn wh Modfed Chnese algorhm and sudy he convergence of Newon omoopy Connuaon mehod. Then, Wu has developed a new formula and compare by usng radonal Adoman Decomposon mehod wh Adoman omoopy Connuaon mehod. Besdes ha, Wu also solve a nonlnear equaon by usng Newon omoopy Connuaon mehod and adjusng an aulary funcon. Then, Wu has modfed he secan mehod wh he omoopy Connuaon echnque o new Secan omoopy Connuaon funcon. All he resul [-6] obaned shown ha omoopy Connuaon mehod s beer han classcal mehod.

Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 In [8], Nor also have develop a new funcon of Quadrac Bezer homoopy funcon for solvng sysem of eponenal equaons. Therefore, we eend from he sudy o develop homoopy funcon for mplemenaon wh sngle eponenal equaon usng Newon omoopy Connuaon mehod o compue he appromae soluons, accuracy of appromae soluons and mamum absolue error of soluons.. CUBIC BEZIER OMOTOPY FUNCTION Frsly, nroduce he Bezer curve based on he De Caseljau algorhm as follow [9], P ( r = ( P r ( + P r + ( r =,..., k =,..., k r ( where P ( = P s he fnal pon n he curve. Then, he curve wll be consruc from he formula n ( defned by four pon P, P, P and P. Then he funcon become ( P = ( P + ( P + ( P +. ( ( P The dea of develop omoopy funcon based on De Caseljau Algorhm. In omoopy, here s a curve ha movng from one curve o anoher curve. Whle n De Caseljau, here s movemen of pon a curve and we consder here s smlary ha can relae beween hem [8]. Therefore, by usng mehod of De Caseljau, becomes easer and sysemac wh recursve consrucon of homoopy funcon s as follow, Fgure : The recursve consrucon of omoopy funcon. Based on he Fg., noe ha A ( = ( G( + ( B = ( ( + ( (

Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 C ( = ( ( + F( D ( = ( A + B E ( = ( B + C hus ( = ( D + E ( Subsue he above equaon, ( = ( [( A + B] + [( B + C] = ( A + ( B + ( B + = ( [( G( + ( ] + ( [( ( + ( ] + ( [( ( + ( ] + [( ( + F( ] = ( G( + ( ( + ( ( + ( ( + ( ( + ( ( + ( ( + C F( = ( G( + ( ( + ( ( + F( (4 = B ( + B + B + B ( = B ( (5 wh = = G( for =, = (, = ( for = (, = F( for =, and B ( s a Bernsen funcon n Bezer curve whch s defned n [] B ( = ( B ( = (( B ( = (( B ( = (( B ( = =!!( (! (6 where =,, [,]

Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6. NEWTON OMOTOPY CONTINUATION METOD The connuaon mehod does no sand alone, we need o combne wh oher mehod and hen s called omoopy Connuaon Mehod (CM. In hs sudy, he CM wll be combned wh he classcal mehod whch s Newon mehods. Therefore, he new name s Newon omoopy Connuaon mehod. Ths mehod wll use o solve he homoopy funcon. Thus, he formula of Newon omoopy Connuaon mehod s as follow, (, + =, =,,,..., k (7 D (, where ( s homoopy funcon. 4. NUMERICAL EXPERIMENT Consder he eamples of sngle eponenal equaon: Eample : Consder he followng eample of sngle eponenal equaon ha was aken from he prevous paper [] f ( = e (8 The aulary funcon s g ( = and he nal value s =. Thus, he equaon (8 wll be mplemened no CBF as follow, ( = ( ( + ( + ( e ( + ( ( (9 The comparson resuls beween Quadrac and homoopy funcon are show n Table, Table 4 and Table 7 wh dfferen numbers of eraons. Eample : Consder he followng eample of sngle eponenal equaon ha are aken from he prevous paper [] f ( = e ( The aulary funcon s g( =. 5 and he nal value s =. 5. Thus, he equaon ( wll be mplemened no CBF as follow, ( = ( (.5 + ( ( e ( + ( ( + ( The comparson resuls beween Quadrac and omoopy funcon are shows n Table, Table 5 and Table 8 wh he dfferen sep szes of. Eample : Consder he followng eample of sngle eponenal equaon ha are aken from he prevous paper [], 4

Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 f ( = + e + ( The aulary funcon s g( = + and he nal value s =. Thus, he equaon ( wll be mplemened no CBF as follow, ( = ( ( + + ( + ( + e + ( + ( ( ( The comparson resul beween Quadrac and homoopy funcon are shows n Table, Table 6 and Table 9 wh he dfferen sep szes of. Table : Comparson of he appromae soluons beween Quadrac and omoopy Funcon for equaon (8 Sep sze, Quadrac Bezer homoopy funcon, ( homoopy funcon, (..7555.877..889.89..89.89 Roo of an equaon, =.8898764 Table : Comparson of he appromae soluons beween Quadrac and omoopy Funcon for equaon ( Sep sze, Quadrac Bezer homoopy funcon, ( homoopy funcon, (..96445.987659..97575.97575..97575.97575 Roo of an equaon, =.975748879 Table : Comparson of he appromae soluons beween Quadrac and omoopy funcon for equaon ( Sep sze, Quadrac Bezer homoopy funcon, ( homoopy funcon, (. -.8796756 -.87694. -.875884 -.87588. -.875889 -.875889 Roo of an equaon, = -.8758894968 Table, Table and Table show he resuls of he appromae soluons for equaon (8, ( and ( respecvely. A Table he real roo of an equaon s.975748879. The value of appromaon soluons a QBF column s approached o roo as he sep sze of decreased. Moreover, he appromae soluons are closer o roo of an equaon when we used CBF. Among hese wo homoopy funcons, he closes appromaons soluon o roo of an equaon s CBF wh sep szes =.. The same resuls are also appeared a Table 5

Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 and Table where he appromae soluons show a beer mprovemen when he sep sze was geng smaller. Table 4: Comparson of he accuracy of appromae soluon beween Quadrac and Cubc Bezer homoopy funcon for equaon (8 Sep szes, Quadrac Bezer homoopy funcon, (. 7 7.666. 6.4. 5 6.49 homoopy funcon, ( 7.657.88.7756 8 7 Table 5: Comparson of he accuracy of appromae soluon beween Quadrac and Cubc Bezer homoopy funcon for equaon ( Sep szes, Quadrac Bezer homoopy funcon, (.695.968. 5. 9.. homoopy funcon, (.55 5.49 4.449 Table 6: Comparson of he accuracy of appromae soluon beween Quadrac and Cubc Bezer homoopy funcon for equaon ( Sep szes, Quadrac Bezer homoopy funcon, (.885.99.9. 5. 9. homoopy funcon, (.548.894 4.449 Table 4, Table 5 and Table 6, show he resuls of he accuracy of appromae soluons for equaon (8, ( and ( respecvely. A Table 4, boh homoopy funcons show an mprovemen of he accuracy of appromae soluons when he sep szes of become smaller. Ths means, he accuracy of he appromae soluons ges beer as he value of decreased. Also from he able, CBF shows he mos precse value of he accuracy of appromae 7 soluons whch are.7756, 4.449 and 4.449 compared o QBF 5 whch are 6.49,., and.9. Ne he same resuls are shown n Table 5 and Table 6 where he value of he accuracy of appromae soluons approached o zero as he sep szes were geng smaller. All of he able showed he Cubc Bezer omoopy funcon are more accurae based on he accuracy of he appromae soluon ha approached o zero compared Quadrac Bezer homoopy funcon. Therefore, shows ha CBF s he more accurae funcon o solve hs knd of problem. 6

Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 Table 7: Comparson of he mamum absolue error of soluons beween Quadrac and homoopy funcon for equaon (8 Sep szes, Error of Soluon, l Quadrac Bezer homoopy funcon, ( homoopy funcon, (. 7 6.897 5.484 5.9.. 5 6.86 9.9476 4.449 8 4 Table 8: Comparson of he mamum absolue error of soluons beween Quadrac and homoopy funcon for equaon ( Sep szes, Error of Soluons, l Quadrac Bezer homoopy funcon, ( homoopy funcon, ( 8.87 9.96.9....94.879. Table 9: Comparson of he mamum absolue error of soluons beween Quadrac and homoopy funcon for equaon ( Sep szes, Error of Soluons, l Quadrac Bezer homoopy funcon, ( homoopy funcon, (.96 4.7 4.855.. 4. 5. 7.687.4 The comparson resuls for mamum absolue errors a Table 7, Table 8 and 9, are dscussed. All of he ables show he mamum absolue error values of soluon become smaller as he value sep szes of decreased. I ndcaes ha he value of mamum absolue errors approached o zero as he sep szes of geng smaller. Besdes, from each able above shows he mprovemen value of CBF whch are 4.449,. and.4 5 4 compared o QBF whch are 5.9,.9 and 4.855. The mamum absolue errors of soluon go smaller and converged o zero were resuled as a good and relable mehod. 7 5. CONCLUSION Ths paper has developed he eended of CBF from QBF for solvng sngle eponenal equaon n order o compue he appromae soluon, he accuracy of he appromae soluon and he mamum absolue error of soluons. From he resul obaned, we can see ha he eended funcon gve a beer, accurae and precse value of solvng hs problem compared o prevous research whch s QBF. Concluded ha, when he appromae soluon 7

Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 approached he real roo of an equaon, he accuracy of he appromae soluon ends o zero and he mamum absolue geng smaller. Thus, shows hs eended funcon s beer and more accurae han prevous funcon. REFERENCES [] Ismal, W. Zura Azan W. "Newon omoopy Connuaon Mehod for Solvng Nonlnear Equaons usng Mahemaca." Journal of Scence and Technology 7, no. (6. [] Rafq, Arf, and Muhammad Awas. "Convergence on he homoopy connuaon mehod." Inernaonal Journal of Applcaons Mahemacs and Mechancs 4, no. 6 (8: 6-7. [] Wu, Tzong-Mou. "A modfed formula of ancen Chnese algorhm by he homoopy connuaon echnque." Appled mahemacs and compuaon65, no. (5: - 5. [4] Wu, Tzong-Mou. "A sudy of convergence on he Newon-homoopy connuaon mehod." Appled Mahemacs and Compuaon 68, no. (5: 69-74. [5] Wu, Tzong-Mou. "A new formula of solvng nonlnear equaons by Adoman and homoopy mehods." Appled mahemacs and compuaon 7, no. (6: 9-97. [6] Wu, Tzong-Mou. "Solvng he nonlnear equaons by he Newon-homoopy connuaon mehod wh adjusable aulary homoopy funcon." Appled mahemacs and compuaon 7, no. (6: 8-88. [7] Wu, Tzong-Mou. "The secan homoopy connuaon mehod." Chaos, Solons & Fracals, no. (7: 888-89. [8] Nor, afzudn Mohamad, Ahmad Izan Md Ismal, and Ahmad Abdul Majd. "Quadrac bezer homoopy funcon for solvng sysem of polynomal equaons." Maemaka 9 (: 59-7. [9] Farn, Gerald. Curves and surfaces for compuer-aded geomerc desgn: a praccal gude. Elsever, 4. [] Agoson, Ma K. Compuer graphcs and geomerc modelng. Vol.. London: Sprnger, 5. [] Chun, Changbum. "A new erave mehod for solvng nonlnear equaons."appled mahemacs and compuaon 78, no. (6: 45-4. [] Fang, Lang, Tao Chen, L Tan, L Sun, and Bn Chen. "A modfed Newon-ype mehod wh sh-order convergence for solvng nonlnear equaons."proceda Engneerng 5 (: 4-8. 8