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

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

6.3 Testing Series With Positive Terms

Some properties of Boubaker polynomials and applications

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 +

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Math 113, Calculus II Winter 2007 Final Exam Solutions

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

Solution of Differential Equation from the Transform Technique

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

wavelet collocation method for solving integro-differential equation.

Find quadratic function which pass through the following points (0,1),(1,1),(2, 3)... 11

Chapter 10: Power Series

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.

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

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

THE TRANSFORMATION MATRIX OF CHEBYSHEV IV BERNSTEIN POLYNOMIAL BASES

Some New Iterative Methods for Solving Nonlinear Equations

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

Additional Notes on Power Series

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.

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

Fall 2013 MTH431/531 Real analysis Section Notes

In this section, we show how to use the integral test to decide whether a series

MATH4822E FOURIER ANALYSIS AND ITS APPLICATIONS

MAT1026 Calculus II Basic Convergence Tests for Series

A multivariate rational interpolation with no poles in R m

Section 11.8: Power Series

1. (25 points) Use the limit definition of the definite integral and the sum formulas 1 to compute

The Numerical Solution of Singular Fredholm Integral Equations of the Second Kind

Archimedes - numbers for counting, otherwise lengths, areas, etc. Kepler - geometry for planetary motion

About the use of a result of Professor Alexandru Lupaş to obtain some properties in the theory of the number e 1

On a class of convergent sequences defined by integrals 1

Lecture 3: Convergence of Fourier Series

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

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

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

IN many scientific and engineering applications, one often

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations

Estimation of Backward Perturbation Bounds For Linear Least Squares Problem

A new iterative algorithm for reconstructing a signal from its dyadic wavelet transform modulus maxima

INFINITE SEQUENCES AND SERIES

Chapter 4. Fourier Series

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

Chapter 6 Infinite Series

SCORE. Exam 2. MA 114 Exam 2 Fall 2016

Section A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics

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

Inverse Nodal Problems for Differential Equation on the Half-line

Solutions to Homework 7

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

PAPER : IIT-JAM 2010

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Computation of Error Bounds for P-matrix Linear Complementarity Problems

Math Solutions to homework 6

Riesz-Fischer Sequences and Lower Frame Bounds

MATH 10550, EXAM 3 SOLUTIONS

Physics 116A Solutions to Homework Set #1 Winter Boas, problem Use equation 1.8 to find a fraction describing

Numerical Method for Blasius Equation on an infinite Interval

It is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.

SCORE. Exam 2. MA 114 Exam 2 Fall 2017

The Choquet Integral with Respect to Fuzzy-Valued Set Functions

Newton Homotopy Solution for Nonlinear Equations Using Maple14. Abstract

PRELIM PROBLEM SOLUTIONS

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

On Some Inverse Singular Value Problems with Toeplitz-Related Structure

Weighted Approximation by Videnskii and Lupas Operators

ON MEAN ERGODIC CONVERGENCE IN THE CALKIN ALGEBRAS

Council for Innovative Research

Math 210A Homework 1

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

NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS

A constructive analysis of convex-valued demand correspondence for weakly uniformly rotund and monotonic preference

REAL ANALYSIS II: PROBLEM SET 1 - SOLUTIONS

Sequences and Series of Functions

SOME TRIGONOMETRIC IDENTITIES RELATED TO POWERS OF COSINE AND SINE FUNCTIONS

-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION

Benaissa Bernoussi Université Abdelmalek Essaadi, ENSAT de Tanger, B.P. 416, Tanger, Morocco

