A Note on the Numerical Solution for Fredholm Integral Equation of the Second Kind with Cauchy kernel

Similar documents
A new Approach for Solving Linear Ordinary Differential Equations

Convexity preserving interpolation by splines of arbitrary degree

Solution for singularly perturbed problems via cubic spline in tension

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

Explicit bounds for the return probability of simple random walk

High-Order Hamilton s Principle and the Hamilton s Principle of High-Order Lagrangian Function

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

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Chapter 3 Differentiation and Integration

Chapter 2 Transformations and Expectations. , and define f

DEGREE REDUCTION OF BÉZIER CURVES USING CONSTRAINED CHEBYSHEV POLYNOMIALS OF THE SECOND KIND

arxiv: v1 [math.co] 12 Sep 2014

Large-Scale Data-Dependent Kernel Approximation Appendix

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

New Liu Estimators for the Poisson Regression Model: Method and Application

Visualization of 2D Data By Rational Quadratic Functions

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

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

The Finite Element Method: A Short Introduction

ENTROPIC QUESTIONING

WHY NOT USE THE ENTROPY METHOD FOR WEIGHT ESTIMATION?

Global Sensitivity. Tuesday 20 th February, 2018

Construction of Serendipity Shape Functions by Geometrical Probability

A FORMULA FOR COMPUTING INTEGER POWERS FOR ONE TYPE OF TRIDIAGONAL MATRIX

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Appendix B. The Finite Difference Scheme

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

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

Beyond Zudilin s Conjectured q-analog of Schmidt s problem

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

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

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

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Inductance Calculation for Conductors of Arbitrary Shape

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON

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

A Hybrid Variational Iteration Method for Blasius Equation

Lecture 12: Discrete Laplacian

Numerical Heat and Mass Transfer

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

p(z) = 1 a e z/a 1(z 0) yi a i x (1/a) exp y i a i x a i=1 n i=1 (y i a i x) inf 1 (y Ax) inf Ax y (1 ν) y if A (1 ν) = 0 otherwise

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

PHZ 6607 Lecture Notes

Maejo International Journal of Science and Technology

A New Recursive Method for Solving State Equations Using Taylor Series

Another converse of Jensen s inequality

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

6.3.4 Modified Euler s method of integration

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

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

A MULTIDIMENSIONAL ANALOGUE OF THE RADEMACHER-GAUSSIAN TAIL COMPARISON

Analytical classical dynamics

risk and uncertainty assessment

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

Report on Image warping

A GENERALIZATION OF JUNG S THEOREM. M. Henk

Numerical Solution of Singular Perturbation Problems Via Deviating Argument and Exponential Fitting

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

APPENDIX A Some Linear Algebra

Haar Wavelet Collocation Method for the Numerical Solution of Nonlinear Volterra-Fredholm-Hammerstein Integral Equations

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

2.29 Numerical Fluid Mechanics

CHAPTER 4d. ROOTS OF EQUATIONS

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

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

Statistics Chapter 4

Kinematics of Fluid Motion

Deriving the X-Z Identity from Auxiliary Space Method

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

1.050 Content overview Engineering Mechanics I Content overview. Outline and goals. Lecture 28

Denote the function derivatives f(x) in given points. x a b. Using relationships (1.2), polynomials (1.1) are written in the form

The Finite Element Method

An Inequality for the trace of matrix products, using absolute values

An efficient algorithm for multivariate Maclaurin Newton transformation

Section 3.6 Complex Zeros

Journal of Mathematical Analysis and Applications

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

Integrals and Invariants of Euler-Lagrange Equations

Time-Varying Systems and Computations Lecture 6

Upwind schemes for the wave equation in second-order form

The Jacobsthal and Jacobsthal-Lucas Numbers via Square Roots of Matrices

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

On a one-parameter family of Riordan arrays and the weight distribution of MDS codes

Numerical modeling of a non-linear viscous flow in order to determine how parameters in constitutive relations influence the entropy production

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

2nd International Conference on Electronics, Network and Computer Engineering (ICENCE 2016)

LETTER A Mathematical Proof of Physical Optics Equivalent Edge Currents Based upon the Path of Most Rapid Phase Variation

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

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

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

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Bounds for Spectral Radius of Various Matrices Associated With Graphs

On the First Integrals of KdV Equation and the Trace Formulas of Deift-Trubowitz Type

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

Lecture 2 Solution of Nonlinear Equations ( Root Finding Problems )

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

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours

Transcription:

