Open AccessArticle

Repeated Derivatives of Hyperbolic Trigonometric Functions and Associated Polynomials

Giuseppe Dattoli

^1,†,

Silvia Licciardi

^1,*,†

Rosa Maria Pidatella

^2,† and

Elio Sabia

^1,†

ENEA—Frascati Research Center, Via Enrico Fermi 45, 00044 Frascati (Rome), Italy

Department of Mathematics and Computer Sciences, University of Catania, Viale A. Doria 6, 95125 Catania, Italy

Author to whom correspondence should be addressed.

Axioms 2019, 8(4), 138; https://doi.org/10.3390/axioms8040138

Submission received: 8 November 2019 / Revised: 4 December 2019 / Accepted: 5 December 2019 / Published: 6 December 2019

(This article belongs to the Special Issue Mathematical Analysis and Applications II)

Download Review Reports Versions Notes

Abstract

Elementary problems as the evaluation of repeated derivatives of ordinary transcendent functions can usefully be treated with the use of special polynomials and of a formalism borrowed from combinatorial analysis. Motivated by previous researches in this field, we review the results obtained by other authors and develop a complementary point of view for the repeated derivatives of

sec (.)

tan (.)

and for their hyperbolic counterparts.

Keywords:

operators theory 44A99, 47B99, 47A62; special functions 33C52, 33C65, 33C99, 33B10, 33B15; Stirling numbers and Touchard polynomials 11B73

1. Introduction

The problem of finding closed forms for the repeated derivatives of trigonometric functions as tangent and secant, even though being an apparently elementary issue, has been solved in relatively recent times in ref. [1]. Inspired by this work, a significant amount of research has been subsequently developed. In ref. [2] the proof of the results of [1] was reformulated in terms of a procedure exploiting the “Zeons” Algebra [3]. In [4,5] the authors addressed this study by employing a class of polynomials (the derivative polynomials (

D P

) introduced in refs. [6,7]) to reformulate the derivation and eventually get a set of fairly simple formulae, providing the successive derivatives of

tan, cot, sec, csc, \dots

along with those of their hyperbolic counterparts.

In this note we develop a point of view not dissimilar from that of ref. [4,5]. We provide straightforward results in terms of a single family of Legendre like polynomials and comment on the two forms of

D P

introduced in [3,4,5].

The repeated derivatives of composite functions

F (x) = f (g (x))

is a well-established topic in calculus. The formulation of a procedure allowing the derivation of a formula comprising all the possible cases was established in the

X I X

century [8] and opened important avenue of research in combinatorics [9] and umbral calculus [10,11] as well.

The problem is particularly interesting, encompasses different topics in analysis, including special polynomials like those belonging to the Touchard family [12] and special numbers like the generalized Stirling forms [13,14] of crucial importance in combinatorial analysis. The repeated derivatives of the Gaussian function are those of a composite function in which

f (.)

is an exponential and

g (.)

a quadratic function. The relevant expression leads to the Hermite polynomials as auxiliary tool, to get a synthetic expression for any order of the derivative [15]. Within the same context, Bell polynomials emerge whenever one is interested in the derivatives of

F (x) = e^{g (x)}

[9,16,17]. The generalization to the case of

f (.)

, provided by a generic infinitely differentiable function, and

g (.)

, a quadratic form, has been discussed in refs. [18,19] where the problem has been solved by the use of generalized nested forms of Hermite polynomials.

The hyperbolic secant is a composite function too in which

f (.) = {(.)}^{- 1}

g (.) = cosh (.)

. In refs. [2,4,5,20] interesting speculations have been presented on the study of the relevant repeated derivatives and of the associated auxiliary polynomials. In this paper we address the same problem, within a different context. Before entering the specific elements of the discussion we review the formalism we are going to exploit in this paper.

A pivotal role within the context of repeated derivatives is played by the Stirling number of second kind [13,14]. They will be introduced using the Touchard polynomials which are defined through the Rodriguez type formula [15,21] where, for

n \geq 0

, we set

T_{n} (x) = e^{- x} {(x {\hat{D}}_{x})}^{n} e^{x}

(1)

where

{\hat{D}}_{x}

is the derivative operator

\frac{d}{d x}

. They can be written in explicit form by the use of the following expression [13]

{(x {\hat{D}}_{x})}^{n} = \sum_{r = 0}^{n} S_{2} (n, r) x^{r} {\hat{C}}_{x}^{r},

