wavelet collocation method for solving integro-differential equation.

Similar documents
Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Numerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets

Chapter 4. Fourier Series

Numerical solution of Bagley-Torvik equation using Chebyshev wavelet operational matrix of fractional derivative

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

μ are complex parameters. Other

Numerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION

Numerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

Modified Decomposition Method by Adomian and. Rach for Solving Nonlinear Volterra Integro- Differential Equations

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

NTMSCI 5, No. 1, (2017) 26

Ma 530 Introduction to Power Series

Chapter 10: Power Series

Solution of Differential Equation from the Transform Technique

Castiel, Supernatural, Season 6, Episode 18

CHEBYSHEV WAVELET COLLOCATION METHOD FOR GINZBURG-LANDAU EQUATION

NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS

Notes on iteration and Newton s method. Iteration

Research Article A New Second-Order Iteration Method for Solving Nonlinear Equations

A) is empty. B) is a finite set. C) can be a countably infinite set. D) can be an uncountable set.

f t dt. Write the third-degree Taylor polynomial for G

Teaching Mathematics Concepts via Computer Algebra Systems

The Phi Power Series

Subject: Differential Equations & Mathematical Modeling -III. Lesson: Power series solutions of Differential Equations. about ordinary points

Two-step Extrapolated Newton s Method with High Efficiency Index

6.3 Testing Series With Positive Terms

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

IN many scientific and engineering applications, one often

Newton Homotopy Solution for Nonlinear Equations Using Maple14. Abstract

x x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0,

THE TRANSFORMATION MATRIX OF CHEBYSHEV IV BERNSTEIN POLYNOMIAL BASES

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations

Example 2. Find the upper bound for the remainder for the approximation from Example 1.

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations

Chapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics:

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method

( a) ( ) 1 ( ) 2 ( ) ( ) 3 3 ( ) =!

Chapter 7: Numerical Series

Quiz. Use either the RATIO or ROOT TEST to determine whether the series is convergent or not.

Power Series: A power series about the center, x = 0, is a function of x of the form

Generating Functions for Laguerre Type Polynomials. Group Theoretic method

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Poincaré Problem for Nonlinear Elliptic Equations of Second Order in Unbounded Domains

Solving a Nonlinear Equation Using a New Two-Step Derivative Free Iterative Methods

MATH301 Real Analysis (2008 Fall) Tutorial Note #7. k=1 f k (x) converges pointwise to S(x) on E if and

Chapter 6 Infinite Series

Some Variants of Newton's Method with Fifth-Order and Fourth-Order Convergence for Solving Nonlinear Equations

The Method of Least Squares. To understand least squares fitting of data.

x x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy

Some New Iterative Methods for Solving Nonlinear Equations

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis

Research Article Two Expanding Integrable Models of the Geng-Cao Hierarchy

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

Polynomials with Rational Roots that Differ by a Non-zero Constant. Generalities

PAPER : IIT-JAM 2010

The Differential Transform Method for Solving Volterra s Population Model

9.3 Power Series: Taylor & Maclaurin Series

THE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.

Diploma Programme. Mathematics HL guide. First examinations 2014

Advanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology

Math 113 Exam 3 Practice

Chapter 8. Euler s Gamma function

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

CHAPTER 10 INFINITE SEQUENCES AND SERIES

Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations

TEACHING THE IDEAS BEHIND POWER SERIES. Advanced Placement Specialty Conference. LIN McMULLIN. Presented by

Subject: Differential Equations & Mathematical Modeling-III

Chapter 8. Uniform Convergence and Differentiation.

Higher-order iterative methods by using Householder's method for solving certain nonlinear equations

-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION

CS321. Numerical Analysis and Computing

MATHEMATICS. 61. The differential equation representing the family of curves where c is a positive parameter, is of

Properties of Fuzzy Length on Fuzzy Set

A GENERALIZED CHEBYSHEV FINITE DIFFERENCE METHOD FOR HIGHER ORDER BOUNDARY VALUE PROBLEMS

SCORE. Exam 2. MA 114 Exam 2 Fall 2016

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

CS537. Numerical Analysis and Computing

FUZZY ALTERNATING DIRECTION IMPLICIT METHOD FOR SOLVING PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS IN THREE DIMENSIONS

Math 113, Calculus II Winter 2007 Final Exam Solutions

Lainiotis filter implementation. via Chandrasekhar type algorithm

Numerical Solution of the First-Order Hyperbolic Partial Differential Equation with Point-Wise Advance

THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT

Numerical Methods in Fourier Series Applications

PC5215 Numerical Recipes with Applications - Review Problems

Seunghee Ye Ma 8: Week 5 Oct 28

1 6 = 1 6 = + Factorials and Euler s Gamma function

The Sumudu transform and its application to fractional differential equations

