scholarly journals Some orthogonal polynomials on the finite interval and Gaussian quadrature rules for fractional Riemann‐Liouville integrals

2020 ◽  
Vol 44 (1) ◽  
pp. 493-516
Author(s):  
Gradimir V. Milovanović
Author(s):  
Walter Gautschi

The connection between orthogonal polynomials and quadrature rules has already been elucidated in §1.4. The centerpieces were the Gaussian quadrature rule and its close relatives—the Gauss–Radau and Gauss–Lobatto rules (§1.4.2). There are, however, a number of further extensions of Gauss’s approach to numerical quadrature. Principal among them are Kronrod’s idea of extending an n-point Gauss rule to a (2n + 1)-point rule by inserting n + 1 additional nodes and choosing all weights in such a way as to maximize the degree of exactness (cf. Definition 1.44), and Turán’s extension of the Gauss quadrature rule allowing not only function values, but also derivative values, to appear in the quadrature sum. More recent extensions relate to the concept of accuracy, requiring exactness not only for polynomials of a certain degree, but also for rational functions with prescribed poles. Gauss quadrature can also be adapted to deal with Cauchy principal value integrals, and there are other applications of Gauss’s ideas, for example, in combination with Stieltjes’s procedure or the modified Chebyshev algorithm, to generate polynomials orthogonal on several intervals, or, in comnbination with Lanczos’s algorithm, to estimate matrix functionals. The present section is to discuss these questions in turn, with computational aspects foremost in our mind. We have previously seen in Chapter 2 how Gauss quadrature rules can be effectively employed in the context of computational methods; for example, in computing the absolute and relative condition numbers of moment maps (§2.1.5), or as a means of discretizing measures in the multiple-component discretization method for orthogonal polynomials (§2.2.5) and in Stieltjes-type methods for Sobolev orthogonal polynomials (§2.5.2). It is time now to discuss the actual computation of these, and related, quadrature rules.


2011 ◽  
Vol 2011 ◽  
pp. 1-21
Author(s):  
Samir A. Ashour ◽  
Hany M. Ahmed

Many algorithms that have been proposed for the numerical evaluation of Cauchy principal value integrals are numerically unstable. In this work we present some formulae to evaluate the known Gaussian quadrature rules for finite part integrals , and extend Clenshow's algorithm to evaluate these integrals in a stable way.


1970 ◽  
Vol 29 ◽  
pp. 117-125
Author(s):  
HT Rathod ◽  
RD Sathish ◽  
Md Shafiqul Islam ◽  
Arun Kumar Gali

Gauss Legendre Quadrature rules are extremely accurate and they should be considered seriously when many integrals of similar nature are to be evaluated. This paper is concerned with the derivation and computation of numerical integration rules for the three integrals: (See text for formulae) which are dependent on the zeros and the squares of the zeros of Legendre Polynomial and is quite well known in the Gaussian Quadrature theory. We have developed the necessary MATLAB programs based on symbolic maths which can compute the sampling points and the weight coefficients and are reported here upto 32 – digits accuracy and we believe that they are reported to this accuracy for the first time. The MATLAB programs appended here are based on symbolic maths. They are very sophisticated and they can compute Quadrature rules of high order, whereas one of the recent MATLAB program appearing in reference [21] can compute Gauss Legendre Quadrature rules upto order twenty, because the zeros of Legendre polynomials cannot be computed to desired accuracy by MATLAB routine roots (……..). Whereas we have used the MATLAB routine solve (……..) to find zeros of polynomials which is very efficient. This is worth noting in the present context. GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 29 (2009) 117-125  DOI: http://dx.doi.org/10.3329/ganit.v29i0.8521 


Sign in / Sign up

Export Citation Format

Share Document