Repeated Derivatives of Hyperbolic Trigonometric Functions and Associated Polynomials

: 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 ﬁeld, 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.


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 (DP) 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, . . .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 DP 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 XIX 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 ≥ 0, we set where Dx is the derivative operator d dx .They can be written in explicit form by the use of the following expression [13] 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 l m ; 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 According to Equations ( 1) and ( 2), the Touchard polynomials are explicitly given by 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.
by setting e x = ξ, we find and in conclusion we obtain Dm in which it is understood that f (r) (e x ) = Dr ξ 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 ∈ R, n ∈ N, a family of two variable polynomials loosely ascribed to the Legendre [22,23] family, then Proof.We start from the computation of the quantity which upon the use of the Laplace transform method can be written as The n-th order derivative inside the integral sign can be explicetely worked in terms of two variable Hermite polynomials [15] namely the auxiliary polynomials for the repeated derivatives of Gaussian functions according to the identity Dn Accordingly we find The proof of the identity ( 9) is readily achieved by noting that which eventually provides ∀n ∈ N 0 .
Corollary 1.If we introduce the further notation an extension of the procedure we have just envisaged allows the further results Proof.∀ν ∈ R Remark 1. Corollary 1 can be exploited to write the repeated derivatives of cos −1 (x) in the form 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 DP will be discussed in the final section of the paper.

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 Dm x (sech(x)), ∀m ∈ N 0 .
then, after setting e x = ξ, we obtain ∀m ∈ N 0 Dm The straightforward applications of the proof of the Theorem 1 yields (we omit the argument of the P n polynomials for simplicity) which, eventually, leads to We can go even further and obtain a closed form for Dm x (sech(x)) ν , ∀m ∈ N, ∀ν ∈ R.
By using the same techniques of the Proposition 1, Corollary 2 and the opportune identity sec(x) = 2 e ix +e −ix , we can move on and provide the higher order derivatives for other trigonometric circular functions.
Remark 2. The use of the Leibniz rule in Equations ( 25) or ( 26) can be avoided by setting and eventually ending up with where δ m,0 is Kronecker's delta, indeed In the forthcoming section we will discuss the comparison with previous papers and present possible developments along the lines we have indicated.

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 DP 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 with P n (z) satisfying the generating function and 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 Dn x (tan(x)) = Π n (tan(x)) with ∏ 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 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 t n n!
, by summing up over the index n and ending up with The use of the Lie derivative identity [26] e t[(1+ξ which is the generating function given in [4,5].
The second family of DP can be defined by the use of the same procedure.According to refs.[4,5] they are implicitly defined by the condition

Dn
x (sec(x)) = sec(x)Q n (tan(x)) which, by the use of the same change of variable leading to Equation (33), yields which can be straightforwardly exploited to derive the relevant properties.The use of the identity (35) eventually yields the associated generating functions which is the same reported in refs.[4,5].By recalling that the Hoppe formula writes for the case of sec(x), we find 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 DP to evaluate the derivatives of the type λ m,j (x) = Dm x (cos j (x)).
Making the ansatz that Dm x (cos j (x)) = (−1) j Λ m,j (tan(x)) and, by using the same procedure leading to Equation (33), we can identify the function which is specified by the generating function It is worth to consider the alternative assumption The two methods we have discussed in this paper, namely the procedure based on the DP 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.