y=gdx=arcsin(tanhx) 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π
It increases monotonically
It is an odd function
The range of tanhx is (−1,1), so the argument of arcsin always falls within its domain. The lines y=±2π are horizontal asymptotes; the values close in on them without ever arriving.
Symmetry and monotonicity
Both tanh and arcsin are odd, so their composition is odd and has point symmetry about the origin. The derivative simplifies remarkably.
So gd′(x)=sechx, and the function may equally be defined by an integral.
gdx=∫0xsechtdt
Since sechx is always positive the function increases monotonically. The slope at the origin is sech0=1, so the graph is tangent there to the line y=x.
Concavity
The second derivative is y′′=−sechxtanhx, whose sign is the opposite of the sign of x. The curve is concave up for x<0 and concave down for x>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.
On the trigonometric side
On the hyperbolic side
sinθ
tanhx
cosθ
sechx
tanθ
sinhx
secθ
coshx
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 gd−1(y)=artanh(siny), which can also be written as follows.
gd−1(y)=ln(secy+tany)
The right-hand side is nothing but the main part of the integral ∫secydy=ln∣secy+tany∣+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 φ is drawn at the vertical coordinate ∫0φsectdt2. That is gd−1(φ), so conversely the latitude corresponding to a vertical coordinate x on the map is gdx.
x, the vertical coordinate
gdx in radians
Latitude
0
0
0∘
1
≈0.8657
≈49.60∘
2
≈1.3018
≈74.58∘
π
≈1.4842
≈85.05∘
Being tangent to y=x near the origin corresponds to the small distortion near the equator, and approaching ±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∘ because it cuts the map at x=±π to make it square.