Documentation Center

  • Trial Software
  • Product Updates

orthpoly::curtz

The Curtz polynomials

Use only in the MuPAD Notebook Interface.

This functionality does not run in MATLAB.

Syntax

orthpoly::curtz(n, x)

Description

orthpoly::curtz(n,x) computes the value of the n-th degree Curtz polynomial at the point x.

These polynomials have rational coefficients.

Evaluation for real floating-point values x is numerically stable. Cf. Example 2.

Examples

Example 1

Polynomials of domain type DOM_POLY are returned, if identifiers or indexed identifiers are specified:

orthpoly::curtz(2, x)

orthpoly::curtz(3, x[1])

However, using arithmetical expressions as input the "values" of these polynomials are returned:

orthpoly::curtz(2, 6*x)

orthpoly::curtz(3, x[1] + 2)

"Arithmetical expressions" include numbers:

orthpoly::curtz(2, sqrt(2)), orthpoly::curtz(3, 8 + I),
orthpoly::curtz(100, 0.3)

If no integer degree is specified, then orthpoly::curtz returns itself symbolically:

orthpoly::curtz(n, x), orthpoly::curtz(1/2, x)

Example 2

If a floating-point value is desired, then a direct call such as

orthpoly::curtz(50, 1.2)

is appropriate and yields a correct result. One should not evaluate the symbolic polynomial at a floating-point value, because this may be numerically unstable:

orthpoly::curtz(50, x): evalp(%, x = 1.2)

Note that only 3 digits are correct due to numerical round-off.

Parameters

n

A nonnegative integer: the degree of the polynomial.

x

An indeterminate or an arithmetical expression. An indeterminate is either an identifier (of domain type DOM_IDENT) or an indexed identifier (of type "_index").

Return Values

If x is an indeterminate, then a polynomial of domain type DOM_POLY is returned. If x is an arithmetical expression, then the value of the Curtz polynomial at this point is returned as an arithmetical expression. If n is not a nonnegative integer, then orthpoly::curtz returns itself symbolically.

Algorithms

The Curtz polynomials are given by the recursion formula

with C(0, x) = 1.

Was this topic helpful?