
Calling Sequence


chebpade(f, x=a..b, [m, n])
chebpade(f, x, [m, n])
chebpade(f, a..b, [m, n])


Parameters


f



procedure or expression representing the function to be approximated

x



the variable appearing in f, if f is an expression

a, b



numerical values specifying the interval of approximation

m, n



desired degree of numerator and denominator, respectively





Description


•

The function chebpade computes a ChebyshevPade approximation of degree $m,n$ for the function f.

•

Specifically, f is expanded in a Chebyshev series on the interval $a..b$ (if $a..b$ is not specified then the interval $\mathrm{1}..1$ is understood), and then the ChebyshevPade rational approximation is computed.

•

If $n=0$ or if the third argument is simply an integer m then the Chebyshev series of degree m is computed.

•

The $m,n$ ChebyshevPade approximation is defined to be the rational function $\frac{p\left(x\right)}{q\left(x\right)}$ with $\mathrm{deg}\left(p\left(x\right)\right)\le m$ and $\mathrm{deg}\left(q\left(x\right)\right)\le n$ such that the Chebyshev series expansion of $\frac{p\left(x\right)}{q\left(x\right)}$ has maximal initial agreement with the Chebyshev series expansion of f. In normal cases, the series expansion agrees through the term of degree $m+n$.

•

If the second argument is a range $a..b$ then the first argument is understood to be a Maple operator, and the result will be returned as an operator. If the second argument is an equation $x=a..b$, or a name x, then the first argument is understood to be an expression in the variable x, and the result will be returned as an expression. In all cases, the numerator and denominator will be expressed in terms of the Chebyshev polynomials $T\left(n\,x\right)$. See orthopoly[T].

•

The method used is based on transforming the Chebyshev series to a power series with the same coefficients, computing a Pade approximation for the power series, and then converting back to the appropriate ChebyshevPade approximation.

•

Note that for the purpose of evaluating a rational function efficiently (i.e. minimizing the number of arithmetic operations), the rational function should be converted to a continuedfraction form. See numapprox[confracform].

•

Various levels of user information will be displayed during the computation if infolevel[chebpade] is assigned values between 1 and 3.

•

The command with(numapprox,chebpade) allows the use of the abbreviated form of this command.



Examples


>

$\mathrm{with}\left(\mathrm{numapprox}\right)\:$

>

$\mathrm{chebpade}\left(\mathrm{exp}\left(x\right)\,x=0..1\,5\right)$

${1.75338765437709}{}{T}{}\left({0}{\,}{2}{}{x}{}{1}\right){+}{0.850391653780811}{}{T}{}\left({1}{\,}{2}{}{x}{}{1}\right){+}{0.105208693630937}{}{T}{}\left({2}{\,}{2}{}{x}{}{1}\right){+}{0.00872210473331555}{}{T}{}\left({3}{\,}{2}{}{x}{}{1}\right){+}{0.000543436831150254}{}{T}{}\left({4}{\,}{2}{}{x}{}{1}\right){+}{0.0000271154349130735}{}{T}{}\left({5}{\,}{2}{}{x}{}{1}\right)$
 (1) 
>

$\mathrm{chebpade}\left(\frac{\mathrm{sin}\left(x\right)}{x}\,x=0..2\,\left[2\,2\right]\right)$

$\frac{{0.771073733750623}{}{T}{}\left({0}{\,}{x}{}{1}\right){}{0.221091073962959}{}{T}{}\left({1}{\,}{x}{}{1}\right){}{0.0421244668861024}{}{T}{}\left({2}{\,}{x}{}{1}\right)}{{T}{}\left({0}{\,}{x}{}{1}\right){+}{0.0836360586596837}{}{T}{}\left({1}{\,}{x}{}{1}\right){+}{0.0336007994536882}{}{T}{}\left({2}{\,}{x}{}{1}\right)}$
 (2) 
>

$\mathrm{chebpade}\left(\mathrm{sin}+\mathrm{cos}\,1..1\,\left[1\,2\right]\right)$

${x}{\mapsto}\frac{{0.865073720371733357}{\cdot}{T}{}\left({0}{\,}{x}\right){+}{1.11718143936440173}{\cdot}{T}{}\left({1}{\,}{x}\right)}{{T}{}\left({0}{\,}{x}\right){+}{0.266549197507840174}{\cdot}{T}{}\left({1}{\,}{x}\right){+}{0.151597635961381677}{\cdot}{T}{}\left({2}{\,}{x}\right)}$
 (3) 
>

$r\u2254\mathrm{confracform}\left(\mathrm{subs}\left(T=\mathrm{orthopoly}\left[T\right]\,\right)\right)$

${r}{\u2254}{y}{\mapsto}\frac{{3.68469281291759376}}{{y}{+}{0.104797831264514252}{+}\frac{{2.71705574090406143}}{{y}{+}{0.774335922429842616}}}$
 (4) 
Check agreement of the Chebyshev series to 7 digits of accuracy.
>

$\mathrm{Digits}\u22547$

${\mathrm{Digits}}{\u2254}{7}$
 (5) 
>

$\mathrm{chebpade}\left(r\,1..1\,4\right)$

${x}{\mapsto}{0.7651975}{\cdot}{T}{}\left({0}{\,}{x}\right){+}{0.8801012}{\cdot}{T}{}\left({1}{\,}{x}\right){}{0.2298070}{\cdot}{T}{}\left({2}{\,}{x}\right){}{0.03912671}{\cdot}{T}{}\left({3}{\,}{x}\right){+}{0.02274805}{\cdot}{T}{}\left({4}{\,}{x}\right)$
 (6) 
>

$\mathrm{chebpade}\left(\mathrm{sin}+\mathrm{cos}\,1..1\,4\right)$

${x}{\mapsto}{0.7651975}{\cdot}{T}{}\left({0}{\,}{x}\right){+}{0.8801012}{\cdot}{T}{}\left({1}{\,}{x}\right){}{0.2298070}{\cdot}{T}{}\left({2}{\,}{x}\right){}{0.03912671}{\cdot}{T}{}\left({3}{\,}{x}\right){+}{0.004953278}{\cdot}{T}{}\left({4}{\,}{x}\right)$
 (7) 


References



Geddes, K.O. "Block Structure in the ChebyshevPade Table." SIAM J. Numer. Anal., Vol. 18(5). (Oct. 1981): 844861.