f(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim

Concavity Solutions of Second-Order Differential Equations

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

Math 132, Fall 2009 Exam 2: Solutions

Solutions to Final Exam Review Problems

Math 113 Exam 3 Practice

Analytic Continuation

L 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!

MAT 271 Project: Partial Fractions for certain rational functions

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

Alternating Series. 1 n 0 2 n n THEOREM 9.14 Alternating Series Test Let a n > 0. The alternating series. 1 n a n.

Math 113 Exam 3 Practice

Løsningsførslag i 4M

PC5215 Numerical Recipes with Applications - Review Problems

Stability of fractional positive nonlinear systems

The Sumudu transform and its application to fractional differential equations

Singular Continuous Measures by Michael Pejic 5/14/10

On the Lebesgue constant for the Xu interpolation formula

The picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION

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

Complex Analysis Spring 2001 Homework I Solution

Transcription:

Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios with cosecat kerel via Semi-trigoometric iterpolatio Xigtia Gog 1 ad Shuwei Yag 1 LogQiao College of Lazhou Uiversity of Fiace ad Ecoomics Lazhou 730101 Gasu Chia. School of Sciece Lazhou Uiversity of Techology Lazhou 730050 Gasu Chia. Abstract I the paper Semi-trigoometric Lagrage iterpolatio method for solvig sigular itegral equatios (SIE) with cosecat kerel is proposed. For solvig SIE difficulties lie i its sigular term. I order to remove sigular term of SIE we choose appropriate weights quadrature formula. Compared with kow ivestigatios its advatages are that the preset method requires the least computatioal effort ad accuracy i umerical computatio is higher. The fial umerical experimet illustrate the method is efficiet. Keywords: sigular itegral equatios; cosecat kerel; Semi-trigoometric; Lagrage iterpolatio; collocatio method 010 AMS Subect Classificatios: 45E99; 30E0

46 Xigtia Gog ad Shuwei Yag 1. INTRODUCTION This paper is cocered with umerical evaluatio for sigular itegral equatios with cosecat kerel of the form b() t t a( t) ( t) csc ( ) d K( t ) ( ) d f ( t) 0 0 t[0 ) (1) where a b H K H is a give costat to fid the solutio H. So far the study o equatio (1) is oly limited to some discussios o its characteristic equatio [3-6]there are few papers devoted to the umerical method of Equatio (1). Du ad She [1] applied the reproducig kerel method to solve this problem while Wei ad Zhog [] adopted Chebyshev polyomial method. I this paper a ew collocatio method is developed by usig Semi-trigoometric Lagrage iterpolatio polyomial. Lagrage basis for the Semi-trigoometric iterpolatio are firstly give meawhile we give umerical quadrature rules for the improper itegral with cosecat kerel usig it. Through compariso the umerical results show the accuracy ad effectiveess of this method. The work is orgaized as follows. Sectio gives basic kowledge ad a useful lemma for the followig cotet. Sectio 3 is devoted to describig the umerical solutio of the itegral equatio with cosecat kerel. I Sectio 4 the umerical results cofirm that the method is accurate ad efficiet. Sectio 5 eds this paper with a brief coclusio.. SEMI-TRIGONOMETRIC INTERPOLATION Theorem.1 The system 1 1 3 3 1 1 si tcos tsi tcos t si k tcos k t () is orthogoal i the space L ([ ]) moreover it is a maximal orthogoal set Ref.[8]. Defiitio.1 For we deote by T 1 the liear space of Semi-trigoometric polyomial 1 1 q ( t) ak si k t bk cos k t k0 (3)

A collocatio method for sigular itegral equatios with cosecat kerel 47 with real (or complex) coefficiets a 0 a ad b 0 b. A Semi-trigoometric polyomial qt 1 () t is said to be of degree if a b 0. If f() t satisfies f ( t ) f ( t) t the we call it ati-periodic fuctio we kow f () t T 1. Theorem. There exists a uique Semi-trigoometric polyomial 1 1 q ( t) ak si k t bk cos k t k0 (4) satisfyig the iterpolatio property q ( t ) y 0 1. Its coefficiets are give by 1 1 1 k ( )si 0 a f t k t 1 1 1 k ( )cos 0 b f t k t where t 0 1 are equidistat subdivisio of the iterval [0 ] with grid poits of eve umber. From this we deduce that the Lagrage basis for Semi-trigoometric iterpolatio have the form 1 1 L( t) cos k t t k 0 (5) for t [0 ] ad 0 1. Lemma.3 For the Semi-trigoometric polyomial we have itegrals 0 t 1 csc si k tdt 0 t 1 csc cos k tdt 0 i the sese of Cauchy pricipal values.

48 Xigtia Gog ad Shuwei Yag Proof. Usig the followig Lagrage trigoometric idetities ad the itegratig it with respect to t o the iterval[0 ]. The sigular part of the itegral equatio is give by t T0 ( t) : csc ( ) d 0 0 t. (6) The Semi-trigoometric basis fuctios i( k 1 ) t f ( t) : e k are eigefuctios i.e. k T f if 0 k k (7) from the Riesz theory for compact operators we have that T0 has a bouded iverse if ad oly if T0 is iective. We ow costruct umerical quadratures for the improper itegral t Q ( t) : csc ( ) d 0 0t (8) by replacig the cotiuous ati-periodic fuctio H by its semi-trigoometric Lagrage iterpolatio polyomial described i Theorem.1. Usig the Lagrage basis we obtai 1 ( ) ( ) ( ) ( ) 0 Q t C t t (9) with the equidistat quadrature poits t ad the quadrature weights ( ) t C ( t) csc L( ) d 0 1. (10) 0 Usig Lemma.3 ad from the form (4) of the Lagrage basis we derive ( ) 1 C ( t) si k t t k0 0 1. (11)

A collocatio method for sigular itegral equatios with cosecat kerel 49 This quadrature is uiformly coverget for all semi-trigoometric polyomials sice by costructioq itegrates Semi-trigoometric polyomial of degree less tha or equal to exactly. Sice K H the Kt ( ) ( ) is periodic fuctio for variable thus we adopt composite trapezoidal rule of periodic fuctio for Kt ( ) ( ) 1 K( t ) ( ) d K( t t ) ( t ) 0 (1) 0 with the equidistat quadrature poits t 0 1. 3. NUMERICAL SOLUTION For the umerical solutio of the itegral equatios we describe the collocatio method based o Semi-trigoometric iterpolatio that is fully discreted via appropriate quadrature approximatio of itegral operators. I the collocatio method we approximate the itegrals by the quadrature formula (6) ad (9) that is we replace the (SIE) (1) by bt () a( t) ( t) C ( t) ( t ) K( t t ) ( t ) f ( t) 1 1 ( ) 0 0 (13) for t [0 ). Now solvig (10) reduces to solvig a fiite dimesioal liear system. I particular for ay solutio of (10) the values ( t ) i 01 1 trivially satisfy the liear system at the quadrature poits i i bt ( ) a t t C t K t t t f t 1 i ( ) ( i) ( i) ( i) ( i ) ( ) ( i) 0 (14) for i01 1. Ad coversely give a solutio i i01 1of the system (11) the the fuctio iterpolatio polyomial is defied by its Semi-trigoometric Lagrage

50 Xigtia Gog ad Shuwei Yag 1 ( t) : L ( t) 0 t [0 ). (15) 4. NUMERICAL EXAMPLES Let us cosider the SIE with cosecat kerel i Ref.[1]: 1 t t 3 ( t) csc ( ) d 5 ( ) d f ( t) 0 0 (16) where right term f ( t) 10( t ) cos( t ) 3si( t ) ad the exact solutio is ( t) si( t ). I order to imply the covergece of the method we take differet the Table 1 illustrates the covergece behaviour through the differece betwee the exact solutio () t ad the approximate solutio () t at five poits. Of course the approximatio () t is obtaied from via (1) by the correspodig Semi-trigoometric Lagrage iterpolatio. covergece is clearly exhibited. From Table 1 we ote that rapid Table 1. Absolutio error betwee the exact solutio ad approximate solutio. t 10 t 5 t t 4 5 t 9 10 8 3.6311E-3 3.768E-3 3.1554E-3.6393E-3.3314E-3 16 8.5054E-4 9.93E-4 8.4496E-4 6.8569E-4 6.405E-4 64 5.7050E-5 6.119E-5 5.5816E-5 4.548E-5 4.1813E-5 56 3.5870E-6 3.8509E-6 3.5388E-6.8745E-6.6475E-6 51 9.0099E-7 9.6455E-7 8.8684E-7 7.034E-7 6.680E-7 048 5.6334E-8 6.0354E-8 5.557E-8 4.5106E-8 4.1497E-8 For the error ad covergece aalysis of this method we poited out that as the approximate solutio coverges uiformly to the exact solutio of itegral equatio ad the covergece order of quadrature errors for (6) ad (9) carries over to

A collocatio method for sigular itegral equatios with cosecat kerel 51 the error. For detailed descriptio of the quadrature error we refer to Kress [7] ad ivolved literatures therei. 5. CONCLUSION I this paper we give the approximate solutio of SIE with cosecat kerel by Semi-trigoometric iterpolatio polyomial. The fial umerical results show that this method is effective. Moreover this method ca be exteded to solve other equatios with ati-periodic itegral kerel. REFERENCES [1] H. Du ad J. She Reproducig kerel method of solvig sigular itegral equatio with cosecat kerel J. Math. Aaly. Appl. 348 (008) pp. 308 314. [] W. Jiag Z. Che C. Zhag Chebyshev collocatio method for solvig sigular itegral equatio with cosecat kerel Iteratioal Joural of Computer Mathematics 89:8(01) pp. 975-98. [3] D. Jiu ad J. Du Sigular itegral equatios with cosecat kerel i solutios with sigularities of order oe Acta Math. Sci. 6B (006) pp. 410 414. [4] D. JiuO the umerieal solutio for sigular itegral equatios with coseeat kerel.ph.d Thesis Wuha Uiversity005 pp.10-30. [5] J. XiaoPiO paratrigoometric Lagrage iterpolatio ad its applicatio i sigular itegral equatio with cosecat kerel.ph.d Thesis Cetral Chia Normal Uiversity 007 pp.0-43. [6] D. Jiyua H. Huili J. Guoxiag. O trigoometric ad paratrigoometric Hermite iterpolatio. Joural of Approximatio Theory 004 131: 74-99. [7] Kress R. Liear Itegral Equatios 3rd. ed. Spriger New York 014. [8] E. Berriochoa About a system of ati-periodic trigoometric fuctios Computers Mathematics with Applicatios56(008) pp.156-1537.

5 Xigtia Gog ad Shuwei Yag