Journal of Mathematcs an Statstcs 7 (): 68-7, ISS 49-3644 Scence Publcatons ote on the umercal Soluton for Freholm Integral Equaton of the Secon Kn wth Cauchy kernel M. bulkaw,.m.. k Long an Z.K. Eshkuvatov Department of Mathematcs an Insttute for Mathematcal Research, Unversty Putra Malaysa, 434 Serang, Selangor, Malaysa bstract: Problem statement: In ths stuy, numercal soluton for the Freholm ntegral equaton of the secon kn wth Cauchy sngular kernel s presente. pproach: The Chebyshev polynomals of the secon kn are use to appromate the unknown functon. Results: umercal results are gven to show the accuracy of the present numercal soluton. Concluson: The present numercal soluton to the Freholm ntegral equaton of the secon kn wth Cauchy kernel s accurate. Key wors: Integral equatons, Cauchy sngular kernel, Chebyshev polynomal appromatons ITRODUCTIO Conser the Freholm ntegral equaton of the secon kn wth Cauchy kernel of the form: ϕ(t) ϕ () +λ t = f (), ( λ s a constant) () wth specfe en contons: ϕ± ( )= () where f s assume to be real-value functons belong to the class of Holer on the nterval [-, ] an ϕ s the unknown functon to be etermne. The sngular ntegral n Eq. beng unerstoo n the sense of Cauchy prncpal value. Many researchers evote ther works on appromatng the sngular ntegrals of Cauchy type. Partcularly, Dagnno an Sant (99) obtane a prouct quarature rules, base on splne nterpolaton, for the numercal evaluaton of Cauchy sngular ntegrals. They also propose an error boun an obtane convergence results for functons f C k [-,], k=, or 3. Ors (99) prove the unform convergence of some quarature formulas base on splne appromaton for Cauchy prncpal value ntegral. They also presente some numercal applcatons. In partcular, they apple ther rules to the well-known Prantl s ntegral equaton. Rabnowtz (99) prove convergence results for prouct ntegraton rules base on appromatng splnes. These results are both for boune an unboune ntegrans. Pontwse an unform convergence results are prove for sequences of Cauchy prncple values of these appromatng sples..dagnno an Sant (99) consere the same rules as n Dagnno an Sant (99) an nvestgate ther convergence for a large class of functons f. They establshe an error boun an some unform convergence results n the case of equally space quarature noes, for functon f, satsfyng a Holer conton of orer μ on [,], < μ. Hasegawa an Tor (99) presente an automatc quarature for computng Cauchy prncpal value ntegrals for smooth functons f(t). They appromate the functon f(t) by a sum of Chebyshev polynomals whose coeffcents are compute usng the Fast Fourer Transform (FFT). Hasegawa an Tor (994) presente an automatc quarature for appromatng Haamar fnte-part (fp) ntegrals of a smooth functon, wth a ouble pole sngularty wthn the range of ntegraton. The quarature rule s rve from the fferentaton of an appromaton to a Cauchy prncpal value ntegral. Dethelm (99) nvestgate the numercal appromaton of the Cauchy prncpal value ntegral. He presente the quarature formula for appromatng the Cauchy prncpal value ntegral. He prove the convergence of the quarature formula an gave the estmaton for the errors. Dagnno an Lambert (996) evaluate the Cauchy prncpal value ntegral by applyng a local splne appromaton metho, efne for any functon f L [,]. They establshe convergence results wth error boun. Dethelm (997) consere the so-calle mofe quarature formulas,.e. formulas obtane by frst subtractng out the sngularty an then applyng a classcal quarature Corresponng uthor: Dr. Mohamma bulkaw, Department of Mathematcs an Insttute for Mathematcal Research, Unversty Putra Malaysa, 434 Serang, Selangor, Malaysa 68