(2)

with

S_{2} (n, r)

being the Stirling number of second kind (we introduce the notation

S_{2} (l, m)

instead of the usually symbol

\{\begin{matrix} l \\ m \end{matrix}\}

; to indicate functions domains, we will use conventional symbols as

N, Z, R

for Natural, Integer or Real numbers respectively including the zero). Ref. [9,12] defined as

S_{2} (l, m) = \frac{1}{m!} \sum_{j = 0}^{m} {(- 1)}^{m - j} (\binom{m}{j}) j^{l} .

(3)

According to Equations (1) and (2), the Touchard polynomials are explicitly given by

T_{n} (x) = \sum_{r = 0}^{n} S_{2} (n, r) x^{r}

(4)

and the numbers

S_{2} (n, k)

are therefore the coefficients of the polynomials.

A fairly straightforward application of the previous computational tool is provided by the evaluation of a closed expression for higher order derivatives.

Theorem 1.

We have that,

\forall m \in N

{\hat{D}}_{x}^{m} f (e^{x}) = \sum_{r = 0}^{m} S_{2} (m, r) e^{x r} f^{(r)} (e^{x})

(5)

where

f^{(r)} (ξ)

denotes the r-order derivative.

Proof.

Let

I_{m} (x) = {\hat{D}}_{x}^{m} f (e^{x}),

(6)

by setting

e^{x} = ξ

, we find

{\hat{D}}_{x}^{m} f (e^{x}) = {(ξ {\hat{D}}_{ξ})}^{m} f (ξ) = \sum_{r = 0}^{m} S_{2} (m, r) ξ^{r} f^{(r)} (ξ)

(7)

and in conclusion we obtain

{\hat{D}}_{x}^{m} f (e^{x}) = \sum_{r = 0}^{m} S_{2} (m, r) e^{x r} f^{(r)} (e^{x})

in which it is understood that

f^{(r)} (e^{x}) = {\hat{D}}_{ξ}^{r} f (ξ) ∣_{ξ = e^{x}}

. □

We next define a family of auxiliary polynomials

P_{n} (x, y)

and show how they can be very useful for the computation of the repeated derivatives of

{tan}^{- 1} (x)

