File Exchange

## Gaussian quadratures for several orthogonal polynomials

version 1.1.0.0 (4.87 KB) by
This function calculates the zeros and weights of several orthogonal polynomials

2 Downloads

Updated 17 Oct 2014

View License

The function calculates the zeros and weights of several orthogonal polynomials to be used in particular numerical integration problems. The quadrature rules implemented are the Hermite (probabilist-type), Hermite (physicist-type), Legendre, Chebyshev and Laguerre.
An interesting contribution is the (probabilist-type) Gauss-Hermite quadrature, which is validated through an example by comparing the results of the numerical integration with the moments of a standard Gaussian variable (see 'examples' section). Furthermore, the function displays two figures, the first shows roots vs. weights, and the second shows the corresponding orthogonal polynomials up to the specified order m.

Finally, it can be seen that other orthogonal polynomials can be easily included in the function (case ...) due to the general implementation of the weight's formula.
-----------------------------------------------------------
1. Input: * m - number of quadrature points
* type - orthogonal polynomial:
'he_prob': Hermite probabilist
'he_phys': Hermite physicist
'legen' : Legendre
'cheby' : Chebyshev
'lague' : Laguerre
2. Output: * xi - zeros
* w - weights
-----------------------------------------------------------

### Cite As

Felipe Uribe (2021). Gaussian quadratures for several orthogonal polynomials (https://www.mathworks.com/matlabcentral/fileexchange/48144-gaussian-quadratures-for-several-orthogonal-polynomials), MATLAB Central File Exchange. Retrieved .

### Comments and Ratings (2)

haisong zhao

fadoua sb

Error using roots (line 27)
Input to ROOTS must not contain NaN
or Inf.

Error in gauss_quad (line 116)
xi = sort(roots(P{m+1}));

##### MATLAB Release Compatibility
Created with R2013a
Compatible with any release
##### Platform Compatibility
Windows macOS Linux

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!