J. Math. & Stat., 7 (): 68-7, formula, for the numercal appromaton of Cauchy prncpal value ntegrals. They have gven new bouns nvolvng the total varaton Var f (s) an L p -norms f (s) p of some ervatve of the ntegran functon. Eshkuvatov et al. (9) constructe a new quarature formulas for evaluatng the sngular ntegral of Cauchy type. The constructon of the quarature formulas s base on the mofcaton of screte vortces metho an lnear splne nterpolaton over the fnte nterval [, ]. They prove that the constructe quarature formula converges for any sngular pont not concng wth the en ponts of the nterval [, ]. They have gven error bouns n the classes of functons H α [, ] an C [, ] of orer O(h α ln h ), < α an O(h ln h ), respectvely. Eshkuvatov et al. () consere the sngular ntegral wth the Cauchy kernel. They constructe new quarature formulas base on the mofcaton of screte vorte metho to appromate the sngular ntegral. They have shown error bouns n the classes of functons H α [,] an C [,] for ether = t or = t where t,,..., are the noe ponts an t =(t + t + ) /. The ntegral equatons wth Cauchy kernel have been wely use n solvng problems assocate wth aeroynamc, hyroynamc an elastcty (Lfanov, 996; Laopoulos, ; bou an aser, 3; Mohankumar an ataraan, 8; Lara an Maragraza, ; Kasoz an Paulsen, a; Kasoz an Paulsen, b; Gan et al., 8; Thukral, ). In ths stuy, we present a numercal soluton for the Eq. wth contons (). MTERILS D METHODS The unknown functon ϕ of Eq. whch satsfes contons () can be represente as: ϕ()= ψ(), (3) where, ψ () s a well behave functon of on the nterval [-, ]. The functon ψ () n Eq. 3 s appromate usng the Chebyshev polynomals of the secon kn, U, as: From (3) an (4), we have: () au (), ϕ (6) Substtutng (6) nto (), yels: au () +λ a t U (t) t =f() (7) It s known that (Kythe an Schaferkotter, ; bulkaw et al., 9) t U (t) t = π T + () (8) where, T are the Chebyshev polynomals of the frst kn whch s efne by: T () = cos cos (), =,,, (9) Usng (8) nto Eq. 7, yels: a U () πλt + () =f() () Multplyng both ses of Eq. by U an ntegratng from - to, we obtan the followng system of lnear equatons; a =, =,,, () where: = B + C () π, = B = t U (t)u (t)t = (3), ψ () where: a U () (4) C = πλ T (t)u (t)t (4) an: + sn[ cos ()] U () =,,, () sn[ cos ()] 69 = f(t)u (t)t ()

J. Math. & Stat., 7 (): 68-7, Solvng the above system for the unknown coeffcents a, =,,, an substtutng the values of a nto the appromate soluton (6), we obtan the numercal soluton of the Eq. wth contons (). RESULTS D DISCUSSIO Eample : Let us conser the sngular ntegral equaton: ϕ(t) (6) ϕ () + t = π + ( π ) + +π wth the contons: ϕ± ( )= (7) It s not ffcult to see that the eact soluton of the Eq.6 s: φ() = ( + ) (8) Due to Eq. -4, we obtan: π = π, =, 3 4 π 4 = π, = π 3 6 π =, = π, =. (9) π a + πa = π+, 3 3 4 4 4 π π a + a + πa = π+, 3 3 6 6 π a+ a = π. () It s easy to see that the soluton to the above system s: a a a = (3) Substtutng (3) nto (6) for we obtan: ϕ() = (+ ) (4) whch s the eact soluton. Eample : Conser the followng equaton: ϕ(t) ϕ () + t = U () πt 6() () wth the contons: ϕ ( ± )= (6) It s clear that the eact soluton of the Eq. s: 3 ϕ() = (3 3 + 6) (7) From Eq. -4, we obtan: From Eq. (), we have: ( π + ( π) ) ( t ) U (t)t = t t t U (t)t + +π whch gves: = π+ π, 3 4 = π+ π, 3 6 = π. () () Thus the system of lnear equatons () for = becomes: 7 =, 3 =, 3 π 4 4 =, 3 8 =, 3 4 =, 3 π 6 3 =, 3 4 4 =, 3 =, 3 =, π 8 4 3 =, 33 43 =, 3 7 9 6 π 4 4 =, 4 34 7 4 44 4 =, =, =, 9 8 π =, 3 =, 4 =. 9 (8)

