Article 778DA exp_q

exp_q

by
John
from John D. Cook on (#778DA)

The function expq(x) is defined by taking the power series for exp(x) and keeping only the terms whose index is a multiple of q. For example, exp2(x) keeps only the even-numbered terms in the exponential power series and so equals cosh(x).

expqx.svg

In general,

expq2.svg

The first sum uses Iverson's bracket notation: a Boolean expression in brackets denotes the function that returns 1 when the expression is true and zero when it is false. Here the bracket equals 1 when q divides n and is zero otherwise.

Closed forms

Let = exp(2i / q). Then

expq3.svg

This lets us find closed-form expressions for expq(x). For example, when q = 4, = i and

expq4.svg

Here's a proof of the identity above:

expq5.svg

In the proof we used the identity

expq6.svg

which is important in deriving the properties of the discrete Fourier transform.

Differential equations

The first time I saw the function expq(x) was in differential equations, though I didn't know at the time the function had a name.

When a course in differential equations gets to power series solutions, a common example or homework problem is to solve

expq7.svg

for k = 3 or 4, i.e. to find a function that equals its third or fourth derivative.

If the initial conditions are

expq8.svg

and

expq9.svg

the unique solution to

expq7.svg

is y(x) = expk(x).

Mathematica and Mittag-Leffler

Mathematica does not have a built-in function implementing expq(x), but it does have an implementation of the Mittag-Leffler function, and so thanks to a relation between this function and expq(x) you can implement the latter as

expq[x_, q_] := MittagLefflerE[q, x^q]
Combinatorics

The first time I saw the notation expq(x) was in combinatorics. I had intended to include an application from that book here, but I make that the topic for the next post.

The post exp_q first appeared on John D. Cook.
External Content
Source RSS or Atom Feed
Feed Location http://feeds.feedburner.com/TheEndeavour?format=xml
Feed Title John D. Cook
Feed Link https://www.johndcook.com/blog
Reply 0 comments