Growth and Form

Volume 1, Issue 1, 2020, Pages 20 - 32

Chebyshev Polynomials, Rhodonea Curves and Pseudo-Chebyshev Functions. A Survey

Paolo Emilio Ricci*
Department of Mathematics, International Telematic University UniNettuno, Corso Vittorio Emanuele II, 39, Roma 00186, Italy
Corresponding Author
Paolo Emilio Ricci
Received 17 December 2019, Accepted 20 December 2019, Available Online 23 February 2020.
10.2991/gaf.k.200124.005How to use a DOI?
Spirals; Grandi curves; pseudo-Chebyshev functions; recurrence relations; differential equations; orthogonality properties

In recent works, starting from the complex Bernoulli spiral and the Grandi roses, sets of irrational functions have been introduced and studied that extend to the fractional degree the polynomials of Chebyshev of the first, second, third and fourth kind. The functions thus obtained are therefore called pseudo-Chebyshev. This article presents a review of the elementary properties of these functions, with the aim of making the topic accessible to a wider audience of readers. The subject is presented as follows. In Section 2 a review of spiral curves is given. In Section 3 the main properties of the classical Chebyshev polynomials are recalled. The Grandi (Rhodonea) curves and possible extensions are introduced in Section 4, and a method for deriving new curves, changing cartesian into polar coordinates, is touched on. The possibility to consider the Grandi curves even for rational indexes allows to introduce in Section 5 the pseudo-Chebyshev functions, which are derived from the Chebyshev polynomials assuming rational values for their degree. The main properties of these functions are shown, including recursions and differential equations. In particular, the case of half-integer degree is examined in Section 6 since, in this case, the pseudo-Chebyshev functions verify even the orthogonality property. As a consequence, new system of irrational orthogonal functions are introduced.

