y=gdx=arcsin(tanhx)y = \operatorname{gd} x = \arcsin(\tanh x)

Graph of the Gudermannian y=gdxy = \operatorname{gd} x

y=gdx=arcsin(tanhx)y = \operatorname{gd} x = \arcsin(\tanh x) is called the Gudermannian function, and it ties the trigonometric and the hyperbolic functions together without the help of complex numbers1. It is named after the nineteenth-century German mathematician Gudermann, and it is also the latitude itself in the Mercator map projection.

Domain and range

  • The domain is all real numbers
  • The range is π2<y<π2-\dfrac{\pi}{2} < y < \dfrac{\pi}{2}
  • It increases monotonically
  • It is an odd function

The range of tanhx\tanh x is (1,1)(-1, 1), so the argument of arcsin\arcsin always falls within its domain. The lines y=±π2y = \pm\dfrac{\pi}{2} are horizontal asymptotes; the values close in on them without ever arriving.

Symmetry and monotonicity

Both tanh\tanh and arcsin\arcsin are odd, so their composition is odd and has point symmetry about the origin. The derivative simplifies remarkably.

ddxarcsin(tanhx)=11tanh2xsech2x=sech2xsechx=sechx\begin{align*} \frac{d}{dx}\arcsin(\tanh x) &= \frac{1}{\sqrt{1 - \tanh^2 x}} \cdot \operatorname{sech}^2 x \\ &= \frac{\operatorname{sech}^2 x}{\operatorname{sech} x} = \operatorname{sech} x \end{align*}

So gd(x)=sechx\operatorname{gd}'(x) = \operatorname{sech} x, and the function may equally be defined by an integral.

gdx=0xsechtdt\operatorname{gd} x = \int_0^{x} \operatorname{sech} t\,dt

Since sechx\operatorname{sech} x is always positive the function increases monotonically. The slope at the origin is sech0=1\operatorname{sech} 0 = 1, so the graph is tangent there to the line y=xy = x.

Concavity

The second derivative is y=sechxtanhxy'' = -\operatorname{sech} x \tanh x, whose sign is the opposite of the sign of xx. The curve is concave up for x<0x < 0 and concave down for x>0x > 0, with the origin as its only inflection point.

A bridge between the two families

The point of this function is a set of identities that hold once one writes θ=gdx\theta = \operatorname{gd} x.

On the trigonometric sideOn the hyperbolic side
sinθ\sin\thetatanhx\tanh x
cosθ\cos\thetasechx\operatorname{sech} x
tanθ\tan\thetasinhx\sinh x
secθ\sec\thetacoshx\cosh x

All of them hold at once. The function is a dictionary that translates hyperbolic functions into trigonometric ones with no recourse to complex numbers.

The inverse and the integral of the secant

The inverse is gd1(y)=artanh(siny)\operatorname{gd}^{-1}(y) = \operatorname{artanh}(\sin y), which can also be written as follows.

gd1(y)=ln(secy+tany)\operatorname{gd}^{-1}(y) = \ln(\sec y + \tan y)

The right-hand side is nothing but the main part of the integral secydy=lnsecy+tany+C\int \sec y\,dy = \ln|\sec y + \tan y| + C. The formula for the integral of the secant, which looks so abrupt in a table of integrals, turns out to be this function.

The Mercator projection

In the Mercator projection a place at latitude φ\varphi is drawn at the vertical coordinate 0φsectdt\int_0^{\varphi} \sec t\,dt2. That is gd1(φ)\operatorname{gd}^{-1}(\varphi), so conversely the latitude corresponding to a vertical coordinate xx on the map is gdx\operatorname{gd} x.

xx, the vertical coordinategdx\operatorname{gd} x in radiansLatitude
000000^\circ
110.8657\approx 0.865749.60\approx 49.60^\circ
221.3018\approx 1.301874.58\approx 74.58^\circ
π\pi1.4842\approx 1.484285.05\approx 85.05^\circ

Being tangent to y=xy = x near the origin corresponds to the small distortion near the equator, and approaching ±π2\pm\dfrac{\pi}{2} at the ends corresponds to the poles being banished to infinity. Web Mercator, used widely by map services, cuts the latitude off at ±85.05\pm 85.05^\circ because it cuts the map at x=±πx = \pm\pi to make it square.

  1. Gudermannian function, Wikipedia
  2. Mercator projection, Wikipedia