(we use the notation

{tan}^{- 1} (x)

for

arctan (x)

Theorem 2.

For any

x, y \in R, n \in N

P_{n} (x, y) = n! \sum_{r = 0}^{⌊ \frac{n}{2} ⌋} \frac{x^{n - 2 r} y^{r} (n - r)!}{(n - 2 r)! r!},

(8)

a family of two variable polynomials loosely ascribed to the Legendre [22,23] family, then

{\hat{D}}_{x}^{n} ({tan}^{- 1} (x)) = \frac{1}{1 + x^{2}} P_{n - 1} (- \frac{2 x}{1 + x^{2}}, - \frac{1}{1 + x^{2}}), \forall n \in N_{0} .

(9)

Proof.

We start from the computation of the quantity

K_{n} (x) = {\hat{D}}_{x}^{n} (\frac{1}{1 + x^{2}}), \forall x \in R, \forall n \in N

(10)

which upon the use of the Laplace transform method can be written as

K_{n} (x) = {\hat{D}}_{x}^{n} \int_{0}^{\infty} e^{- s (1 + x^{2})} d s = \int_{0}^{\infty} e^{- s} {\hat{D}}_{x}^{n} e^{- s x^{2}} d s .

(11)

The n-th order derivative inside the integral sign can be explicetely worked in terms of two variable Hermite polynomials [15]

H_{n} (x, y) = n! \sum_{r = 0}^{⌊ \frac{n}{2} ⌋} \frac{x^{n - 2 r} y^{r}}{(n - 2 r)! r!},

(12)

namely the auxiliary polynomials for the repeated derivatives of Gaussian functions according to the identity

{\hat{D}}_{x}^{n} e^{- a x^{2}} = H_{n} (- 2 a x, - a) e^{- a x^{2}} .

(13)

Accordingly we find

\begin{matrix} K_{n} (x) & = \int_{0}^{\infty} e^{- s} H_{n} (- 2 x s, - s) e^{- s x^{2}} d s = \int_{0}^{\infty} e^{- σ} H_{n} (- 2 \frac{x σ}{1 + x^{2}}, - \frac{σ}{1 + x^{2}}) \frac{1}{1 + x^{2}} d σ = \\ = \frac{n!}{1 + x^{2}} \sum_{r = 0}^{⌊ \frac{n}{2} ⌋} \frac{{(\frac{- 2 x}{1 + x^{2}})}^{n - 2 r} {(\frac{- 1}{1 + x^{2}})}^{r}}{(n - 2 r)! r!} \int_{0}^{\infty} e^{- σ} σ^{n - r} d σ = \frac{1}{1 + x^{2}} P_{n} (- \frac{2 x}{1 + x^{2}}, - \frac{1}{1 + x^{2}}) . \end{matrix}

(14)

The proof of the identity (9) is readily achieved by noting that

{\hat{D}}_{x} ({tan}^{- 1} (x)) = (\frac{1}{1 + x^{2}}),

(15)

which eventually provides

{\hat{D}}_{x}^{n} ({tan}^{- 1} (x)) = K_{n - 1} (x) = \frac{1}{1 + x^{2}} P_{n - 1} (- \frac{2 x}{1 + x^{2}}, - \frac{1}{1 + x^{2}}), \forall n \in N_{0} .

□

Corollary 1.

If we introduce the further notation

K_{n}^{ν} (x) : = {\hat{D}}_{x}^{n} (\frac{1}{{(1 + x^{2})}^{ν}}),

(16)

an extension of the procedure we have just envisaged allows the further results

\begin{matrix} K_{n}^{ν} (x) = \frac{1}{{(1 + x^{2})}^{ν}} P_{n}^{ν} (- \frac{2 x}{1 + x^{2}}, - \frac{1}{1 + x^{2}}), \\ P_{n}^{ν} (w, z) = \frac{n!}{Γ (ν)} \sum_{r = 0}^{⌊ \frac{n}{2} ⌋} \frac{w^{n - 2 r} z^{r} Γ (ν + n - r)}{(n - 2 r)! r!}, \forall ν \in R . \end{matrix}

(17)

Proof.

\forall ν \in R

\begin{matrix} K_{n}^{ν} (x) & = {\hat{D}}_{x}^{n} (\frac{1}{{(1 + x^{2})}^{ν}}) = {\hat{D}}_{x}^{n} \frac{1}{Γ (ν)} \int_{0}^{\infty} s^{ν - 1} e^{- s (1 + x^{2})} d s = \frac{1}{Γ (ν)} \int_{0}^{\infty} s^{ν - 1} e^{- s} {\hat{D}}_{x}^{n} e^{- s x^{2}} d s = \\ = \frac{1}{Γ (ν)} \int_{0}^{\infty} s^{ν - 1} e^{- s} H_{n} (- 2 x s, - s) e^{- s x^{2}} d s = \\ = \frac{1}{Γ (ν)} \int_{0}^{\infty} {(\frac{σ}{1 + x^{2}})}^{ν - 1} e^{- σ} H_{n} (- 2 \frac{x σ}{1 + x^{2}}, - \frac{σ}{1 + x^{2}}) \frac{1}{1 + x^{2}} d σ = \\ = \frac{n!}{Γ (ν) {(1 + x^{2})}^{ν}} \sum_{r = 0}^{⌊ \frac{n}{2} ⌋} \frac{{(\frac{- 2 x}{1 + x^{2}})}^{n - 2 r} (\frac{- 1}{1 + x^{2}})}{(n - 2 r)! r!} \int_{0}^{\infty} e^{- σ} σ^{ν + n - r - 1} d σ = \\ = \frac{1}{{(1 + x^{2})}^{ν}} P_{n}^{ν} (- \frac{2 x}{1 + x^{2}}, - \frac{1}{1 + x^{2}}) . \end{matrix}

□

Remark 1.

Corollary 1 can be exploited to write the repeated derivatives of

{cos}^{- 1} (x)

in the form

{\hat{D}}_{x}^{n} {cos}^{- 1} (x) = - {\hat{D}}_{x}^{n - 1} (\frac{1}{\sqrt{1 - x^{2}}}) = - \frac{1}{\sqrt{1 - x^{2}}} P_{n - 1}^{\frac{1}{2}} (\frac{2 x}{1 - x^{2}}, \frac{1}{1 - x^{2}}), \forall n \in N_{0} .

(18)

The previous examples have shown the interplay between special polynomials of the

P_{n} (x, y)

and the derivatives of the inverse of trigonometric functions.

In the following we will see that the same polynomials are of central importance for the evaluation of the repetead derivatives of trigonometric functions. They play the role of auxiliary polynomials for the lorentzian type functions and inverse trigonometric functions as well. The relevant elements of contact with the

D P

will be discussed in the final section of the paper.

2. Higher Order Derivatives of Trigonometric Functions

The starting point of the discussion of this section is the derivation of a closed form for the quantity

{\hat{D}}_{x}^{m} (sec h (x))

\forall m \in N_{0}

Proposition 1.

\forall m \in N_{0}, \forall x \in [0, 2 π]

\begin{matrix} {\hat{D}}_{x}^{m} (sec h (x)) = & sec h (x) \sum_{k = 0}^{m} S_{2} (m, k) e^{(k - 1) x} (e^{x} P_{k} (- sec h (x), - \frac{1}{2 e^{x}} sec h (x)) + \\ + k P_{k - 1} (- sec h (x), - \frac{1}{2 e^{x}} sec h (x))) . \end{matrix}

(19)

Proof.

Let

sec h (x) = \frac{2}{e^{x} + e^{- x}}

(20)

then, after setting

e^{x} = ξ

, we obtain

\forall m \in N_{0}

{\hat{D}}_{x}^{m} (sec h (x)) = {\hat{D}}_{log ξ}^{m} (\frac{2 ξ}{1 + ξ^{2}}) = 2 {(ξ {\hat{D}}_{ξ})}^{m} (\frac{ξ}{1 + ξ^{2}}) .

(21)

The straightforward applications of the proof of the Theorem 1 yields (we omit the argument of the

P_{n}

polynomials for simplicity)

2 {(ξ {\hat{D}}_{ξ})}^{m} (\frac{ξ}{1 + ξ^{2}}) = \frac{2}{1 + ξ^{2}} \sum_{k = 0}^{m} S_{2} (m, k) ξ^{k} (ξ P_{k} (., .) + k P_{k - 1} (., .))

(22)

which, eventually, leads to

\begin{matrix} {\hat{D}}_{x}^{m} (sec h (x)) = & sec h (x) \sum_{k = 0}^{m} S_{2} (m, k) e^{(k - 1) x} (e^{x} P_{k} (- sec h (x), - \frac{1}{2 e^{x}} sec h (x)) + \\ + k P_{k - 1} (- sec h (x), - \frac{1}{2 e^{x}} sec h (x))) . \end{matrix}

□

We can go even further and obtain a closed form for

{\hat{D}}_{x}^{m} {(sec h (x))}^{ν}

\forall m \in N, \forall ν \in R

Corollary 2.

\forall m \in N, \forall ν \in R

{\hat{D}}_{x}^{m} (sec h^{ν} (x)) = sec h^{ν} (x) \sum_{r = 0}^{m} S_{2} (m, r) e^{x r} \sum_{s = 0}^{r} (\binom{r}{s}) \frac{n!}{(n - r + s)!} e^{- x (r - s)} P_{s}^{ν} (- sec h (x), - \frac{1}{2 e^{x}} sec h (x))

(23)

where

P_{m}^{k} (., .)

are the polynomials in Equation (17).

By using the same techniques of the Proposition 1, Corollary 2 and the opportune identity

sec (x) = \frac{2}{e^{i x} + e^{- i x}}

, we can move on and provide the higher order derivatives for other trigonometric circular functions.

Example 1.

\forall m \in N_{0}

\forall x \in [0, 2 π]

\begin{matrix} {\hat{D}}_{x}^{m} (sec (x)) = & i^{m} sec (x) \sum_{r = 0}^{m} S_{2} (m, r) e^{i x (r - 1)} (e^{i x} P_{r} (- sec (x), - \frac{1}{2 e^{i x}} sec (x)) + \\ + r P_{r - 1} (- sec (x), - \frac{1}{2 e^{i x}} sec (x))) . \end{matrix}

(24)

Example 2.

\forall m \in N

\forall x \in [0, 2 π] : x \neq \frac{π}{2} + k π, k \in Z

{\hat{D}}_{x}^{m} (tan (x)) = \sum_{s = 0}^{m} (\binom{m}{s}) sin (x + (m - s) \frac{π}{2}) {\hat{D}}_{x}^{s} (sec (x)) .

(25)

Example 3.

\forall m \in N

\forall x \in [0, 2 π] : x \neq k π, k \in Z

{\hat{D}}_{x}^{m} (cot (x)) = \sum_{s = 0}^{m} (\binom{m}{s}) cos (x + (m - s) \frac{π}{2}) {\hat{D}}_{x}^{s} (sec (x - \frac{π}{2})) .

(26)

Remark 2.

The use of the Leibniz rule in Equations (25) or (26) can be avoided by setting

tan (x) = i \frac{1 - ξ^{2}}{1 + ξ^{2}} = i (\frac{2}{1 + ξ^{2}} - 1), ξ = e^{i x}

(27)

and eventually ending up with

{\hat{D}}_{x}^{m} (tan (x)) = i^{m + 1} (sec (x) \sum_{r = 0}^{m} S_{2} (m, r) e^{i x (r - 1)} P_{r} (- sec (x), - \frac{1}{2 e^{i x}} sec (x)) - δ_{m, 0}),

(28)

where

δ_{m, 0}

is Kronecker’s delta, indeed

\begin{matrix} {\hat{D}}_{x}^{m} (tan (x)) & = i^{m} {(ξ {\hat{D}}_{ξ})}^{m} (i (\frac{2}{1 + ξ^{2}} - 1)) = 2 i^{m + 1} {(ξ {\hat{D}}_{ξ})}^{m} (\frac{1}{1 + ξ^{2}} - \frac{1}{2}) = \\ = 2 i^{m + 1} \sum_{r = 0}^{m} S_{2} (m, r) ξ^{r} {(\frac{1}{1 + ξ^{2}} - \frac{1}{2})}^{(r)} = \\ = 2 i^{m + 1} \sum_{r = 0}^{m} S_{2} (m, r) ξ^{r} ({(\frac{1}{1 + ξ^{2}})}^{(r)} + {(- \frac{1}{2})}^{(r)}) = \\ = 2 i^{m + 1} ((\sum_{r = 0}^{m} S_{2} (m, r) ξ^{r} \frac{1}{1 + ξ^{2}} P_{r} (., .)) + (\sum_{r = 0}^{m} S_{2} (m, r) ξ^{r} {(- \frac{1}{2})}^{(r)})) = \\ = i^{m + 1} (sec (x) \sum_{r = 0}^{m} S_{2} (m, r) e^{i x (r - 1)} P_{r} (- sec (x), - \frac{1}{2 e^{i x}} sec (x)) - δ_{m, 0}) . \end{matrix}

In the forthcoming section we will discuss the comparison with previous papers and present possible developments along the lines we have indicated.

3. Final Comments

In the previous two sections we have dealt with a general procedure useful to derive closed formulae for the repeated derivatives of circular and hyperbolic functions and of their inverse. We have underscored the importance of the polynomials

P_{n} (x, y)

which play a role analogous to the

D P

introduced in refs. [6,7] and used in [4,5] for analogous occurrences.

The

P_{n} (x, y)

are two variable polynomials defined in terms of the Laplace transform of the Hermite Kampé dé Fériét family [24]. They have been loosely defined Legendre-like and can be reduced to more familiar forms, by noting that

\begin{matrix} P_{n} (x, y) = y^{\frac{n}{2}} P_{n} (\frac{x}{\sqrt{- y}}), \\ P_{n} (z) = n! \sum_{r = 0}^{⌊ \frac{n}{2} ⌋} \frac{{(- 1)}^{r} z^{n - 2 r} (n - r)!}{(n - 2 r)! r!} \end{matrix}

(29)

with

P_{n} (z)

satisfying the generating function

\sum_{n = 0}^{\infty} \frac{t^{n}}{n!} P_{n} (z) = \frac{1}{1 - t z + t^{2}}, \forall t, z \in R : ∣ t^{2} - t z ∣ < 1

(30)

and

\frac{1}{n!} P_{n} (2 x, - 1) = \sum_{k = 0}^{⌊ \frac{n}{2} ⌋} \frac{{(- 1)}^{k} (n - k)! {(2 x)}^{n - 2 k}}{(n - 2 k)! k!} = U_{n} (x) .

(31)

The link with previous papers addressing the same problem treated here, in particular with that developed in refs. References [4,5] can be obtained by using the same steps suggested by Cvijović or Boyadzhiev. We assume that an identity of the type

{\hat{D}}_{x}^{n} (tan (x)) = Π_{n} (tan (x))

(32)

with

\prod_{n} (.)

being not yet specified polynomials, (in refs. [4,5] they are denoted by

P_{n} (x)

, we have used the capital Greek letter to avoid confusion with the

P_{n} (x, y)

polynomials defined in this paper) does holds. The polynomials

Π_{n} (.)

can be determined by the change of variable

tan (x) = ξ

which allows to transform Equation (32) into

{[(1 + ξ^{2}) {\hat{D}}_{ξ}]}^{n} (ξ) = Π_{n} (ξ)

(33)

which is a kind of Rodrigues type relation [25] defining the polynomials

Π_{n} (ξ)

The relevant generating function can be obtained by multiplying both sides of Equation (33) by

\frac{t^{n}}{n!}

, by summing up over the index n and ending up with

e^{t [(1 + ξ^{2}) {\hat{D}}_{ξ}]} ξ = \sum_{n = 0}^{\infty} \frac{t^{n}}{n!} Π_{n} (ξ) .

(34)

The use of the Lie derivative identity [26]

e^{t [(1 + ξ^{2}) {\hat{D}}_{ξ}]} f (ξ) = f (\frac{ξ cos (t) + sin (t)}{cos (t) - ξ sin (t)})

(35)

finally yields

\sum_{n = 0}^{\infty} \frac{t^{n}}{n!} Π_{n} (ξ) = \frac{ξ + tan (t)}{1 - ξ tan (t)}

(36)

which is the generating function given in [4,5].

The second family of

D P

can be defined by the use of the same procedure. According to refs. [4,5] they are implicitly defined by the condition

{\hat{D}}_{x}^{n} (sec (x)) = sec (x) Q_{n} (tan (x))

(37)

which, by the use of the same change of variable leading to Equation (33), yields

Q_{n} (ξ) = \frac{1}{\sqrt{1 + ξ^{2}}} {((1 + ξ^{2}) {\hat{D}}_{ξ})}^{n} \sqrt{1 + ξ^{2}}

(38)

which can be straightforwardly exploited to derive the relevant properties. The use of the identity (35) eventually yields the associated generating functions

\sum_{n = 0}^{\infty} \frac{t^{n}}{n!} Q_{n} (ξ) = \frac{sec (t)}{1 - ξ tan (t)}

(39)

which is the same reported in refs. [4,5].

By recalling that the Hoppe formula writes

\begin{matrix} {\hat{D}}_{t}^{m} f (g (t)) = \sum_{k = 0}^{m} \frac{1}{k!} f^{(k)} (σ) ∣_{σ = g (t)} A_{m, k}, \\ A_{m, k} = \sum_{j = 0}^{k} (\binom{k}{j}) {(- 1)}^{k - j} g {(t)}^{k - j} {\hat{D}}_{t}^{m} (g {(t)}^{j}), \forall m, k \in N, \end{matrix}

(40)

for the case of

sec (x)

, we find

{\hat{D}}_{x}^{m} sec (x) = \sum_{k = 0}^{m} {sec}^{k + 1} (x) \sum_{j = 0}^{k} (\binom{k}{j}) {(- 1)}^{j} {cos}^{k - j} (x) {\hat{D}}_{x}^{m} ({cos}^{j} (x))

(41)

which appears impractical since it needs the repeated derivatives of an integer power of the cosine function.

Before closing the paper we consider worth to underscore the possible use of

D P

to evaluate the derivatives of the type

λ_{m, j} (x) = {\hat{D}}_{x}^{m} ({cos}^{j} (x)) .

(42)

Making the ansatz that

{\hat{D}}_{x}^{m} ({cos}^{j} (x)) = {(- 1)}^{j} Λ_{m, j} (tan (x))

(43)

and, by using the same procedure leading to Equation (33), we can identify the function

Λ_{m, j} (ξ) = {[(1 + ξ^{2}) {\hat{D}}_{ξ}]}^{m} (\frac{1}{{(\sqrt{1 + ξ^{2}})}^{j}})

(44)

which is specified by the generating function

\sum_{m = 0}^{\infty} \frac{t^{m}}{m!} Λ_{m, j} (ξ) = \frac{{[1 - ξ tan (t)]}^{j}}{{[(1 + ξ^{2}) (1 + {tan}^{2} (t))]}^{\frac{j}{2}}} .

(45)

It is worth to consider the alternative assumption

{\hat{D}}_{x}^{m} ({cos}^{j} (x)) = {(- 1)}^{m} Δ_{m, j} (cos (x))

(46)

which, after the change of variable

ξ = cos (x)

, yields for the function

Δ_{m, j} (.)

the operational definition

Δ_{m, j} (ξ) = {(- \sqrt{1 - ξ^{2}} {\hat{D}}_{ξ})}^{m} ξ^{j}

(47)

The use of the Lie derivative identity [26]

e^{- t (\sqrt{1 - ξ^{2}} {\hat{D}}_{ξ})} f (ξ) = f (ξ cos (t) - \sqrt{1 - ξ^{2}} sin (t))

(48)

yields, for

Δ_{m, j} (ξ)

, the generating function

\sum_{m = 0}^{\infty} \frac{t^{m}}{m!} Δ_{m, j} (ξ) = {(ξ - \sqrt{1 - ξ^{2}} tan (t))}^{j} {cos}^{j} (t) .

(49)

The two methods we have discussed in this paper, namely the procedure based on the

D P

of refs. [4,5] or on the support polynomials

P_{n} (x, y)

, are complementary. There are no prevailing reasons to prefer one or the other method. The formulae associated with the first are more synthetic but

Π_{n} (ξ), Q_{n} (ξ)

are not given explicitly and should be evaluated recursively (which becomes cumbersome for larger order of the derivative). On the other side, the use of the other procedure leads to less appealing formulae in terms of polynomials which are, however, explicitly given.

Author Contributions

Conceptualization: G.D.; methodology: G.D., S.L.; data curation: S.L.; validation: G.D., S.L., R.M.P.; formal analysis: G.D., S.L., R.M.P.; writing—original draft preparation: G.D., E.S.; writing—review and editing: S.L.

Funding

This research received no external funding.

Acknowledgments

The work of S. Licciardi was supported by an Enea Research Center individual fellowship. R.M. Pidatella wants to thank the fund of University of Catania “Metodi gruppali e umbrali per modelli di diffusione e trasporto” included in “Piano della Ricerca 2016/2018 Linea di intervento 2” for partial support of this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

Adamchik, V.S. On the Hurwitz function for rational arguments. Appl. Math. Comput. 2007, 187, 3–12. [Google Scholar] [CrossRef]
Neto, A.F. Higher Order Derivatives of Trigonometric Functions, Stirling Numbers of the Second Kind and Zeon Algebra. J. Integer Seq. 2014, 17, 14.9.3. [Google Scholar]
Feinsilver, P. Zeon algebra, Fock space, and Markov chains. Commun. Stoch. Anal. 2008, 2, 6. [Google Scholar] [CrossRef]
Cvijović, D. Derivative Polynomials and Closed-Form Higher Derivative Formulae. Appl. Math. Comput. 2009, 215, 3002–3006. [Google Scholar] [CrossRef] [Green Version]
Boyadzhiev, K.N. Derivative polynomials for tanh, tan, sech and sec in explicit form. Fibonacci Q. 2007, 45, 291–303. [Google Scholar]
Hoffman, M.E. Derivative polynomials for tangent and secant. Am. Math. Mon. 1995, 102, 23–30. [Google Scholar] [CrossRef]
Hoffman, M.E. Derivative polynomials, Euler polynomials and associated integer sequences. Electron. J. Comb. 1999, 6, R21. [Google Scholar]
Johnson, W.P. The Curious History of Faa di Bruno’s Formula. Am. Math. Mon. 2002, 109, 217–234. [Google Scholar]
Comtet, L. Advanced Combinatorics: The Art of Finite and Infinite Expansions; Springer: Berlin/Heidelberg, Germany, 1974. [Google Scholar]
Roman, S. The Umbral Calculus; Dover Publications: New York, NY, USA, 2005. [Google Scholar]
Licciardi, S. Umbral Calculus, a Different Mathematical Language. Ph.D. Thesis, Department of Mathematics and Computer Sciences, XXIX Cycle, University of Catania, Catania, Italy, 2018. [Google Scholar]
Dattoli, G.; Germano, B.; Martinelli, M.R.; Ricci, P.E. Touchard like polynomials and generalized Stirling numbers. Appl. Math. Comput. 2012, 218, 6661–6665. [Google Scholar] [CrossRef] [Green Version]
Graham, R.L.; Knuth, D.E.; Patashnik, O. Concrete Mathematics: A Foundation for Computer Science, 2nd ed.; Addison-Wesley Longman Publishing Co., Inc.: Boston, MA, USA, 1994. [Google Scholar]
Mansour, T.; Schork, M. Commutation Relations, Normal Ordering, and Stirling Numbers; Chapman and Hall/CRC: Boca Raton, FL, USA, 2015; p. 504. ISBN 9781466579880. [Google Scholar]
Dattoli, G. Generalized polynomials, operational identities and their applications. J. Comp. Appl. Math. 2000, 118, 111–123. [Google Scholar] [CrossRef] [Green Version]
Riordan, J. An Introduction to Combinatorial Analysis; Wiley: New York, NY, USA, 1980. [Google Scholar]
Bell, E.T. Exponential Polynomials. Ann. Math. 1934, 35, 258–277. [Google Scholar] [CrossRef]
Dattoli, G.; Srivastava, H.M.; Sacchetti, D. Mixed and higher-order generating functions: An extension to new families. Int. J. Math. Stat. Sci. 1999, 8, 79–87. [Google Scholar]
Babusci, D.; Dattoli, G.; Gorska, K.; Penson, K. Repeated derivatives of composite functions and generalizations of the Leibniz rule. Appl. Math. Comput. 2014, 241, 193. [Google Scholar] [CrossRef] [Green Version]
Wintucky, E.G. Formulas for nth Order Derivatives of Hyperbolic and Trigonometric Functions; NASA Lewis Research Center: Cleveland, OH, USA, 1971. [Google Scholar]
Andrews, L.C. Special Functions For Engeneers and Applied Mathematicians; Mc Millan: New York, NY, USA, 1985. [Google Scholar]
Dattoli, G.; Germano, B.; Martinelli, M.R.; Ricci, P.E. A novel theory of Legendre polynomials. Math. Comput. Modell. 2011, 54, 80–87. [Google Scholar] [CrossRef]
Dattoli, G.; Germano, B.; Martinelli, M.R.; Ricci, P.E. Sheffer and Non-Sheffer Polynomial Families. Int. J. Math. Math. Sci. 2012, 2012, 323725. [Google Scholar] [CrossRef] [Green Version]
Appél, P.; Kampé de Fériét, J. Fonctions Hypergeometriques and Hyperspheriques. Polynomes d’Hermite; Gauthiers-Villars: Paris, France, 1926. [Google Scholar]
Rainville, E.D. Special Functions; Macmillan: New York, NY, USA, 1960. [Google Scholar]
Dattoli, G.; Ottaviani, P.L.; Torre, A.; Vazquez, L. Evolution operator equations-integration with algebraic and finite-difference methods-applications to physical problems in classical and quantum mechanics and quantum field theory. La Rivista del Nuovo Cimento 1997, 20, 3–133. [Google Scholar] [CrossRef]

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Dattoli, G.; Licciardi, S.; Pidatella, R.M.; Sabia, E. Repeated Derivatives of Hyperbolic Trigonometric Functions and Associated Polynomials. Axioms 2019, 8, 138. https://doi.org/10.3390/axioms8040138

AMA Style

Dattoli G, Licciardi S, Pidatella RM, Sabia E. Repeated Derivatives of Hyperbolic Trigonometric Functions and Associated Polynomials. Axioms. 2019; 8(4):138. https://doi.org/10.3390/axioms8040138

Chicago/Turabian Style

Dattoli, Giuseppe, Silvia Licciardi, Rosa Maria Pidatella, and Elio Sabia. 2019. "Repeated Derivatives of Hyperbolic Trigonometric Functions and Associated Polynomials" Axioms 8, no. 4: 138. https://doi.org/10.3390/axioms8040138

APA Style

Dattoli, G., Licciardi, S., Pidatella, R. M., & Sabia, E. (2019). Repeated Derivatives of Hyperbolic Trigonometric Functions and Associated Polynomials. Axioms, 8(4), 138. https://doi.org/10.3390/axioms8040138

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

Repeated Derivatives of Hyperbolic Trigonometric Functions and Associated Polynomials

Abstract

1. Introduction

2. Higher Order Derivatives of Trigonometric Functions

3. Final Comments

Author Contributions

Funding

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI