E408: Conversion of Chebyshev to Power and Power to Chebyshev Series

Author(s): K.S. Kölbig Library: MATHLIB
Submitter: Submitted: 15.02.1994
Language: Fortran Revised:

Subroutine subprograms RCHPWS, RPWCHS and DCHPWS, DPWCHS perform the conversion of a finite Chebyshev series to a finite power series (i.e. a polynomial) and vice versa.

Thus, given the coefficients tex2html_wrap_inline98 , tex2html_wrap_inline100 of a finite Chebyshev series, RCHPWS and DCHWPS calculate the coefficients tex2html_wrap_inline102 , tex2html_wrap_inline104 of the equivalent polynomial:

displaymath106

Conversely, given the coefficients tex2html_wrap_inline108 , tex2html_wrap_inline110 of a power series, RPWCHS and DPWCHS calculate the coefficients tex2html_wrap_inline112 , tex2html_wrap_inline114 of the equivalent finite Chebyshev series:

displaymath116

In both cases, tex2html_wrap_inline118 is the Chebyshev polynomial of degree j.

Note that sometimes the constant term in the Chebyshev series is defined differently, i.e. tex2html_wrap_inline122 instead of tex2html_wrap_inline124 . Here, the definition of Ref. 1 is used.

On computers other than CDC or Cray, only the double-precision versions DCHPWS and DPWCHS are available. On CDC and Cray computers, only the single-precision versions RCHPWS and RPWCHS are available.

Structure:

SUBROUTINE subprograms
User Entry Names: RCHPWS, RPWCHS, DCHPWS, DPWCHS
Files referenced: Unit 6

Usage:

For tex2html_wrap_inline126 (type REAL), tex2html_wrap_inline128 (type DOUBLE PRECISION),

    CALL tCHPWS(N,C,A)
N
(INTEGER) Degree n of last Chebyshev polynomial in the expansion.
C
(type according to t) One-dimensional array of dimension (0:d), where tex2html_wrap_inline132 . On entry, C must contain the coefficients tex2html_wrap_inline134 , tex2html_wrap_inline136 of the Chebyshev expansion.
A
(type according to t) One-dimensional array of dimension (0:d), where tex2html_wrap_inline138 . On exit, A contains the coefficients tex2html_wrap_inline140 , tex2html_wrap_inline142 of the power series expansion.
    CALL tPWCHS(N,A,C)
N
(INTEGER) Degree n of the polynomial.
A
(type according to t) One-dimensional array of dimension (0:d), where tex2html_wrap_inline146 . On entry, A must contain the coefficients tex2html_wrap_inline148 , tex2html_wrap_inline150 of the polynomial.
C
(type according to t) One-dimensional array of dimension (0:d), where tex2html_wrap_inline152 . On exit, C contains the coefficients tex2html_wrap_inline154 , tex2html_wrap_inline156 of the Chebyshev expansion.

Error handling:

Error E408.1: tex2html_wrap_inline158 or tex2html_wrap_inline160 .
A message is written on Unit 6, unless subroutine MTLSET (N002) has been called.

References:

  1. Y.L. Luke, Mathematical functions and their approximations, (Academic Press, New York 1975)
tex2html_wrap_inline162

Michel Goossens Wed Jun 5 03:45:51 METDST 1996