Due to Eq., we have: π 6 = ( t U (t) T (t))u (t)t (9) whch gves: =, =, 3 4 3 4 =, =. 9 J. Math. & Stat., 7 (): 68-7, (3) Thus the corresponng system of lnear equatons for = s: a =, =,,, (3) where,,, =,,, are efne n (8) an, =,,, are efne n (3). It s not ffcult to see that the soluton to the above system s: a =,a =,a =,a =,a =,a = (3) 3 4 Substtutng the values of a, =,,, nto (6) where =, we obtan the numercal soluton of Eq. whch s entcal to the eact soluton (7). COCLUSIO The truncate seres nvolvng the Chebyshev polynomal appromaton of the secon kn s use to appromate the unknown functon for solvng the Freholm ntegral equaton of the secon kn wth Cauchy kernel an constant coeffcents. The metho of appromaton llustrate here gves a goo way for obtanng the numercal soluton avong the complcate ntegratons. umercal results show the accuracy of the metho presente whch, for some functons f (), gves the eact soluton. CKOWLEDGEMET Ths stuy was supporte by Unversty Putra Malaysa uner Postoctoral Fellowshp an proect o. -3--988RU. REFERECES bou, M.,. an.. aser, 3. On the numercal treatment of the sngular ntegral equaton of the secon kn. pple Math. Computaton, 46: 373-38. DOI:.6/S96-33()87-8 7 bulkaw, M., Z.K. Eshkuvatov an.m... Long, 9. numercal treatment for solvng Cauchy type sngular ntegral equaton. FJM., 3: 7-8. Dagnno, C. an E. Sant, 99. Splne prouct quarature rules for cauchy sngular ntegrals. J. Comput. pple Math., 33: 33-4. DOI:.6/377-47(9)9363- Dagnno, C. an E. Sant, 99. On the convergence of splne prouct quaratures for Cauchy prncpal value ntegrals. J. Comput. pple Math., 36: 8-87. DOI:.6/377-47(9)9-F Dagnno, C. an P. Lambert, 996. umercal evaluaton of Cauchy prncpal value ntegrals base on local splne appromaton operators. J. Comput. pple Math., 76: 3-38. DOI:.6/S377-47(96)-7 Dethelm, K., 99. Gaussan quarature formulae of the thr kn for Cauchy prncpal value ntegrals: Basc propertes an error estmates. J. Comput. pple Math., 6: 97-4. DOI:.6/377-47(9)3-4 Dethelm, K., 997. ew error bouns for mofe quarature formulas for Cauchy prncpal value ntegrals. J. Comput. pple Mathe., 8: 93-4. DOI:.6/S377-47(97)4-9 Eshkuvatov, Z.K.,.M... Long an M. bulkaw, 9. Quarature formula for appromatng the sngular ntegral of Cauchy type wth unboune weght functon on the eges. J. Comput. pple Math., 33: 334-34. DOI:.6/.cam.9.7.34 Eshkuvatov, Z.K.,.M... Long an M. bulkaw,. umercal evaluaton for Cauchy type sngular ntegrals on the nterval. J. Comput. pple Math., 33: 99-. DOI:.6/.cam.9.9.34 Gan, D.D., M. mn an. Kolahooz, 8. nalytcal nvestgaton of Hyperbolc equatons va he s methos. m. J. Eng. pple Sc., : 399-47. DOI:.3844/aeassp.8.399.47 Hasegawa, T. an T. Tor, 99. n automatc quarature for Cauchy prncpal value ntegrals. Mathe. Comput., 6: 74-4. Hasegawa, T. an T. Tor, 994. Hlbert an haamar transforms by generalze chebyshev epanson. J. Comput. pple Math., : 7-83. DOI:.6/377-47(9)3-R Kasoz, J. an J. Paulsen, a. Flow of vens uner a constant force of nterest. m. J. pple Sc., : 389-394. DOI:.3844/aassp.. 389.394 Kasoz, J. an J. Paulsen, b. umercal ultmate run probabltes uner nterest force. J. Math. Stat., : 46-. DOI:.3844/mssp..46.

J. Math. & Stat., 7 (): 68-7, Kythe, K.P. an R.M. Schaferkotter,. Han-book of computatonal Methos for ntegraton. st En., Chapman an Hall /CRC Press, Lonon, ISB-: 848848, pp: 98. Laopoulos, E.,. Sngular ntegral equatons: Lnear an on-lnear Theory an ts applcatons n Scence an Engneerng. st En., Sprnger Verlag, US., ISB-: 97834673, pp: 84. Lara, F. an G. Maragraza,. Parsmonous var moels for ar polluton ynamc analyss. J. Math. Stat., :8-67. DOI:.3844/mssp..8.67 Lfanov, I. K., 996. Sngular Integral Equaton an Dscrete Vortces. st En., Walter e Gruyter, Croata, ISB-: 967647X, pp: 47. Mohankumar,. an. ataraan, 8. On the numercal soluton of cauchy sngular ntegral equatons n neutron transport. nnals uclear Energy, 3: 8-84. DOI:.6/.anucene.8..7 Ors,.P., 99. Splne appromaton for Cauchy prncpal value ntegrals. J. Comput. pple Math., 3: 9-. DOI:.6/377-47(9)97-W Rabnowtz, P., 99. umercal ntegraton base on appromatng splnes. J. Comput. pple Math., 33: 73-83. DOI:.6/377-47(9)97-Z Thukral, R.,. Introucton to the new type of algorthms for acceleratng convergence of sequence. J. Math. Stat., : -3. DOI:.3844/mssp...3 7