^{1}

^{1}

^{1}

^{1}

^{1}

In this paper, we present formulas that turn finite power series into series of shifted Chebyshev polynomials of the first kind. Thereafter, we derive formulas for coefficients of economized power series obtained by truncating the resulting Chebyshev series. To illustrate the utility of our formulas, we apply them to the solution of first order ordinary differential equations via Taylor methods and to solving the Schrödinger equation (SE) for a spherically symmetric hyperbolic potential via the Fröbenius method. In each of the two applications, we show that the use of our formulas makes it possible to reduce the computing time, while preserving the accuracy of the results.

In this paper, we consider the formulas

∑ i = 0 N a i x i = ∑ i = 0 N c i T i * ( x ) (1)

and

∑ i = 0 N − k c i T i * ( x ) = ∑ i = 0 N − k a i [ k ] x i (2)

where the a i coefficients in the truncated power series intervening in the left-hand side of Equation (1) are known, N > 2 , N 2 ≤ k < N , the c i and a i [ k ] coefficients are to be determined. The shifted Chebyshev polynomials T i * ( x ) are defined in the interval [ 0,1 ] by [

T 0 * ( x ) = 1 / 2 , T n * ( x ) = T n ( 2 x − 1 ) , n = 1 , 2 , 3 , ⋯ , (3)

where the first kind Chebyshev polynomials T n ( x ) are given by the recurrence relation

T n + 1 ( x ) = 2 x T n ( x ) − T n − 1 ( x ) , n = 1 , 2 , 3 , ⋯ , (4)

with the starting values T 0 ( x ) = 1 , T 1 ( x ) = x .

Explicit expressions for the five first shifted Chebyshev polynomials are T 0 * ( x ) = 1 / 2 , T 1 * ( x ) = 2 x − 1 , T 2 * ( x ) = 8 x 2 − 8 x + 1 , T 3 * ( x ) = 32 x 3 − 48 x 2 + 18 x − 1 , and T 4 * ( x ) = 128 x 4 − 256 x 3 + 160 x 2 − 32 x + 1 .

The power of a variable x appeared originally purely in algebraic problems [

The power series method is one of the most powerful analytical methods that physicists have for solving linear differential equations. It gives solutions in the form of a power series expansion about any desired point x = x 0 , x being the independent variable. The FM generalizes the power series one. It gives solutions in the form of a power series multiplied by a logarithm term ln | x − x 0 | or a fractional power ( x − x 0 ) r [

We have to emphasize that the economization technique for power series using Chebyshev polynomials is based on the minimax property [

Consequently, it is for considerable interest to find formulas that turn truncated power series of the form ∑ i = 0 N a i x i into economized power series ∑ i = 0 n e i x i of a smaller degree n such that, in the range of interest, the absolute error introduced by the replacement of the original power series by the economized one is less than some acceptable value E [

| ∑ j = 0 N a j x j − ∑ j = 0 n e j x j | < E , x 0 < x < x 1 (5)

where x 0 = − 1 and x 1 = 1 when first kind Chebyshev polynomials are used in the economization procedure, x 0 = 0 and x 1 = 1 in the case of shifted Chebyshev polynomials of the first kind.

The main purpose of this paper is to find general formulas that directly transform any Nth-degree truncated power series defined in the interval [ 0,1 ] into an economized power series of a smaller degree n, using shifted Chebyshev polynomials of the first kind.

The rest of this paper is organized as follows. In Section 2, we establish formulas to transform any finite power series of a certain degree N into a series of first kind shifted Chebyshev polynomials and vice versa. Thereafter, we derive formulas for coefficients of truncated power series associated with the obtained Chebyshev series truncated to the N − k + 1 first terms, where N > 2 and N 2 < k < N . In Section 3, we illustrate the utility of our formulas by means of two applications: first, we apply them to the solution of the Initial Value Problem (IVP) of the form

y ′ ( t ) ≡ d y d t = f ( t , y ( t ) ) , a ≤ t ≤ b , y ( a ) = y a (6)

via Taylor methods [

It is useful and convenient in various applications to be able to express power series explicitly in terms of power of Chebyshev polynomials, and vice versa. In this section, we concentrate on such expressions in the case of shifted Chebyshev polynomials of the first kind.

We begin by relations between powers of x and shifted Chebyshev polynomials, which are useful to reduce the degree of a polynomial in the interval [ 0,1 ] . We found it clarifying to use the formulas established by Ĺopez, Garcia and Caraveo [

[ 1 2 ( 4 x ) 0 1 2 ( 4 x ) 1 1 2 ( 4 x ) 2 ⋮ 1 2 ( 4 x ) N ] = S [ T 0 * ( x ) T 1 * ( x ) T 2 * ( x ) ⋮ T N * ( x ) ] (7)

and

[ T 0 * ( x ) T 1 * ( x ) T 2 * ( x ) ⋮ T N * ( x ) ] = P [ 1 2 ( 4 x ) 0 1 2 ( 4 x ) 1 1 2 ( 4 x ) 2 ⋮ 1 2 ( 4 x ) N ] (8)

where:

・ T 0 * ( x ) = 1 2 and thus x 0 = 2 T 0 * ( x )

・ S = [ S i j ] i , j = 0 N is a ( N + 1 ) × ( N + 1 ) lower triangular matrix such that [

S i i = 1 , i = 0 , 1 , 2 , ⋯ , N , S i j = 0 , j > i , (9)

S i 0 = 2 ( S i − 1 , 0 + S i − 1 , 1 ) S i j = S i − 1 , j − 1 + 2 S i − 1 , j + S i − 1 , j + 1 } i , j = 1 , ⋯ , N − 1 (10)

S N 0 = 2 ( S N − 1 , 0 + S N − 1 , 1 ) , S N j = S N − 1 , j − 1 + 2 S N − 1 , j + S N − 1 , j + 1 , j ∈ [ 1 , N − 1 ] (11)

P = S − 1 (12)

Since S is a lower triangular matrix with ones on the main diagonal, elements of the matrix P defined by Equation (12) can be computed recursively from [

P j j = 1 , P i j = − ∑ k = j i − 1 S i k P k j , i = j + 1 , ⋯ , N (13)

From Equations (7) and (8), we obtain―after expressing the first four powers of x (resp. the first four T i ⋆ ( x ) polynomials) in terms of T i ⋆ ( x ) polynomials (resp. in terms of powers of x), and then generalizing the resulting formulas―the two following expressions:

x i = 2 4 i ∑ j = 0 i S i j T j * ( x ) , i = 0 , 1 , 2 , ⋯ , N (14)

and

T i * ( x ) = 1 2 ∑ j = 0 i P i j 4 j x j , i = 0 , 1 , 2 , ⋯ , N . (15)

Non-vanishing elements of the first nine lines of matrices S and P are given in the two following equations:

S = [ 1 2 1 6 4 1 20 15 6 1 70 56 28 8 1 252 210 120 45 10 1 924 792 495 220 66 12 1 3432 3003 2002 1001 364 91 14 1 12870 11440 8008 4368 1820 560 120 16 1 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋱ ] (16)

and

P = [ 1 − 2 1 2 − 4 1 − 2 9 − 6 1 2 − 16 20 − 8 1 − 2 25 − 50 35 − 10 1 2 − 36 105 − 112 54 − 12 1 − 2 49 − 196 294 − 210 77 − 14 1 2 − 64 336 − 672 660 − 352 104 − 16 1 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋱ ] (17)

It can be easily verified that

S i , i − 1 = 2 i = − P i , i − 1 , i = 1 , 2 , ⋯ , N ; (18)

S i , i − 2 = i 2 + i ( i − 1 ) ; P i , i − 2 = i 2 + i ( i − 3 ) , i = 2 , 3 , ⋯ , N . (19)

In order to transform the finite power series ∑ i = 0 N a i x i into the Chebyshev series of the form ∑ i = 0 N c i T i * , we first replace each power of x by its expression in terms of first kind Chebyshev polynomials, using relations given in Equation (14). Then we collect all the coefficients with the same Chebyshev polynomial T i * ( x ) ; i = 0 , 1 , 2 , ⋯ , N . We find, after some manipulations, that

∑ i = 0 N a i x i = ∑ i = 0 N c i T i * ( x )

with

c i = 2 ∑ j = i N S j i 4 j a j , i = 0 , 1 , 2 , ⋯ , N . (20)

We consider here the problem of transforming the finite Chebyshev series ∑ i = 0 N c i T i * ( x ) into a finite power series of the same degree, i.e., ∑ i = 0 N a i x i To solve this problem, we begin by expressing each polynomial T i * ( x ) , i = 0 , 1 , ⋯ , N , in terms of powers of x, using relations given in Equation (15). Thereafter, we collect all the coeffcients with the same power of x. We obtain the following expression not found explicitely in the literature to the best of our knowledge.

a i = 4 i 2 ∑ j = i N P j i c j , i = 0 , 1 , 2 , ⋯ , N . (21)

As stated earlier, the economization algorithm is lengthy because it has many distinct phases. More precisely, the economization of a power series, such as ∑ i = 0 ∞ a i x i has four basic steps which are as follows when the economization is based on the use of shifted Chebyshev polynomials of the first kind:

Step 1 Truncate the given power series to obtain a polynomial

P N ( x ) = a 0 + a 1 x + ⋯ + a N x N (22)

of degree N.

Step 2 Expand P N ( x ) in a Chebyshev series as follows:

P N ( x ) = c 0 T 0 * ( x ) + c 1 T 1 * ( x ) + ⋯ + c N T N * ( x ) (23)

making use of relations given in Equation (14).

Step 3 Truncate this Chebyshev series to a smaller number of terms by retaining the first n terms, choosing n so that the maximum error given by

| P N ( x ) − M n ( x ) | ≤ E + | c n + 1 | + ⋯ + | c N | (24)

is acceptable, where M n ( x ) designs the resulting small Chebyshev series, i.e.,

M n ( x ) = c 0 T 0 * ( x ) + c 1 T 1 * ( x ) + ⋯ + c n T n * ( x ) , (25)

whereas E represents the absolute error introduced by the replacement of P N ( x ) with M n ( x ) .

Step 4 Replace T j * ( x ) ( j = 0,1, ⋯ , n ) by its polynomial form, which leads to

P N ( x ) ≈ e 0 + e 1 x + ⋯ + e n x n , (26)

using Equation (15).

If necessary in step 1, i.e., when we have an interval [ a , b ] other than [ 0,1 ] , make a transformation of independent variables so that the economization is valid on that interval, by means of the expression

y = x − a b − a (27)

In this case, it is necessary to change variable back to x after step 4, which makes the computing time more considerable.

To reduce the duration of the economization procedure, it is convenient to have a general formula for the coefficients a i [ k ] , i = 0,1,2, ⋯ , N − k , defined by

∑ i = 0 N − k c i T i * ( x ) = ∑ i = 0 N − k a i [ k ] x i , 1 < k < N (28)

where the superscript [k] means that the finite power series of degree N − k whose the a i [ k ] are coefficients is just the Chebyshev form of P N ( x ) (see equations (23)), truncated to the first N − k + 1 terms, i.e., the last k terms of the right hand side of Equation (23) are neglected.

In the rest of this section, we concentrate on establishing a general formula for the a i [ k ] coefficients. For this purpose, we first derive expressions for the a i [ 1 ] , a i [ 2 ] and a i [ 3 ] coefficients successively, in terms of the a i coefficients, i.e., coefficients of the original Nth-degree truncated power series to be economized. Thereafter, we generalize these expressions to establish a formula for a i [ k ] , i = 0 , 1 , ⋯ , N − k ; 1 < k < N .

We start with the expression (1), where the a i coefficients are known and the c i coefficients are given by Equation (20). From Equation (1), we have

∑ i = 0 N − 1 c i T i * ( x ) + c N T N * ( x ) = ∑ i = 0 N − 1 a i x i + a N x N (29)

where

c N = 2 S N N 4 N a N = 2 4 N a N (30)

according to Equation (20) and to the fact that S N N = 1 . Using Equation (15), we obtain

T N * ( x ) = 1 2 ∑ j = 0 N − 1 P N j 4 j x j . (31)

Combining equations (29), (30) and (31), we have:

∑ i = 0 N − 1 c i T i * ( x ) = ∑ i = 0 N − 1 a i x i + a N x N − a N 4 N P N N 4 N x N − a N 4 N ∑ j = 0 N − 1 P N j 4 j x j (32)

= ∑ i = 0 N − 1 a i x i − a N 4 N ∑ j = 0 N − 1 P N j 4 j x j (33)

since P N N = 1 . This can be rewritten as

∑ i = 0 N − 1 c i T i * ( x ) = ∑ i = 0 N − 1 a i [ 1 ] x i (34)

where

a i [ 1 ] = a i − P N i 4 i 4 N a N , i = 0 , 1 , 2 , ⋯ , N − 1. (35)

From Equation (34), we have

∑ i = 0 N − 2 c i T i * ( x ) + c N − 1 T N − 1 * ( x ) = ∑ i = 0 N − 2 a i [ 1 ] x i + a N − 1 [ 1 ] x N − 1 (36)

where

c N − 1 = 2 ∑ j = N − 1 N S j , N − 1 4 j a j = 2 ( S N − 1 , N − 1 4 N − 1 a N − 1 + S N , N − 1 4 N a N ) (37)

and

T N − 1 * ( x ) = 1 2 [ ∑ j = 0 N − 2 P N − 1 , j 4 j x j + P N − 1 , N − 1 4 N − 1 x N − 1 ] . (38)

Combination of Equations (36), (37) and (38) yields

∑ i = 0 N − 2 c i T i * ( x ) = ∑ i = 0 N − 2 a i [ 1 ] x i + a N − 1 [ 1 ] x N − 1 − ( S N − 1 , N − 1 4 N − 1 a N − 1 + S N , N − 1 4 N a N ) × P N − 1 , N − 1 4 N − 1 x N − 1 − ( S N − 1 , N − 1 4 N − 1 a N − 1 + S N , N − 1 4 N a N ) ∑ j = 0 N − 2 P N − 1 , j 4 j x j . (39)

Collecting all coefficients of x N − 1 and taking into account Equation (18) and the expression of a N − 1 [ 1 ] derived from (35), i.e.,

a N − 1 [ 1 ] = a N − 1 − P N , N − 1 4 N − 1 4 N a N , (40)

we find that the term of degree N − 1 in Equation (39) vanishes. Therefore, this equation can be rewritten as

∑ i = 0 N − 2 c i T i * ( x ) = ∑ i = 0 N − 2 a i [ 2 ] x i (41)

with

a i [ 2 ] = a i − P N − 1 , i 4 i 4 N − 1 a N − 1 − [ S N , N − 1 P N − 1 , i 4 i 4 N + S N N P N i 4 i 4 N ] a N . (42)

In this case, we start with Equation (42) which can be rewritten as follows:

∑ i = 0 N − 3 c i T i * ( x ) + c N − 2 T N − 2 * ( x ) = ∑ i = 0 N − 3 a i [ 2 ] x i + a N − 2 [ 2 ] x N − 2 (43)

where

c N − 2 = 2 ∑ j = N − 2 N S j , N − 2 4 j a j (44)

and

T N − 2 * ( x ) = 1 2 ∑ j = 0 N − 2 P N − 2 , j 4 j x j . (45)

Inserting the last two expressions into Equation (43), we get:

∑ i = 0 N − 3 c i T i * ( x ) = ∑ i = 0 N − 3 a i [ 2 ] x i + a N − 2 [ 2 ] x N − 2 − [ S N − 2 , N − 2 4 N − 2 a N − 2 + S N − 1 , N − 2 4 N − 1 a N − 1 + S N , N − 2 4 N a N ] [ ∑ j = 0 N − 3 P N − 2 , j 4 j x j + P N − 2 , N − 2 4 N − 2 x N − 2 ] (46)

The term of degree N − 2 vanishes. This can be shown by inserting in the above equation expressions of S N − 2, N − 2 ; S N − 1, N − 2 ; S N , N − 2 and a N − 2 [ 2 ] obtained from Equations (9), (18), (19) and (42) respectively, not forgetting to replace the two matrix elements P N − 1, N − 2 and P N , N − 2 intervening in a N − 2 [ 2 ] with their expressions derived from (18) and (19). Therefore, Equation (46) becomes

∑ i = 0 N − 3 c i T i * ( x ) = ∑ i = 0 N − 3 a i [ 3 ] x i (47)

where

a i [ 3 ] = a i − P N − 2 , i 4 i 4 N − 2 a N − 2 − P N − 1 , i 4 i + S N − 1 , N − 2 P N − 2 , i 4 i 4 N − 1 a N − 1 − S N , N − 2 P N − 2 , i 4 i + S N , N − 1 P N − 1 , i 4 i + P N i 4 i 4 N a N , i = 0 , ⋯ , N − 3. (48)

A careful analysis of Equations (35), (42) and (48) suggests the following general formula for the a i [ k ] coefficients:

a i [ k ] = a i − ∑ j = 1 k b j i a N − k + j , i = 0 , 1 , ⋯ , N - k (49)

where

b j i = 4 i ∑ q = 0 j − 1 S N − k + j , N − k + j − q P N − k + j − q , i 4 N − k + j (50)

To the best of our knowledge, this formula has not yet been reported in the literature.

For k = 4 , Equations (49) and (50) yield

a i [ 4 ] = a i − 4 i ( P N − 3 , i 4 N − 3 a N − 3 − S N − 2 , N − 3 P N − 3 , i + P N − 2 , i 4 N − 2 a N − 2 ) − 4 i ( S N − 1 , N − 3 P N − 3 , i + S N − 1 , N − 2 P N − 2 , i + P N − 1 , i 4 N − 1 a N − 1 ) − 4 i ( S N , N − 3 P N − 3 , i + S N , N − 2 P N − 2 , i + S N , N − 1 P N − 1 , i + P N i 4 N a N ) . (51)

To illustrate the utility of our formulas, we apply them to two problems via two distinct methods:

・ The Initial Value Problem (IVP) of the form

y ′ ( t ) = d y d t = f ( t , y ( t ) ) , a ≤ t ≤ b , y ( a ) = y a (52)

via the Taylor method of a certain order,

・ the Boundary Value Problem (BVP) constituted by the SE for a spherically symmetric hyperbolic potential given by [

V ( r ) = C − 2 coth 2 ( α r ) + C − 1 coth ( α r ) + ∑ i = 1 3 C i tanh i ( α r ) + ∑ j = 1 K B j sech 2 j ( α r ) , r ∈ ] 0 , + ∞ [ (53)

and the boundary conditions ψ ( 0 ) = 0 and ψ ( ∞ ) = 0 , the method to be used being the FM.

We here consider the IVP of type (52) that possesses a unique solution on some specified interval. The object of Taylor methods is, in this situation, to obtain approximations to the solution of the IVP in question at N particular equally spaced points

t j = a + j h , j ∈ [ 1 , N ] , N = ⌊ b − a h ⌋ , t 0 = a , (54)

approximations to the numbers y ( t 1 ) , y ( t 2 ) , ⋯ , y ( t N ) , rather than to the curve of y ( t ) . Taylor method of order n is, in the case of the IVP (52), defined by [

w 0 = y a , w i + 1 = w i + h T ( n ) ( t i , w i , h ) , i = 0 , 1 , ⋯ , N − 1 (55)

where

T ( n ) ( t i , w i , h ) = f ( t i , w i ) + h 2 f ′ ( t i , w i ) + ⋯ + h n − 1 n ！ f ( n − 1 ) ( t i , w i ) , (56)

w i being the approximation of y i ≡ y ( t i ) . It is remarquable that Taylor methods are explicit and that the explicit Euler’s method is Taylor’s method of order one. Note that Equation (55) becomes

w 0 = y a , w i + 1 = w i + ∑ j = 1 n f ( n − 1 ) n ! ( t i , w i ) h n (57)

if T ( n ) ( t i , w i , h ) is replaced by its expression given by Equation (56).

The existence of the truncated power series of h in Equation (57) immediately suggests that the economization process can be used to improve the accuracy and speed of Taylor’s methods.

y ′ ( t ) = y − t 2 + 1 , y ( 0 ) = 0.5 (58)

using Taylor methods of orders 14 and 20 with two different step sizes, i.e., h = 0.2 and h = 0.9 . The exact solution, which is [

y ( t ) = ( t + 1 ) 2 − 0.5 e t (59)

is given in colunn 3, global errors for the fourteenth- and twentieth-orders Taylor methods in colunns 4 and 5 respectively.

It is obvious that the accuracy of Taylor’s methods decreases with h and increases with the order of the method.

In

h | t i | y i | w i − y i | |
---|---|---|---|---|

N = 14 | N = 20 | |||

0.2 | 0.0 | 0.5000000000000000000000000000000000000000 | 0.0 | 0.0 |

0.2 | 0.8292986209199150830394640026801629148462 | 1.26875901781491456 × 10^{−23} | 2.07119 × 10^{−}^{35} | |

0.4 | 1.2140876511793648410875735235813888596784 | 3.09933152759944941 × 10^{−23} | 5.05952 × 10^{−}^{35} | |

0.6 | 1.6489405998047455125623161659185677433089 | 5.67829811439411006 × 10^{−23} | 9.26957 × 10^{−}^{35} | |

0.8 | 2.1272295357537661977102312343024616214732 | 9.24731863810221013 × 10^{−23} | 1.509584 × 10^{−}^{34} | |

1.0 | 2.6408590857704773823198562643236687511213 | 1.411837561270498104 × 10^{−2}^{2} | 2.304763 × 10^{−}^{34} | |

1.2 | 3.1799415386317262552346162851991778399632 | 2.069306749691896987 × 10^{−2}^{2} | 3.378052 × 10^{−}^{34} | |

1.4 | 3.7324000165776627063879455523856898739162 | 2.948699800145329006 × 10^{−2}^{2} | 4.813622 × 10^{−}^{34} | |

1.6 | 4.2834837878024425981728568217880178717935 | 4.116057221581538750 × 10^{−2}^{2} | 6.719282 × 10^{−}^{34} | |

1.8 | 4.8151762677935269581344880234861373309194 | 5.655784098582873846 × 10^{−2}^{2} | 9.232818 × 10^{−}^{34} | |

2.0 | 5.3054719505346748863847862697124960934098 | 7.675544775075057786 × 10^{−2}^{2} | 1.2529988 × 10^{−}^{33} | |

2.2 | 5.7274932502829395367641114165556679851399 | 1.0312424714524411566 × 10^{−2}^{1} | 1.6834579 × 10^{−}^{33} | |

2.4 | 6.0484118096791991738810301151660989957414 | 1.3740680715970056325 × 10^{−2}^{1} | 2.2431055 × 10^{−}^{33} | |

2.6 | 6.2281309824991548011245873337079413776030 | 1.8181480769274424322 × 10^{−2}^{1} | 2.9680466 × 10^{−}^{33} | |

2.8 | 6.2176766144514750642509919945374922181378 | 2.3915134663568473533 × 10^{−2}^{1} | 3.9040403 × 10^{−}^{33} | |

0.9 | 0.0 | 0.5000000000000000000000000000000000000000 | 0.0 | 0.0 |

0.9 | 2.3801984444215251680999367181987646522891 | 8.339913295594717 × 10^{−}^{14} | 1.116413 × 10^{−2}^{1} | |

1.8 | 4.8151762677935269581344880234861373309194 | 4.10257533772465599 × 10^{−}^{13} | 5.491866 × 10^{−2}^{1} | |

2.7 | 6.2501341375635829440655034902658021096556 | 1.51360605966344928 × 10^{−}^{12} | 2.0261716 × 10^{−2}^{0} | |

3.6 | 2.8608827781610061237026170504081713635556 | 4.96382689788547359 × 10^{−}^{12} | 6.6447708 × 10^{−2}^{0} | |

4.5 | −14.75856565026090677505772837278718042396 | 1.526130510160405894 × 10^{−}^{11} | 2.0429374 × 10^{−}^{19} | |

5.4 | −69.74320810209354351254734005713950522379 | 4.504410420986339807 × 10^{−}^{11} | 6.0297781 × 10^{−}^{19} | |

6.3 | −218.9959550629645165296943338665826849488 | 1.2925572199616312967 × 10^{−}^{10} | 1.73026711 × 10^{−}^{18} | |

7.2 | −602.4753821972089148436757576493594419345 | 3.6333460109324551554 × 10^{−}^{10} | 4.86373756 × 10^{−}^{18} | |

8.1 | −1564.424037641920666544064178264125187657 | 1.0053662796448713144 × 10^{−}^{9} | 1.3458222046 × 10^{−}^{17} | |

9.0 | −3951.541963787692003854998344716379982506 | 2.7475578102963661256 × 10^{−}^{9} | 3.6779872018 × 10^{−}^{17} |

h | t i | y i | w i − y i |
---|---|---|---|

0.2 | 0.0 | 0.5000000000000000000000000000000000000000 | 0.0 |

0.2 | 0.8292986209199150830394640026801629148462 | −2.777268398886 × 10^{−2}^{2} | |

0.4 | 1.2140876511793648410875735235813888596784 | −6.784326565101 × 10^{−2}^{2} | |

0.6 | 1.6489405998047455125623161659185677433089 | −1.2429592768310 × 10^{−2}^{1} | |

0.8 | 2.1272295357537661977102312343024616214732 | −2.0242051853363 × 10^{−2}^{1} | |

1.0 | 2.6408590857704773823198562643236687511213 | −3.0904622455648 × 10^{−2}^{1} | |

1.2 | 3.1799415386317262552346162851991778399632 | −4.5296389328673 × 10^{−2}^{1} | |

1.4 | 3.7324000165776627063879455523856898739162 | −6.4545990670861 × 10^{−2}^{1} | |

1.6 | 4.2834837878024425981728568217880178717935 | −9.0099029752651 × 10^{−2}^{1} | |

1.8 | 4.8151762677935269581344880234861373309194 | −1.23803103878373 × 10^{−2}^{0} | |

2.0 | 5.3054719505346748863847862697124960934098 | −1.68014947273151 × 10^{−2}^{0} | |

2.2 | 5.7274932502829395367641114165556679851399 | −2.25735312012718 × 10^{−2}^{0} | |

2.4 | 6.0484118096791991738810301151660989957414 | −3.00778617497978 × 10^{−2}^{0} | |

2.6 | 6.2281309824991548011245873337079413776030 | −3.97986152424936 × 10^{−2}^{0} | |

2.8 | 6.2176766144514750642509919945374922181378 | −5.23493798456861 × 10^{−2}^{0} | |

0.9 | 0.0 | 0.5000000000000000000000000000000000000000 | 0.0 |

0.9 | 2.3801984444215251680999367181987646522891 | −4.380236415433 × 10^{−2}^{3} | |

1.8 | 4.8151762677935269581344880234861373309194 | −2.1547286230004 × 10^{−2}^{2} | |

2.7 | 6.2501341375635829440655034902658021096556 | −7.9496658372463 × 10^{−2}^{2} | |

3.6 | 2.8608827781610061237026170504081713635556 | −2.60706971012654 × 10^{−2}^{1} | |

4.5 | −14.75856565026090677505772837278718042396 | −8.01544596253786 × 10^{−2}^{1} | |

5.4 | −69.74320810209354351254734005713950522379 | −2.365777899212222 × 10^{−2}^{0} | |

6.3 | −218.9959550629645165296943338665826849488 | −6.788687128076857 × 10^{−2}^{0} | |

7.2 | −602.4753821972089148436757576493594419345 | −1.9082829692444543 × 10^{−}^{19} | |

8.1 | −1564.424037641920666544064178264125187657 | −5.2803210691366666 × 10^{−}^{19} | |

9.0 | −3951.541963787692003854998344716379982506 | −1.44305490328401511 × 10^{−}^{18} |

Comparing

The time-independent SE for a particule of mass M that moves in three dimemsional space under the effect of a spherically symmetric potential is given by

− ℏ 2 2 M [ Δ r + V ( r ) ] ψ ( r → ) = E ψ ( r → ) (60)

Taking ψ ( r → ) = r − 1 ψ ( r ) Y l m ( θ , φ ) and considering potential (53), it is found that ψ ( r ) must satisfy the radial SE

where

To reduce Equation (61) into a differential equation solvable by means of the FM, we appeal to the following approximation for the centrifugal term 1 / R 2 [

1 R 2 ≈ α ˜ 2 [ s sinh 2 ( α ˜ R ) + t cosh 2 ( α ˜ R ) ] (63)

where α ˜ , s and t are adjustable dimensionless parameters. Substituting this approximation into (61) and making the change of variable ξ = − tanh ( α ˜ R ) , such that the domain 0 < r < ∞ maps to − 1 < ξ < 0 , we find:

ξ 2 ( ξ 2 − 1 ) d 2 ψ d ξ 2 + 2 ξ 3 ( ξ 2 − 1 ) d ψ d ξ + 2 M ℏ 2 α 2 [ − ( C − 2 + l ( l + 1 ) α 2 s 2 M ) + C − 1 ξ + ( E + l ( l + 1 ) α 2 s 2 M ) ξ 2 + C 1 ξ 3 − C 2 ξ 4 + C 3 ξ 5 − ∑ j = 1 K B ˜ j ( 1 − ξ 2 ) j ξ 2 ] ψ ( ξ ) = 0 (64)

where the B ˜ j coefficients are such that

B ˜ 1 = l ( l + 1 ) α 2 t 2 M + B 1 , B ˜ i = B i , i = 2 , 3 , ⋯ , K . (65)

If we make the substitution η = ξ + 1 in Equation (64), such that the domain − 1 < ξ < 0 maps to 0 < η < 1 , and then apply the FM to the resulting differential equation by writing

ψ ( η ) = η δ ∑ n = 0 ∞ a n η n , (66)

we find that the recursion relation between successive coefficients is:

[ 4 + d 0 + ε + 4 δ ( 2 + δ ) + 8 n ( 1 + δ ) + 4 n 2 ] a n + 1 + [ d 1 − 2 ε − 2 δ ( 1 + 6 δ ) − 2 n ( 1 + 12 δ ) − 12 n 2 ] a n + [ 8 + d 2 + ε − δ ( 21 − 13 δ ) − n ( 21 − 26 δ ) + 13 n 2 ] a n − 1 + [ − 16 + d 3 + 2 δ ( 10 − 3 δ ) + 4 n ( 5 − 3 δ ) − 6 n 2 ] a n − 2 + [ 6 + d 4 − δ ( 5 − δ ) − n ( 5 − 2 δ ) + n 2 ] a n − 3 + d 5 a n − 4 + ∑ p = 6 2 K + 2 b p a n − p + 1 = 0 (67)

where

ε = 2 M ℏ 2 α 2 E , d 0 = − 2 M ℏ 2 α 2 [ C − 2 + C − 1 + C 1 + C 2 + C 3 ] (68)

d 1 = 2 M ℏ 2 α 2 [ − 2 l ( l + 1 ) α 2 s 2 M + C − 1 + 3 C 1 + 4 C 2 + 5 C 3 − β 1 ] (69)

d 2 = 2 M ℏ 2 α 2 [ 2 l ( l + 1 ) α 2 s 2 M − 3 C 1 − 6 C 2 − 10 C 3 − β 2 ] (70)

d 3 = 2 M ℏ 2 α 2 ( C 1 + 4 C 2 + 10 C 3 + β 3 ) , d 4 = − 2 M ℏ 2 α 2 ( C 2 + 5 C 3 + β 4 ) (71)

d 5 = 2 M ℏ 2 α 2 ( C 3 − β 5 ) ; b j = − 2 M ℏ 2 α 2 β j , j = 6 , 7 , ⋯ , 2 K + 2. (72)

Here, d i refers to the coefficient of η i in the expression between brackets which intervene in the left-hand side of Equation (64).

It is worth noting that the β j quantities that appear in Equations (69)-(72) are defined as follows:

β 1 = B ^ 1 ; β 2 = − 2 B ^ 1 + B ^ 2 ; β j = B ^ j − 2 − 2 B ^ j − 1 + B ^ j ; j = 3 , 4 , ⋯ , 2 K ; (73)

β 2 K + 1 = B ^ 2 K − 1 − 2 B ^ 2 K ; β 2 K + 2 = B ^ 2 K (74)

where the B ^ j quantities are calculated from the two following equations :

B ^ 2 p + 1 = ∑ j = p + 1 j ≤ K 2 p + 1 ( − 1 ) j + 1 ( j 2 p + 1 − j ) 2 2 j − 2 p − 1 B ˜ j ; p = 0 , ⋯ , K − 1 ; (75)

B ^ 2 p = ∑ j = p j ≤ K 2 p ( − 1 ) j ( j 2 p − j ) 2 2 j − 2 p B ˜ j ; p = 1 , 2 , ⋯ , K . (76)

Let us mention in passing that the potential (53) has been studied in [

Setting n = − 1 in the above recursion relation, we obtain the indicial equation ( 4 δ 2 + d 0 + ε ) a 0 = 0 which is solved by δ 1 = − − d 0 − ε / 2 and δ 2 = − d 0 − ε / 2 . We can therefore conclude that the two solutions of the radial SE obtained as generalized series, one with δ = δ 1 and the other with δ = δ 2 , are linearly independent. Clearly, only δ > 0 is acceptable for bound states since in this case ψ ( r = ∞ ) = ψ ( η = 0 ) = 0 . Such a solution to the SE contains only the series with δ = δ 2 , and in the following will be denoted by

ψ ( η , ε ) = η δ 2 ∑ n = 0 ∞ a n η n (77)

where the dependence on the energy eigenvalue is explicitly marked. For arbitrary ε , the wave function ψ ( η , ε ) obtained from (77) is not square integrable. If, for a particular value of ε , ψ ( η = 0 , ε ) = 0 , the wave function is square integrable. In other words, this condition determines the energies of the discrete spectrum. Substituting the calculated bound state energy into the recursion relation (67), the coefficients of the generalized series (77) can be successively determined in order to obtain the wave function as a sum of the series.

We have to emphasize that in practice, the function ψ ( η , ε ) has to be approximated by truncating the series in Equation (77) at a suitably high order N, which requires high precision computing since some of the elements of matrices S and P needed to compute the a i [ k ] coefficients (se relations (49) and (50)) are extremely big for high values of i. To deal with this situation, we appeal to a software package for arbitrary precision computation, named MPFUN2015 and developed by David H. Bailey [^{−9}. The number of terms in the series (77) is 602, which means that N = 601. One has to add that the truncated series is a function of ε whose zeros correspond to bound state values of ε. If for two neighboring values of ε the wave function at η = 0 ( r → ∞ ) takes values of different signs, we can deduce that one bound state reduced energy (ε) lies between these two values. The bound state reduced energy in question can be computed numerically by means of the dichotomy method [

and which is a special case of the potential (53) since this one reduces to MRP by choosing K = 1 , C − 2 = a ( a − 1 ) ℏ 2 / ( 8 M b 2 ) , C − 1 = − [ a ( a − 1 ) + A ] ℏ 2 / ( 4 M b 2 ) , C 1 = 0 , C 2 = B 1 = [ a ( a − 1 ) + 2 A ] ℏ 2 / ( 8 M b 2 ) and α = 1 / ( 2 b ) .

V eff ( r , l ) = l ( l + 1 ) ℏ 2 2 M r 2 + V ( r ) , (79)

where the dependence on l is explicitly marked, for the first four values of l and given values of a , b and A.

In our calculations, l = 1 , the value of a is 0.75, that of 1/b is 0.025, A = 2 b and the quantum number n ≤ 8 . The parameter b has been chosen such that α ˜ = 0.0001 , which implies that _{21} = −0.12054171328 a.u., E_{31} = −0.04588384531 a.u., E_{41} = −0.02081525239 a.u., E_{51} = −0.00981695557 a.u., E_{61} = −0.00437316621 a.u., E_{71} = −0.00161118952 a.u. and E_{81} = −0.00034181936 a.u.

In order to show that our formulas can help to reduce the computation time, we consider the domain 0 < r < 1000 and then compute the wave functions at r j = 0 + j r s t e p for j = 0,1,2, ⋯ , J where the step size r s t e p = 0.001 a . u . and J is the number of the mesh points minus 1. In our case, J = 100 , 000 , which means that for each bound state, it is necessary to evaluate the sum of the first N + 1 terms of the series intervening in the right-hand side of Equation (77) 100,001 times. Using a (N-k) th-degree economized power series allows to reduce the number of terms in this (N + 1)-term truncated series, which saves us a little time during the computation of the sum. In our calculations, we have chosen k = 200 . We found that the total duration of the whole computation is 9287 seconds (s) when the economization procedure is used to calculate the wave function and 11,540 s in the absence of the economization. The use of the economization procedure thus helps us save 253 s due to the fact that the time spent on computing the a i [ k ] ( i = 0 , 1 , ⋯ , N − k ) coefficients (see Equations (49) and (50)) is compensated by the gain of time when calculating the (N + 1)-term sum mentioned above.

New transformation formulas between finite power series and series of first-kind shifted Chebyshev polynomials are proposed. These formulas have been used to establish formulas that directly turn any Nth-degree truncated power series into an economized one whose degree n is lower than N. We also have shown that the application of these formulas to the Fröbenius and Taylor’s methods for ordinary differential equations can, in some situations, help to reduce the computing time while preserving and even improving the accuracy of the results.

The author declares no conflicts of interest regarding the publication of this paper.

Nyengeri, H., Nizigiyimana, R., Mutankana, J.-P., Bayaga, H. and Bayubahe, F. (2021) Power and Chebyshev Series Transformation Formulas with Applications to Solving Ordinary Differential Equations via the Fröbenius and Taylor’s methods. Open Access Library Journal, 8: e7142. https://doi.org/10.4236/oalib.1107142