orthopolynom (version 1.0-6.1)

jacobi.g.weight: Weight function for the Jacobi polynomial

Description

This function returns the value of the weight function for the order \(k\) Jacobi polynomial, \(G_k \left( {p,q,x} \right)\).

Usage

jacobi.g.weight(x,p,q)

Value

The value of the weight function

Arguments

x

the function argument which can be a vector

p

the first polynomial parameter

q

the second polynomial parameter

Author

Frederick Novomestky

Details

The function takes on non-zero values in the interval \( \left( 0,1 \right) \). The formula used to compute the weight function is as follows.

\(w\left( x \right) = \left( {1 - x} \right)^{p - q} \;x^{q - 1}\)

References

Abramowitz, M. and I. A. Stegun, 1968. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover Publications, Inc., New York.

Courant, R., and D. Hilbert, 1989. Methods of Mathematical Physics, John Wiley, New York, NY.

Press, W. H., S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, 1992. Numerical Recipes in C, Cambridge University Press, Cambridge, U.K.

Szego, G., 1939. Orthogonal Polynomials, 23, American Mathematical Society Colloquium Publications, Providence, RI.

Examples

###
### compute the Jacobi G weight function for argument values
### between 0 and 1
### parameter p is 3 and q is 2
###
x <- seq( 0, 1, .01 )
y <- jacobi.g.weight( x, 3, 2 )