f(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction.

EE / EEE SAMPLE STUDY MATERIAL. GATE, IES & PSUs Signal System. Electrical Engineering. Postal Correspondence Course

MATH 324 Summer 2006 Elementary Number Theory Solutions to Assignment 2 Due: Thursday July 27, 2006

Research Article Laplace Transform Collocation Method for Solving Hyperbolic Telegraph Equation

Abstract Vector Spaces. Abstract Vector Spaces

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

8. Applications To Linear Differential Equations

Ma 4121: Introduction to Lebesgue Integration Solutions to Homework Assignment 5

Average Number of Real Zeros of Random Fractional Polynomial-II

Transcription:

IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma Uiversity of Techology Applied Sciece Departmet Abstract: - Wavelet collocatio method for umerical solutio th order Volterra itegro diferetial equatios (VIDE) by epadig the ukow fuctios, as series i terms of chebyshev wavelets secod kid with ukow coefficiets. The aim of this paper is to state ad prove the uiform covergece theorem ad accuracy estimatio for series above. Fially, some illustrative eamples are give to demostrate the validity ad applicability of the proposed method. Keywords: chebyshev wavelets secod kid; itegro-differetial equatio; operatioal matri of itegratios; uiform covergece; accuracy estimatio. I. INTRODUCTION Basic wavelet theory is a atural topic. By ame, wavelets date back oly to the 98s. o the boudary betwee mathematics ad egieerig, wavelet theory shows studets that mathematics research is still thrivig, with importat applicatios i areas such as image compressio ad the umerical solutio of differetial equatios [], itegral equatio [],ad itegro differetial equatios [3,8]. The author believes that the essetials of wavelet theory are sufficietly elemetary to be taught successfully to advaced udergraduates [4]. Orthgoal fuctios ad polyomials series have received cosiderable attetio i dealig with various problems wavelets permit the accurate represetatio of a fuctios ad operators. Special attetio has bee give to applicatio of the Legedre wavelets [5], Harr wavelets [6] ad Sie-Cosie wavelets [7]. The solutio of itegro-differetial equatios have a major role i the fields of sciece ad igieerig whe aphysical system is modeled uder the differetial sese, it fially gives a differetial equatio, a itegral equatio or a itegro-differetial equatios mostly appear i the last equatio [8]. I this paper the operatioal matri of itegratio for Hermite wavelets is derived ad used it for obtaiig approimate solutio of the followig th order VIDE. u () ()g()+ k, t u (s) t dt () where s k(,t) ad g() are kow fuctios, ad u() is a ukow fuctio. II.SOE PROPERTIES OF SECOND CHEBYSHEV WAVELETS Wavelets costitute a family of fuctios costructed from dilatio ad traslatio of a sigle fuctio t called the mother wavelet.whe the dilatio parameter a ad the traslatio parameter b vary cotiuously we have the followig family of cotiuous wavelets as [9]. a,b t a tb, a, b R, a. a The secod chebyshev wavelets,m t (k,, m, t) ivolve four argumets,,,, k is assumed ay positive iteger, m is the degree of the secod chebyshev polyomials ad t is the ormalized time. They are defied o the iterval,) as,m t k U m k t +, k t < () oterwise where U m t U m t m,,, (3) here U m t are the secod chebyshev polyomialsof degree m with respect to the weight fuctio w t t o the iterval [-,] ad satisfy the followig recursive formula U t, U t t, U m + t tu m t U m t, m,, The set of chebyshev wavelets are a orthoirmal set with respect to the weight fuctio w t w ( k t + ). Iteratioal orgaizatio of Scietific Research P a g e

III.FUNCTION APPROXIAT A fuctio f(t) defied over [, ) may be epaded as wavelet collocatio method for solvig itegro-differetial equatio. f t m C m m (t) (4) were C m (f t, m t ) I which.,. deoted the ier product i L w [,). If the ifiite series equatio (.34) is trucated, the it ca be writte f t k m C m m t C T (t) (5) C C, C,, C (), C,, C (),, C,, C T (t) (t), (t),, (t), (t),, for Chebyshev wavelets of Secod kid. A fuctio f t defied over, may be epaded as: f t m f m m t (t), (6) (t), (t) T Covergece Aalysis where f m f t, m t (7) I eq.(5)..,. deotes the ier product with weight fuctio w t. If the ifiite series i eq.(4) is tracated the eq.(4) ca be writte as: f t f t f m m m F T m where F ad m are k matrices give by F f, f,, f, f,, f,, f k,, f T (8) (t) t, t,, t, t,, t, t,, t T Theorem(): (Covergece Aalysis theorm) Assume that a fuctio f(t) L w,, w t with t L, ca be epaded as ifiite series of secod kid chebyshev wavelets, the the series coverges uiformly to f(t). Proof: sice f m f t, m (t) the f m f(t) m (t)w t dt k+ f(t)u m k t + w k t + dt If we make use of the substitutio k t + cos θ i (5), yields cos θ + f m k f k si θ si m + θ By usig the itegratio by parts, The eq.() becomes cos θ + si mθ si m + θ f m k f k m m + cos θ + + 3k f k si θ cos θ + f m m 3k f k si θ si mθ m + 3k f After performig itegratio by parts agai yields, si mθ m (9) () si m + θ m + cos θ + k si θ si m + θ () Iteratioal orgaizatio of Scietific Research P a g e

