keisoku

Gauss-Jacobi Quadrature Nodes and Weights

Compute the nodes and weights of Gauss-Jacobi quadrature for the weight (1-x)^a (1+x)^b from the degree n and parameters a and b, shown in a table.

Input

Enter the degree n and weight parameters a and b to compute the Gauss-Jacobi quadrature nodes and weights for the weight (1-x)^a (1+x)^b.

Integer from 1 to 128 (the number of nodes)

Exponent of (1-x)^a. Must be greater than minus 1

Exponent of (1+x)^b. Must be greater than minus 1

Result

Degree n

5

Weight (1-x)^a (1+x)^b with a=0, b=0

Number of nodes

5

Sum of weights mu0

2

Weight function

(1-x)^a (1+x)^b

Nodes and weights

Nodes x_i are sorted ascending within the interval from minus 1 to 1, with the matching quadrature weights w_i.

Index iNode x_iWeight w_i
1-0.90617984590.2369268851
2-0.53846931010.4786286705
3-2.2967128660e-170.5688888889
40.53846931010.4786286705
50.90617984590.2369268851

How it works

  • Computes the nodes and weights of degree n Gauss-Jacobi quadrature for the weight (1-x)^a (1+x)^b on the interval from minus 1 to 1.
  • Nodes and weights are obtained by the Golub-Welsch method. The three-term recurrence coefficients of the Jacobi polynomials form a symmetric tridiagonal matrix whose eigenvalues are the nodes and whose first eigenvector components give the weights.
  • The eigendecomposition uses the implicitly shifted QL iteration for symmetric tridiagonal matrices.
  • The zeroth moment mu0 equals the sum of the weights, given by 2^(a+b+1) times Gamma(a+1)Gamma(b+1)/Gamma(a+b+2).
  • Both a and b must be greater than minus 1 so the weight is integrable. When a equals b equals 0 the rule reduces to Gauss-Legendre quadrature.

Reviews

Tell us what you think of this calculator.

Write a review

  1. Home
  2. Gauss-Jacobi Quadrature Nodes and Weights