© 2020 The Authors. Published by Atlantis Press SARL
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (


In writing an article on Growth and Shape, one cannot help but link the treatment to geometrical entities that translate those concepts into mathematical terms. They are the logarithmic spiral of Bernoulli, the curves of Lamé, the roses of Grandi, the lemniscate of Bernoulli and their generalizations.

The spiral has always been associated with growth phenomena, starting with that of the shell Nautilus widely studied in the book by Thompson [1] and in many subsequent works.

The Lamé’s curves have been generalized by J. Gielis in the 2D and 3D case in works [2,3] that have had wide international resonance [47].Grandi’s roses (also called Rhodoneas) and Bernoulli’s lemniscate have polar equations that lend themselves to being generalized, as is done here in Section 4.1. All these curves (or surfaces) of the plane (of space) lend themselves to creating mathematical forms that model natural forms [3].

In this article, starting from the spiral of Bernoulli, in the complex form, we make the obvious connection with the first and second kind Chebyshev polynomials, and with the roses of Grandi.

After that, having observed that roses also exist for rational index values, extensions of that polynomials are introduced in the case of fractional degree. Thus, irrational functions are found which are called pseudo-Chebyshev of first and second kind, because they continue to verify many of the properties of the corresponding Chebyshev polynomials.

Subsequently, using the links with Chebyshev polynomials of third and fourth kind and the good work [8], the pseudo-Chebyshev functions of third and fourth kind are also introduced and studied.

Particular importance is given, in Section 6, to the case of the half-integer degree, because, in this case, the pseudo-Chebyshev functions verify not only the corresponding recurrence relations and differential equations, but also the orthogonality properties, in the interval [−1, 1], with respect to the same weights of the classical polynomials.

In this survey we limited ourselves to considering only the most elementary properties of the pseudo-Chebyshev functions, which can be proven starting from trigonometric identities, that are known to secondary school students, so as to make the treatment usable to a wide audience. Moreover, the use of higher tools seems to be unessential in the context of this study, which deals with functions of elementary nature, connected in a simple way to trigonometric functions.


The spiral symbol is found in every ancient culture, all over the world (see e.g. Figures 1, 2). The spiral is a sacred symbol, possibly reminding us the evolution of our life.

Figure 1

Ancient Crete island vases.

Figure 2

A a well of Nazca culture.

The first attempt to describe a spiral is due to Theodore of Cyrene, a mathematician from the school of Pythagoras, in the 5th century bc.

By the mathematical point of view spirals are described by polar equations.

Many information on this subject can be found in Lockwood [9] and in Thompson [1], where applications to natural shapes (see e.g. Figure 3) are deeply analyzed. In a recent article [10] a Bernoulli spiral in complex form has been related to the Grandi (Rhodonea) curves and Chebyshev polynomials.

Figure 3

Spirals - natural shapes [29]

Connection with curvature can be found in Gielis et al. [11].

2.1. Archimedes vs Bernoulli Spiral

The Archimedes (Figure 4) spiral [12] (Figure 5) has the polar equation:

ρ=aθ,(a>0,θIR+). (1)

Figure 4

Archimedes (traditional) and his death by N. Barabino.

Figure 5

Archimedes vs Bernoulli spiral.

If θ > 0 the spiral turns counter-clockwise, if θ < 0 the spiral turns clockwise. Bernoulli’s (logarithmic) spiral [13] (Figure 5) has the polar equation

ρ=abθ,(a,b,IR+),θ=log(ρa). (2)

Varying the parameters a and b one gets different types of spirals.

The size of the spiral depends on a, while the term b controls the verse of rotation and how it is “narrow”.

Being a and b positive costants, there are some interesting cases. The most popular logarithmic spiral is the harmonic spiral, in which the distance between the spires is in harmonic progression, with ratio ϕ=512 , that is the “Golden ratio” of the unit segment.

The logarithmic spiral was discovered by René Descartes in 1638, and studied by Jakob Bernoulli (1654–1705) (Figure 6).

Figure 6

René Descartes by F. Hals and Jakob Bernoulli.

Pierre Varignon (1654–1722) called it spiral equiangular, because:

  1. 1.

    There is a constant angle between the tangent at a given point and the polar radius passing through the same point.

  2. 2.

    The inclination angle with respect to concentric circles centered at the origin is also constant.

It is a first example of a fractal. As it is written on J. Bernoulli’s tomb: Eadem Mutata Resurgo (but the spiral represented there is of Archimedes type).

2.2. Fermat Spiral, Fibonacci and Other Types of Spirals

The Fermat (parabolic) spiral (Figure 7) has polar equation:

ρ=±aθ1/2. (3)

Fermat’s spiral suggests the possibility of introducing other kind of spiral graphs.

In fact there is a straightforward correspondence, between cartesian and polar systems of coordinates, which transforms y = f(x) functions of the (x, y) plane into polar curves ρ = f(θ) of the (ρ, θ) plane.

In this planar transformation, the Archimedes spiral ρ = corresponds to the straight line y = ax, the Bernoulli spiral ρ = abθ to the exponential function y = abx, and the Fermat spiral to the parabolic function y=ax .

Then, putting:

ρ=aθm/n,(m,npositiveintegers,n0), (4)
one gets a parametric family of spirals, at varying m and n.

Notice that, if m > n, so that the exponent is >1, the coils of spiral are widening (Figure 7), while if m < n being the exponent <1, the coils of spiral are shrinking (as in Fermat’s case).

Other possibilities are:

  1. 1.

    To assume θm/n with m/n < 0; in this case the coils are wrapped around the origin.

  2. 2.

    To use a graph with horizontal asymptotes, in order to get an asymptotic spiral (Figure 8).

In what follows, we consider a “canonical form” of the Bernoulli spirals assuming a = 1, b = en, that is, the simplified polar equation:

ρ=enθ,(nN). (5)

2.3. The Complex Bernoulli Spiral

We now introduce the complex case, putting

ρ=ρ+iρ, (6)
and considering a Bernoulli spiral of the type:
ρ=einθ=cosnθ+isinnθ. (7)

Therefore, we have:

ρ1=ρ=cosnθ,ρ2=ρ=sinnθ. (8)

The curves with polar equation:

ρ=cos(nθ) (9)
are known as Rhodonea or Grandi curves, in honour of G. G. Grandi (1671–1742) (Figure 12), who communicated his discovery to G. W. Leibniz (1646–1716), in 1713.

Curves with polar equation: ρ = sin() are equivalent to the preceding ones, up to a rotation of π/(2n) radians.

The Grandi roses display

  • n petals, if n is odd.

  • 2n petals, if n is even.

By using Eq. (9) it is impossible to obtain, roses with 4n + 2 (n ∈ N ∪ {0};) petals.

Roses with 4n + 2 petals can be obtained by using the Bernoulli lemniscate and its extensions. More precisely,

  • for n = 0, a two petals rose comes from the equation ρ = cos1/2(2 θ) (the Bernoulli lemniscate),

  • for n ≥ 1, a 4n + 2 petals rose comes from the equation ρ = cos1/2[(4n + 2) θ].

Further very general extensions of the Bernoulli lemniscate are presented in Section 4.1.


P. Butzer and F. Jongmans, in their biography of Chebyshev [14], assert that Pafnuty Lvovich Chebyshev (Figure 9) was the creator in St. Petersburg of the greatest Russian mathematical school before the revolution.

Figure 9

Pafnuty Lvovich Chebyshev (1821–1894).

Starting from the equations:

(eit)n=eint,(cost+isint)n=cos(nt)+isin(nt), (10)
and using the binomial expansion, we find:

By comparing these equations with (10), we find:

cos(nt)=h=0[n2](1)h(n2h)cosn2ht(1cos2t)h (11)
sin(nt)sint=h=0[n12](1)h(n2h+1)cosn2h1t(1cos2t)h. (12)

Putting x = cost, in Eqs. (11) and (12) we find two polynomials, in the x variable, of degrees respectively n and n − 1, which are the first and second kind Chebyshev polynomials [15,16]:

Tn(x):=cos(narccosx)=h=0[n2](1)h(n2h)xn2h(1x2)h, (13)

Un1(x):=sin(narccosx)sin(arccosx)=h=0[n12](1)h(n2h+1)xn2h1(1x2)h. (14)

3.1. Basic Properties of the Chebyshev Polynomials of the First Kind

The trigonometric equation

gives the recurrence relation:
Tn+1(x)=2xTn(x)Tn1(x). (15)

By using the initial values:

the subsequent polynomials are found:

Note that:

  • The leading coefficient of Tn(x) is equal to 2n−1.

  • The polynomials T2n(x) are even functions and the T2n+1(x) are odd functions.

  • nN,Tn(1)=1,Tn(1)=(1)n.

  • The n zeros of Tn(x) are real, distinct and internal to the interval [−1, 1].

More precisely, they are given by:

xn,k=cos((2k+1)nπ2),(k=0,1,,n1), (16)

In fact, we have:

  • The polynomials {Tn(x)} are orthogonal in the interval [−1, 1], with respect to the weight: (1 − x2)−1/2.

In fact, from the orthogonality of the cosine functions:

by the change of variable: t = arc cosx, we find:
11Tn(x)Tm(x)1x2dx=0,(ifmn). (17)

We furthermore have:

1111x2dx=π (18)
11Tn2(x)1x2dx=0πcos2(nt)dt=π2,(nN). (19)

3.2. Basic Properties of the Chebyshev Polynomials of the Second Kind

In a similar way, the same recurrence relation holds for the second kind polynomials:

Un+1(x)=2xUn(x)Un1(x). (20)

By using the initial values:

the subsequent polynomials are found:

The polynomials {Un(x)} are orthogonal in the interval [−1, 1], with respect to the weight: (1 − x2)1/2:

11Un(x)Um(x)1x2dx=0,(ifmn). (21)

We furthermore have:

11Un2(x)1x2dx=π2,(nN). (22)

Connections with the polynomials of the first kind

(1x2)Un1(x)=xTn(x)Tn+1(x),Tn(x)=Un(x)xUn1(x). (23)

The second kind Chebyshev polynomials play an important role in representing the powers of a 2 × 2 non singular matrix [17,18]. Extension of this polynomial family to the multivariate case has been considered for representing the powers of a r × r (r ≥ 3) non-singular matrix (see [18,19]).

Remark 3.1.

Chebyshev polynomials are a particular case of the Jacobi polynomials Pn(α,β)(x) , which are orthogonal in the interval [−1, 1] with respect to the weight (1 − x)α(1 + x)β. More precisely, the following equations hold:


Therefore, properties of the Chebyshev polynomials could be deduced in a more general framework of the hypergeometric functions. However, in this approach the connection with trigonomeric function disappears. In this article, dealing with elementary functions, we use only elementary methods, in order to make the topic accessible to a wider audience of readers. In such a way, we avoid to shoot flies with cannons.

Remark 3.2.

In connection with interpolation and quadrature problems, another couple of Chebyshev polynomials have been considered. They correspond to different choices of weights:


These were called by Gautschi [20] the third and fourth kind Chebyshev polynomials, and will be considered in what follows.

3.3. Basic Properties of the Chebyshev Polynomials of Third and Fourth Kind

The third and fourth kind Chebyshev polynomials have been studied and applied by several scholars (see e.g. [8,21]), because they are useful in quadrature rules, when the singularities occur only at one of the end points (+1 or −1) (see [22]). Furthermore, recently they have been applied in Numerical Analysis for solving high odd-order boundary value problems with homogeneous or nonhomogeneous boundary conditions [21].

The third and fourth kind Chebyshev polynomials are defined in [−1, 1] as follows:

Vn(x)=cos[(n+12)arccosx]cos[(arccosx)2],Wn(x)=sin[(n+12)arccosx]sin[(arccosx)2]. (24)

Since Wn(x) = (−1)nVn(−x), as it can be see by their graphs (Figures 10 and 11), the third and fourth kind Chebyshev polynomials are essentially the same polynomial set, but interchanging the ends of the interval [−1, 1].

Figure 10

Vk(x), k = 1, 2, 3, 4. 1. violet, 2. brown, 3. red, 4. blue

Figure 11

Wk(x), k = 1, 2, 3, 4. 1. violet, 2. brown, 3. red, 4. blue.

The orthogonality properties hold [8]:


One of the explicit advantages of Chebyshev polynomials of third and fourth kind is to estimate some definite integrals as

with the precision degree 2n − 1, by using the n interpolatory points xk=cos(2k1)π2n+1 , (k = 1, 2, ..., n), in the interval [−1, 1] [7,21,22].


A few graphs of Rhodonea curves are shown in Figures 1314.

Figure 12

G. G. Grandi and G. W. Leibniz.

Figure 13

Rhodonea cos(pqθ) .

Figure 14

Rodoneas cos(18θ)andcos(38θ) .

4.1. Cartesian vs Polar Coordinates

The study of curves is of great importance for the modelling of natural objects and has attracted many scholars [9,23]. In particular, the so-called superformula by Gielis [2,3,5] has made it possible to construct the most diverse figures by varying a few parameters (see also [4,6,24]).

By exploiting the correspondence, mentioned in Section 2.2, between curves in the cartesian or polar form, we can find further graphs both of symmetric and asymmetric type.

In what follows we show some symmetric graphs, which does not coincide with the shapes of Grandi curves, and several others which are not symmetric. These graphs generalize in a wide way both the Grandi curves and the Bernoulli lemniscate recalled in Section 2.3.

To this aim, we consider first circular functions of the type:


The corresponding polar curves are:

and, by choosing m, n, h, k as particular positive integers, we find the few graphs depicted in Figures 1517.

Of course we have more examples by changing sine by cosine or increasing the number of sine/cosine factors in the above trigonometric products.

The curves obtained in this way can be interpreted as a generalization of the Grandi roses.

A first example of this kind of figures can be found in Thompson [1] (Figure 504, page 1047), where the polar equation ρ = sin(θ/2)sin() can be found.

A infinite number of non-symmetric possibilities (generalizing in particular the Bernoulli lemniscate), are obtained by starting from the above circular functions multiplied by powers of the x variable (of polynomials Pq(x) of degree q), and taking positive rational numbers for the parameters m, n:

therefore, finding the polar equations:

The graphic results obtained in these cases recall the inflatable balloons of children’s games. Note that in this case, the resulting curve may depend on the considered interval.

A few examples of this type are shown in Figures 1822.

Figure 20

ρ = cos(θ/2)sin2(3θ) in [0, 2π] & in [0, 4π].

Figure 21

ρ = θ sin(θ/2)cos2(2θ) & ρ = θ sin(θ/2)cos2(3θ) in [0, 2π].

Figure 22

ρ = (θ2 + 1)sin(θ/2)cos(2θ) in [0, 2π] & in [0, 4π].


5.1. Basic properties of the First and Second Kind

We put, by definition:

Tpq(x)=cos(pqarccos(x)), (25)
1x2Upq(x)=sin(pqarccos(x)), (26)
where p and q are integer numbers, (q ≠ 0).

Note that definitions (25) and (26) hold even for negative indexes, that is for p/q < 0, according to the parity properties of the trigonometric functions.

5.2. Pseudo-Chebyshev Functions of the First Kind

The following theorems hold:

Theorem 5.1.

The pseudo-Chebyshev functions Tp/q(x) satisfy the recurrence relation

Tpq+1(x)=2xTpq(x)Tpq1(x). (27)

Proof - Write Eq. (27) in the form:

then use definition (25) and the trigonometric identity:

Theorem 5.2.

The first kind pseudo-Chebyshev functions Tp/q(x) satisfy the differential equation:

(1x2)yxy+(pq)2y=0. (28)

Proof - Note that

so that Eq. (28) follows.

5.3. Pseudo-Chebyshev Functions of the Second Kind

The following theorems hold:

Theorem 5.3.

The pseudo-Chebyshev functions Upq(x) satisfy the recurrence relation

Upq+1(x)=2xUpq(x)Upq1(x). (29)

Proof - Write Eq. (29) in the form:

then use definition (26) and the trigonometric identity:

Theorem 5.4.

The pseudo-Chebyshev functions y(x)=Upq(x) satisfy the differential equation

(1x2)y3xy+[(pq)21]y=0. (30)

The proof is obtained in a similar way to that of the first kind functions.

5.4. Basic Relations of the Pseudo-Chebyshev Functions of Third and Fourth Kind

According to Eq. (24), we put by definition:

Vpq(x)=cos[(pq+12)arccosx]cos[(arccosx)2],Wpq(x)=sin[(pq+12)arccosx]sin[(arccosx)2]. (31)

Theorem 5.5.

The third and fourth kind pseudo-Chebyshev functions are related to the 1st and 2nd kind ones by the equations:

Vpq(x)=Tpq(x)(1x)Upq(x),Wpq(x)=Tpq(x)+(1+x)Upq(x). (32)

Proof - It is sufficient to use the addition formulas for the cosine and sine functions.

Therefore, we can derive the equations:

Wpq(x)Vpq(x)=2Upq(x),Wpq(x)+Vpq(x)=2Tpq(x)+2xUpq(x). (33)

5.5. Some General Formulas

By using cosine addition formulas, putting:

mn=pq+rs, (34)
we find:
Tmn(x)=Tpq(x)Trs(x)(1x2)U(pq)1(x)U(rs)1(x), (35)
and by using the sine addition formulas:
Umn(x)=U(pq)1(x)Trs(x)+U(rs)1(x)Tpq(x). (36)

5.5.1. Particular results

T1(x)=T13(x)T23(x)(1x2)U23(x)U13(x). (37)
T1(x)=cos[313arccos(x)]=4T133(x)3T13(x). (38)
T2(x)=cos[323arccos(x)]=4T233(x)3T23(x). (39)
T23(x)=cos[213arccos(x)]=12sin2[13arccos(x)]=12(1x2)U23(x). (40)
U13(x)=sin[213arccos(x)]1x2=21x2sin[13arccos(x)]cos[13arccos(x)]=2U23(x)T13(x). (41)
U23(x)=sin[13arccos(x)]1x2=1T132(x)1x2. (42)

Combining the above equations, we find:

T1(x)=T13(x)T23(x)2T13(x)(1T132(x))=T13(x)(2T132(x)+T23(x)2). (43)

5.6. Links with the Pseudo-Chebyshev Functions

Actually, the definitions of the third and fourth kind Chebyshev polynomials are as follows:

Vn(x)=cos[(n+12)arccosx]cos[(arccosx)2]=21+xTn+12(x)=T121(x)Tn+12(x), (44)
Wn(x)=sin[(n+12)arccosx]sin[(arccosx)2]=21+x2Un12(x)=2T12(x)Un12(x). (45)

Therefore, we find the equations:

Vpq(x)=21+xTpq+12(x)=T121(x)Tpq+12(x), (46)
Wpq(x)=21+x2Upq12(x)=2T12(x)Upq12(x). (47)


In what follows, we consider the case of the half-integer degree, which seems to be the most interesting one, since the resulting pseudo-Chebyshev functions satisfy the orthogonality properties in the interval [−1, 1] with respect to the same weights of the corresponding Chebyshev polynomials [25].

Definition: Let, for any integer k:

Tk+12(x)=cos((k+12)arccos(x)),1x2Uk12(x)=sin((k+12)arccos(x)),1x2Vk+12(x)=cos((k+12)arccos(x)),Wk+12(x)=sin((k+12)arccos(x)). (48)

Note that the above definition holds even for k + 1/2 < 0, taking into account the parity properties of the circular functions.

The pseudo-Chebyshev functions Tk+1/2 (x), Uk−1/2 (x), Vk+1/2 (x) and Wk+1/2 (x) can be represented, in terms of the third and fourth kind Chebyshev polynomials as follows:

Tk+12(x)=1+x2Vk(x),1x2Uk12(x)=12(1+x)Wk(x),1x2Vk+12(x)=12(1x)Vk(x),Wk+12(x)=1x2Wk(x). (49)

We will show that, in the case of half-integer degree, the pseudo-Chebyshev functions satisfy not only the recurrence relations and differential equations analogues to the classical ones, but even the orthogonality properties.

6.1. Orthogonality of the Tk+1/2 (x) and Uk+1/2 (x) Functions

A few graphs of the Tk+12 functions are shown in Figure 23.

Figure 23

Tk+1/2 (x), k = 1, 2, 3, 4. 1. green, 2. red, 3. blue, 4. orange.

Theorem 6.1.

The pseudo-Chebyshev functions Tk+1/2 (x) satisfy the orthogonality property:

11Th+12(x)Tk+12(x)11x2dx=0,(hk), (50)
where h, k are integer numbers,
11Tk+122(x)11x2dx=π2. (51)

A few graphs of the Uk+1/2 functions are shown in Figure 24.

Figure 24

Uk+1/2 (x), k = 1, 2, 3, 4. 1. green; 2. red; 3. blue; 4. orange.

Theorem 6.2.

The pseudo-Chebyshev functions Uk+12(x) satisfy the orthogonality property:

11Uh+12(x)Uk+12(x)1x2dx=0,(mn), (52)
where h, k are integer numbers,
11Uk+122(x)1x2dx=π2. (53)

Proof - We prove only Theorem 6.1, since the proof of Theorem 6.2 is similar.

From the Werner formulas, we have:


6.2. The Third and Fourth Kind Pseudo-Chebyshev Functions

The results of this section are based on the excellent survey by Aghigh et al. [8]. By using that article, it is possible to derive, in an almost trivial way, the links among the pseudo-Chebyshev functions and the third and fourth kind Chebyshev polynomials.

We recall here only the principal properties, without proofs. Proofs and other properties are reported in Cesarano et al. [27]. In Figures 25 and 26, we show the graphs of the first few third and fourth kind pseudo-Chebyshev functions.

Figure 25

Vk+1/2 (x), k = 1, 2, 3, 4, 5. 1. grey, 2. red, 3. blue, 4. orange, 5. violet.

Figure 26

Wk+1/2 (x), k = 1, 2, 3, 4, 5. 1. red, 2. blue, 3. orange, 4. violet, 5. grey.

6.3. The Third Kind Pseudo-Chebyshev Vk+1/2

The third kind pseudo-Chebyshev functions satisfy the recurrence relation:

{Vk+12(x)=2xVk12(x)Vk32(x),V±12(x)=12(1x). (54)

Theorem 6.3.

The pseudo-Chebyshev functions Vk+1/2(x) verify the differential equation:

(1x2)y k3xy k+[(k+12)21]yk=0. (55)

The orthogonality property holds:

11Vh+12(x)Vk+12(x)1x2dx=0,(hk), (56)
11Vk+122(x)1x2dx=π2. (57)

6.4. The Fourth Kind Pseudo- Chebyshev Wk+1/2

The fourth kind pseudo-Chebyshev functions satisfy the recurrence relation:

{Wk+12(x)=2xWk12(x)Wk32(x),W±12(x)=±1x2. (58)

Theorem 6.4.

The pseudo-Chebyshev functions Wk+1/2 (x) verify the differential equation:

(1x2)ykxyk+(k+12)2yk=0, (59)

The orthogonality property holds:

11Wh+12(x)Wk+12(x)11x2dx=0,(hk), (60)
11Wk+122(x)11x2dx=π2. (61)

6.5. Explicit Forms

Theorem 6.5.

It is possible to represent explicitly the pseudo-Chebyshev functions as follows:

Tk+12(x)=1+x2h=0k(1)h(2k+12h)(1x2)h(1+x2)kh,Uk12(x)=12(1+x)h=0k(1)h(2k+12h+1)(1x2)h(1+x2)kh,Vk+12(x)=12(1x)h=0k(1)h(2k+12h)(1x2)h(1+x2)kh,Wk+12(x)=1x2h=0k(1)h(2k+12h+1)(1x2)h(1+x2)kh. (62)

6.6. Location of Zeros

By Eq. (48), the zeros of the pseudo-Chebyshev functions Tk+1/2(x) and Vk+1/2 (x) are given by

xk,h=cos((2h1)π2k+1),(h=1,2,,k), (63)
and the zeros of the pseudo-Chebyshev functions Uk+1/2(x) and Wk+1/2 (x) are given by
xk,h=cos(2hπ2k+1),(h=1,2,...,k), (64)
furthermore, the Wk+1/2 (x) functions always vanish at the end of the interval [−1, 1].

Remark 6.6.

More technical properties as the Hypergeometric representations and the Rodrigues-type formulas are reported in Cesarano et al. [26].

6.7. Links with First and Second Kind Chebyshev Polynomials

Theorem 6.7.

The pseudo-Chebyshev functions are connected with the first and second kind Chebyshev polynomials by means of the equations:

Tk+12(x)=T2k+1(1+x2)=T2k+1(T1/2(x)),Uk12(x)=11+x12(1x)U2k(1+x2),Vk+12(x)=11x2T2k+1(1+x2)=11x2Tk+12(x),Wk+12(x)=1x2U2k(1+x2)=(1x2)Uk12(x). (65)

Proof - The results follow from the equations:

Vk(x)=21+xT2k+1(1+x2),Wk(x)=U2k(1+x2), (66)
(see [8]), by using definitions (49).

Remark 6.8.

Note that the first equation in (65), extends the known nesting property verified by the first kind Chebyshev polynomials:

Tm(Tn(x))=Tmn(x). (67)

This property, already considered in Brandi and Ricci [27] for the first kind pseudo-Chebyshev functions, actually holds in general, as a consequence of the definition Tk (x) = cos(k arccos(x)). Note that this composition identity even holds for the first kind Chebyshev polynomials in several variables, as it was proven in Ricci [28].


The growth of living organisms is often described by the Bernoulli’s logarithmic spiral, which is one of the first examples of fractals. The study of natural forms (see e.g. Bini et al. [29]) is commonly associated with mathematical entities like extensions of Lamé’s curves, Grandi’s roses, Bernoulli’s lemniscate, etc.

In this article, it has been shown that innumerable plane forms can be described by means of polar equations that extend some of the above mentioned geometrical entities. Moreover, the consideration of Grandi’s roses in the case of fractional indexes gives rise, in a natural way, to mathematical functions that generalize to the case of fractional degree the classical Chebyshev polynomials. The resulting functions, called of pseudo-Chebyshev type, verify many properties of the corresponding polynomials and, in the case of half-integer degree, also the orthogonality properties.


[1]DW Thompson, On growth and form, Cambridge, Cambridge University Press, 1942, pp. 116.
[3]J Gielis, The geometrical beauty of plants, Atlantis Press, Springer Nature, Basingstoke, UK, 2017.
[5]J Gielis, S Haesen, and L Verstraelen, Universal natural shapes: from the supereggs of piet hein to the cosmic egg of georges lemaître, Kragujevac J Math, Vol. 28, 2005, pp. 55-67.
[9]EH Lockwood, A book of curves, Cambridge University Press, Cambridge, 1961.
[11]J Gielis, D Caratelli, P Shi, and PE Ricci, A note on spirals and curvature, 2019. (This issue).
[12]TL Heath, The works of Archimedes, Cambridge University Press, Cambridge, 1897, pp. 524.
[13]RC Archibald, Notes on the logarithmic spiral, golden section and the Fibonacci series, J Hambidge (editor), Dynamic symmetry, Yale University Press, New Haven, 1920, pp. 146-57.
[15]TJ Rivlin, The Chebyshev polynomials, John Wiley & Sons, New York, 1974.
[16]HM Srivastava and HL Manocha, A treatise on generating functions, Halsted Press (Ellis Horwood Limited, Chichester), John Wiley and Sons, New York, Chichester, Brisbane and Toronto, 1984.
[17]PE Ricci, Alcune osservazioni sulle potenze delle matrici del secondo ordine e sui polinomi di Tchebycheff di seconda specie, Atti Accad Sci Torino, Vol. 109, 1975, pp. 405-10.
[18]PE Ricci, Sulle potenze di una matrice, Rend Mat, Vol. 9, 1976, pp. 179-94.
[19]PE Ricci, I polinomi di Tchebycheff in più variabili, Rend Mat, Vol. 11, 1978, pp. 295-327.
[20]W Gautschi, On mean convergence of extended Lagrange interpolation, J Comput Appl Math, Vol. 43, 1992, pp. 19-35.
[22]JC Mason and DC Handscomb, Chebyshev polynomials, Chapman & Hall/CRC, Boca Raton, London, New York, Washington D.C., 2003.
[27]P Brandi and PE Ricci, Some properties of the pseudo-Chebyshev polynomials of half-integer degree, Tbilisi Math J, Vol. 12, 2019, pp. 111-21.
[28]PE Ricci, Una proprietà iterativa dei polinomi di Chebyshev di prima specie in più variabili, Rend Mat Appl, Vol. 6, 1986, pp. 555-63.
[29]D Bini, C Cherubini, S Filippi, A Gizzi, and PE Ricci, On the universality of spiral waves, Commun Comput Phys, Vol. 8, 2010, pp. 610-22.
Growth and Form
1 - 1
20 - 32
Publication Date
ISSN (Online)
10.2991/gaf.k.200124.005How to use a DOI?
© 2020 The Authors. Published by Atlantis Press SARL
Open Access
This is an open access article distributed under the CC BY-NC 4.0 license (

Cite this article

AU  - Paolo Emilio Ricci
PY  - 2020
DA  - 2020/02/23
TI  - Chebyshev Polynomials, Rhodonea Curves and Pseudo-Chebyshev Functions. A Survey
JO  - Growth and Form
SP  - 20
EP  - 32
VL  - 1
IS  - 1
SN  - 2589-8426
UR  -
DO  - 10.2991/gaf.k.200124.005
ID  - Ricci2020
ER  -