Cosider: f m m 3 5k wavelet collocatio method for solvig itegro-differetial equatio. cos θ + cos m θ cos mθ cos mθ cos m + θ k m m + cos θ + cos m θ cos m + θ m + 3 5k k m + cos m + θ cos m + 4 θ m + 3 () cos θ + cos m θ cos mθ cos mθ cos m + θ k m m + cosθ + cos m θ cos mθ cos mθ cos m + θ k m m + cosθ + m cos m + θ m cos m θ + m + cos m θ k m m + m cos m + θ + 4m cos mθ + m + cos m < L m m + L m m + m + 4m + m + L m m + 3m + Thus, we get cos θ + cos m θ cos mθ cos mθ cos m + θ k m m + < L 3m + (3) m m + Similarly, cos θ + cos m θ cos m + θ cos m + θ cos m + 4 θ k m + m + 3 cos θ + k m + cos m + 4 θ (m + ) cos m + θ + m + 3 cos m + θ + m + 3 cos m θ (m + ) m + 3 cos θ + k Thus we get < L Usig eqs.(3)ad(4), oe ca get m cos m + 3 θ m cos m + θ + m + θ cos m θ m m + m + cos m + 4 θ + m + 4 cos m + θ + m + 3 cos m θ m + m + 3 L m + m + 3 6m + 4m + 6 L m + m + 3 3m + m + 3 cos θ + cos m θ cos m + θ k m + < L m + m + 3 3m + m + 3 cos m + θ cos m + 4 θ m + 3 (4) Iteratioal orgaizatio of Scietific Research 3 P a g e

5k f m < 3 m + wavelet collocatio method for solvig itegro-differetial equatio. f m < 5k3 L 3m + m m m + L 3m + m + 3 m + m + m + 3 m+ 5k L m m Fially sice k,the f m < L + 5 m Accuracy Estimatio of m (): If the fuctio f() is epaded iterms of fourth kid chebyshev wavelets, f m f m m () (5) It is ot possible to perform computatio a ifiite umber of terms, therfore we must trucate the series i (5). I place of (5), we take f k m f m m () (6) so that f f + + m f m m () or f f r() (7) where r() is the residual fuctio r + m f m m () (8) we must select coefficiets i eqs.(7) ad (8) such that the orm of the residual fuctio r() is less tha some covergece criterio,that is for all greater tha some value. f f w d < Theorem () Let f() be a cotiuaus fuctio defied o [,), ad f () < L, the we have the followig accuracy estimatio c k, < L + m (9) + 5 m 4 where Proof:- c k, r w d Sice r w d + r w d f m w d m + m f m m () w d or f m w d + m f k + m U m k + k + d Let t k +, f U m t t d we have, the + < m U m t t d m f + L m + 5 + m. 8 m 4 Iteratioal orgaizatio of Scietific Research 4 P a g e

wavelet collocatio method for solvig itegro-differetial equatio. wavelet collocatio method for VIDE with mth order: I this sectio the itroduced wavelets collocatio will be applied to solve VIDE with mth order, ( u ) i g i + K i,j, t u (s) i t dt, s () s With the followig coditios u i a is i,,, l s,,,, Afuctio u i which is defied o the iterval, ca be epaded ito the secod chebyshev wavelet series u i i c i i (t) () Where c i are the wavelet coefficiets. Itegrate eq.() m times,yields m j u i c i i t dt + j a j! m j () Usig the followig formula i t dt t ( )! i (t)dt therefore eq.() becomes u c i t i (t)dt + i j (3) This leads to Let K, t i u c i L i + j I similar way, we ca get t! j j! a j ()! ad L i K, t i t dt s j j! a j i,,, u (s) i c i L s i + j a j! sj (4) Substitutig eqs () ad (4) i (), yield or i c i i (t) g i + K i,j, t c i L s i + i c i i (t) A i g i + where A i K, t L i s t dt B j K, t t j dt s j i j a j! sj dt s a sj B j! j () (5) j (6) i,,,, j,,,,-s- (7) Net the iterval, is devided i to l ad itroduce the collocatio poits l k k, k,,,l eq() is satisfied oly at the collocatio poits we get asystem of liear equatios l s a sj i c i [ i () A i ] g i + j j! B j () (8) The matri form of this system s a sj is C FG+ j B j () where F (), Gg() j!.desig of the matri A:- Whe chebyshev wavelets secod kid are itegrated m times, the followig itegral must be evaluated. L i K, t i t dt, i,,,, L i t k! 3 4 3 6 4 () Therefore the matri A i ca be costructed as follows s Sice A i K, t L i t dt i,,,, l k < l k Iteratioal orgaizatio of Scietific Research 5 P a g e

wavelet collocatio method for solvig itegro-differetial equatio. A i s K, t L i t dt i K i, t L i s t dt i > IV. WAVELETS ETHOD FOR VIDE WITH NTH ORDER For solvig VIDE with mth order the matri L i i sectio above will be followed to get s a sj i c i [ i ( L A L )] g L + j B j ( L ) L a, b But A i L L L j! K L, t L i s t dt where i,, B j ( L ) K L, t t s dt s where L i t as i eq(7),(8) that is A i L A L, F i L i L A i L F L Numerical Results: I this sectio VIDE is cosidered ad solved by the itroduced method. parameters k ad are cosidered to be ad 3 respectively. Eample : Cosider the followig VIDE: U e e (t) U t dt Iitial coditios U(), U (). The eact solutio U e e +. Table shows the umerical results for this eample with k, 3 with error -3 ad k, 4, with error -4 are compared with eact solutio graphically i fig. Table :some umerical results for eample Eact solutio Approimat solutio k,3 Approimat solutio k,4.....87779.8.87.4.9458.945544.94544.6.7548.586756.786756.8.554898.5433957.5533957..99999995.99999998.9.8.7.6.5.4.3.....3.4.5.6.7.8.9 Fig :Approimate solutio for eample Eample : Cosider the followig VIDE : U (5) si + cos + ( t)u (3) t dt Iitial coditios U(), U (), U"(), U 3 (). The eact solutio U cos. Table shows the umerical results for this eample with k, 3 with error -3 ad k, 4, with error -4 are compared with eact solutio graphically i fig,. Iteratioal orgaizatio of Scietific Research 6 P a g e

wavelet collocatio method for solvig itegro-differetial equatio. Table :some umerical results for eample Eact solutio Approimat solutio k,3 Approimat solutio k,4..99835.99999875..986658.9847.98554.4.9699.95899.9436.6.853356.84798.8535367.8.696767.6968963.69678976.5433.543879.5435879.95.9.85.8.75.7.65.6.55.5...3.4.5.6.7.8.9 Fig :Approimate solutio for eample V. CONCLUSION This work proposes a powerful techique for solvig VIDE secod kid usig wavelet i collocatio method compariso of the approimate solutios ad the eact solutios shows that the proposed method is more faster algorithms tha ordiary oes. The covergece ad accuracy estimatio of this method was eamied for several umerical eamples. REFERENCES [] Asmaa A. A, 4, Numerical, Solutio of Optimal Cotrol Problems Usig New Third kid Chebyshev Wavelets Operatioal atri of Itegratio, Eg&Tech joural,vol 3,part(B),:45-56. [] Tao. X. ad Yua. L.. Numerical Solutio of Fredholm Itegral Equatio of Secod kid by Geeral Legedre Wavelets, It. J. I. Comp ad Cot. 8(): 799-85. [3] A.Barzkar,.K.Oshagh,, Numerical solutio of the oliear Fredholm itegro-differetial equatio of secod kid usig chebysheve wavelets, World Applied Scices Joural (WASJ),Vol.8, N():774-78. [4] E.Johasso, 5, Wavelet Theory ad some of its Applicatios. [5] Jafari. H ad Hosseizadeh. H.. Numerical Solutio of System of Liear Itegral Equatios by usig Legedre Wavelets, It. J. Ope Problems Compt. ath., 3(5): 998-66. [6] Shihab. S. N. ad ohammed. A.. A Efficiet Algorithm for th Order Itegro- Differetial Equatios Usig New Haar Wavelets atri Desigatio, Iteratioal Joural of Emergig & Techologies i Computatioal ad Applied Scieces (IJETCAS). (9): 3-35. [7].Razzaghi & S.Yousefi.. Si-Cosie wavelets operatioal matri of itegratio ad it s applicatios i the calculus of variatios. Vol 33. No : 85-8. [8] A.Arikoglu & I.Ozkol. 8. Solutio of itegral itegro-differetial equatio systems by usig differetial trasform method. Vol 56,Issue 9.4-47. [9] Arsalai.. & Vali.. A., Numerical Solutio of Noliear Problems With ovig Boudary Coditios by Usig Chebyshev Wavelets, Applied athematical Scieces,Vol.5(): 947-964. Iteratioal orgaizatio of Scietific Research 7 P a g e