Fact-checked by Grok 2 weeks ago

Gudermannian function

The Gudermannian function, denoted \gd(x) or \gamma(x), is an odd special function in mathematics that provides a direct relationship between hyperbolic and trigonometric functions without invoking complex numbers, defined as \gd(x) = \int_0^x \sech t \, dt for -\infty < x < \infty. It can also be expressed in closed form as \gd(x) = 2 \arctan(\tanh(x/2)), \gd(x) = \arctan(\sinh x), or \gd(x) = 2 \arctan(e^x) - \pi/2. Named after the German mathematician Christoph Gudermann (1798–1852), who explored connections between circular and hyperbolic functions in publications starting in 1830, the function serves as a bridge between these two classes of functions through identities such as \sin(\gd(x)) = \tanh x, \cos(\gd(x)) = \sech x, and \tan(\gd(x)) = \sinh x. Its derivative is \gd'(x) = \sech x, and it maps the real line to (-\pi/2, \pi/2), making it strictly increasing and bijective. One of its most notable applications is in cartography, where it appears in the inverse Mercator projection to relate latitude \phi to the vertical map coordinate y via \phi = \gd(y), facilitating the conformal representation of the Earth's surface on a cylinder. The inverse Gudermannian function, \gd^{-1}(x) = \int_0^x \sec t \, dt for -\pi/2 < x < \pi/2, satisfies \gd^{-1}(x) = \ln(\sec x + \tan x) and extends these relations symmetrically. Additionally, the function connects to other special functions, such as elliptic integrals and the dilogarithm in its antiderivative, and features a Maclaurin series expansion x - \frac{1}{6}x^3 + \frac{1}{24}x^5 - \frac{61}{5040}x^7 + \cdots.

Definition and Basic Properties

Definition

The Gudermannian function, often denoted as \mathrm{gd}(x) or \gamma(x), is an odd real-valued function that provides a bridge between hyperbolic and trigonometric measures without invoking complex numbers. Its primary definition is given by the definite integral \mathrm{gd}(x) = \int_0^x \sech(t) \, dt, where \sech(t) = 1 / \cosh(t). This integral form originates from the accumulated arc length or area considerations in hyperbolic geometry, reflecting the function's role in mapping hyperbolic parameters to angular ones. Equivalent closed-form expressions for the Gudermannian function include \mathrm{gd}(x) = \arctan(\sinh x) and \mathrm{gd}(x) = 2 \arctan(e^x) - \frac{\pi}{2}. These representations facilitate direct computation and reveal the function's smooth, monotonically increasing behavior from hyperbolic inputs to bounded angular outputs. Geometrically, the Gudermannian function interprets \mathrm{gd}(x) as the circular angle \phi corresponding to a hyperbolic angle \psi = x under a stereographic projection, which equates the areas of associated circular and hyperbolic sectors while preserving angles between the two geometries. The function is defined for all real x \in \mathbb{R}, with its range being the open interval (- \pi/2, \pi/2), ensuring it maps the entire real line bijectively onto this angular domain.

Inverse Gudermannian function

The inverse Gudermannian function, denoted \operatorname{gd}^{-1}(x) or \operatorname{arcgd}(x), is the inverse of the Gudermannian function \operatorname{gd}(x), providing a bijection from the open interval (-\pi/2, \pi/2) to \mathbb{R}. It arises naturally in contexts bridging trigonometric and hyperbolic functions, such as in the where it relates latitude \phi to the vertical coordinate y. One standard definition is given by the integral \operatorname{gd}^{-1}(x) = \int_0^x \sec t \, dt, valid for -\pi/2 < x < \pi/2. This form emphasizes its role as an area function under the secant curve. Equivalent closed-form expressions include \operatorname{gd}^{-1}(x) = \sinh^{-1}(\tan x) = \tanh^{-1}(\sin x), \operatorname{gd}^{-1}(x) = \frac{1}{2} \ln \left( \frac{1 + \sin x}{1 - \sin x} \right), and \operatorname{gd}^{-1}(x) = \ln (\sec x + \tan x) = \ln \left[ \tan \left( \frac{\pi}{4} + \frac{x}{2} \right) \right]. These identities highlight its connections to inverse hyperbolic functions and logarithmic forms, facilitating computations in hyperbolic geometry and projections. The derivative of the inverse Gudermannian function is \frac{d}{dx} \operatorname{gd}^{-1}(x) = \sec x, which follows directly from the fundamental theorem of calculus applied to its integral definition. In applications, it is particularly valuable for converting between spherical latitudes and Cartesian coordinates in cartography; for instance, in the Mercator projection, y = \operatorname{gd}^{-1}(\phi) yields the northward distance from the equator for a given latitude \phi. The function's Maclaurin series expansion is \operatorname{gd}^{-1}(x) = x + \frac{1}{6} x^3 + \frac{1}{24} x^5 + \frac{61}{5040} x^7 + O(x^9), useful for small-angle approximations in numerical evaluations.

Fundamental properties

The Gudermannian function \operatorname{gd}(x) is defined for all real numbers x \in \mathbb{R} and is an odd function, satisfying \operatorname{gd}(-x) = -\operatorname{gd}(x) for all x. Its range is the open interval (-\pi/2, \pi/2), and it provides a bijection from \mathbb{R} onto this interval. The function is continuous and strictly increasing on \mathbb{R}, with \operatorname{gd}(0) = 0. This monotonicity follows from the integral definition \operatorname{gd}(x) = \int_0^x \operatorname{sech} t \, \mathrm{d}t, as the integrand \operatorname{sech} t > 0 for all real t. As x \to \infty, \operatorname{gd}(x) \to \pi/2, and as x \to -\infty, \operatorname{gd}(x) \to -\pi/2. These horizontal asymptotes reflect the bounded nature of the function despite its unbounded domain. The Gudermannian function is infinitely differentiable, with its expansion around zero given by \operatorname{gd}(x) = x - \frac{1}{6}x^3 + \frac{1}{24}x^5 - \frac{61}{5040}x^7 + \cdots, converging for all real x.

Relations to Trigonometric and Hyperbolic Functions

Circular-hyperbolic identities

The Gudermannian function \gd(x) provides a bridge between circular and through a set of fundamental identities that express of \gd(x) directly in terms of of x. These relations arise from the definition \gd(x) = \int_0^x \sech t \, \mathrm{d}t and its equivalent expressions using . Specifically, the sine identity is \sin(\gd(x)) = \tanh x, which follows from \gd(x) = \arcsin(\tanh x). Similarly, \cos(\gd(x)) = \sech x, derived from \gd(x) = \arccos(\sech x). The tangent relation is \tan(\gd(x)) = \sinh x, consistent with \gd(x) = \arctan(\sinh x). These identities hold for x \in (-\infty, \infty), with \gd(x) \in (-\pi/2, \pi/2). Additional identities extend to reciprocal functions: \cot(\gd(x)) = \csch x, \sec(\gd(x)) = \cosh x, and \csc(\gd(x)) = \coth x. These can be obtained by taking reciprocals of the primary identities or using the expressions \gd(x) = \arccsc(\coth x) and \gd(x) = \arcsec(\cosh x). Such relations facilitate conversions between hyperbolic and spherical geometries, notably in cartographic projections. A complementary set of identities involves half-arguments: \tanh(x/2) = \tan(\gd(x)/2). Furthermore, \gd(x) = 2 \arctan(e^x) - \pi/2, linking it to the and underscoring its role in unifying circular and hyperbolic . These identities are essential for applications in and special function theory.

Symmetries and periodicity

The Gudermannian function \operatorname{gd}(x) is an odd function, satisfying \operatorname{gd}(-x) = -\operatorname{gd}(x) for all real x, which implies point symmetry (or of 180 degrees) about the . This property follows directly from its integral definition \operatorname{gd}(x) = \int_0^x \operatorname{sech} t \, dt, as the integrand \operatorname{sech} t is even. In the complex domain, the function exhibits additional symmetry relations, such as \operatorname{gd}(ix) = i \operatorname{gd}^{-1}(x), linking it to its and highlighting its role in bridging trigonometric and without invoking complex exponentials explicitly. The Gudermannian function is not periodic over the real numbers, as it is strictly increasing and bounded, with \lim_{x \to \infty} \operatorname{gd}(x) = \pi/2 and \lim_{x \to -\infty} \operatorname{gd}(x) = -\pi/2, approaching horizontal asymptotes without repetition. No quasi-periodic or higher-order periodic behaviors are observed in its real-valued form, consistent with its monotonic nature derived from the positive derivative \operatorname{gd}'(x) = \operatorname{sech} x > 0.

Evaluation Methods

Specific values

The Gudermannian function satisfies \operatorname{gd}(0) = 0, as follows directly from its definition or any of the equivalent closed-form expressions. It is an odd function, with \operatorname{gd}(-x) = -\operatorname{gd}(x) for all real x, a property inherited from the odd nature of \sinh x and \arctan y. The function maps the real line to the open interval (-\pi/2, \pi/2), approaching these bounds asymptotically: \lim_{x \to \infty} \operatorname{gd}(x) = \pi/2 and \lim_{x \to -\infty} \operatorname{gd}(x) = -\pi/2. The identity \operatorname{gd}(x) = \arctan(\sinh x) enables exact evaluation at points where \sinh x = \tan \theta for known angles \theta \in (-\pi/2, \pi/2). For example:
  • At x = \sinh^{-1}(1/\sqrt{3}) = \frac{1}{2} \ln 3, \operatorname{[gd](/page/GD)}(x) = \arctan(1/\sqrt{3}) = \pi/6.
  • At x = \sinh^{-1}(1) = \ln(1 + \sqrt{2}), \operatorname{[gd](/page/GD)}(x) = \arctan(1) = \pi/4.
  • At x = \sinh^{-1}(\sqrt{3}) = \ln(2 + \sqrt{3}), \operatorname{[gd](/page/GD)}(x) = \arctan(\sqrt{3}) = \pi/3.
These evaluations highlight the function's role in bridging and circular geometries, with arguments expressible via the standard formula \sinh^{-1} y = \ln(y + \sqrt{y^2 + 1}).

Taylor series expansion

The expansion of the Gudermannian function \mathrm{gd}(x) about x=0 is an infinite series consisting solely of odd powers of x, reflecting its odd nature. This expansion is derived by term-by-term of the Maclaurin series for \mathrm{gd}'(x) = \sech x, which itself involves the Euler numbers E_n. The general form is \mathrm{gd}(x) = \sum_{k=0}^{\infty} (-1)^k \frac{E_{2k}}{(2k+1)!} x^{2k+1}, where the Euler numbers E_{2k} (with E_0 = 1, E_2 = 1, E_4 = 5, E_6 = 61, E_8 = 1385, and so on) are the absolute values of the secant (or zigzag) Euler numbers, which appear in the expansion of \sech x. The first several terms of the series are \mathrm{gd}(x) = x - \frac{1}{6} x^3 + \frac{5}{120} x^5 - \frac{61}{5040} x^7 + \frac{1385}{362880} x^9 - \cdots, or equivalently, \mathrm{gd}(x) = x - \frac{1}{6} x^3 + \frac{1}{24} x^5 - \frac{61}{5040} x^7 + \frac{277}{72576} x^9 - \cdots. These coefficients correspond directly to the integrated form of the \sech x series \sech x = \sum_{k=0}^{\infty} (-1)^k \frac{E_{2k}}{(2k)!} x^{2k}. The series converges for |x| < \pi/2, limited by the radius of convergence of the \sech x expansion, beyond which poles in the complex plane (at x = i(\pi/2 + k\pi) for integer k) affect analytic continuation. This expansion provides a useful approximation for small |x|, connecting the Gudermannian function to classical special function series.

Numerical computation

The Gudermannian function \gd(x) is typically computed using closed-form expressions involving elementary functions, which are highly efficient and supported in standard numerical libraries such as those in Python's NumPy, MATLAB, or C++'s . One primary expression is \gd(x) = \arctan(\sinh x), which directly relates the hyperbolic sine to the inverse tangent and avoids explicit integration. This form is suitable for moderate values of x, where \sinh x remains within the representable range of floating-point arithmetic, typically up to |x| \approx 710 in double precision to prevent overflow. For improved , especially for large |x|, the equivalent expression \gd(x) = 2 \arctan(\tanh(x/2)) is preferred, as \tanh(x/2) saturates to \pm 1 without , ensuring accurate evaluation near the asymptotic limits \gd(x) \to \pm \pi/2. Another variant, \gd(x) = 2 \arctan(e^x) - \pi/2, appears in classical references but can introduce for large positive x due to e^x, making it less ideal for general-purpose computation. These expressions exploit the bounded nature of \arctan and the saturation of , minimizing cancellation errors and achieving full machine precision for most inputs. Although the defining integral \gd(x) = \int_0^x \sech t \, dt provides a conceptual basis, direct numerical quadrature (e.g., via Gaussian methods or adaptive ) is rarely used in practice due to higher computational cost compared to the closed forms; it may be relevant only for custom high-precision needs or validation. In multiprecision environments, such as those using , the closed forms extend naturally to evaluate \gd(x) to hundreds of decimal places by applying series expansions or iterative refinements to the component functions if needed. Implementations in libraries like John Burkardt's POLPAK collection follow the $2 \arctan(\tanh(x/2)) approach for robustness across real arguments.

Calculus

Derivatives

The first derivative of the Gudermannian function \gd(x) is \frac{d}{dx} \gd(x) = \sech x. This follows directly from its integral definition \gd(x) = \int_0^x \sech t \, dt. For the inverse Gudermannian function \gd^{-1}(\phi), the derivative is \frac{d}{d\phi} \gd^{-1}(\phi) = \sec \phi. This arises from the integral representation \gd^{-1}(\phi) = \int_0^\phi \sec t \, dt. The form of the first derivative \sech x provides an interpretive link to the function's integral definition, where \sech x serves as the integrand in the context of area relations between hyperbolic and circular sectors. The second derivative is \gd''(x) = -\sech x \tanh x. This expression is obtained by differentiating \sech x using standard hyperbolic differentiation rules. Higher-order derivatives of \gd(x) can be derived explicitly through repeated application of hyperbolic identities to the derivatives of \sech x, or recursively using general formulas for the nth derivatives of \sech x. Applications of the chain rule to compositions, such as \gd(f(x)), yield the \sech(f(x)) f'(x), facilitating analysis in contexts like differential equations or parametric relations between trigonometric and .

Integral representations

The Gudermannian function \gd(x) is fundamentally defined by its integral representation as the of the : \gd(x) = \int_0^x \sech t \, dt = \int_0^x \frac{dt}{\cosh t}, where -\infty < x < \infty and \gd(x) maps to (-\pi/2, \pi/2). This form arises from the geometric interpretation connecting and circular sectors through , where the measures the accumulated "" corresponding to the hyperbolic parameter x. This can be evaluated through substitutions that yield closed-form expressions in terms of . For instance, substituting u = \sinh t transforms the integral to \gd(x) = \arctan(\sinh x), while u = \tanh(t/2) leads to \gd(x) = 2 \arctan(\tanh(x/2)). These derivations preserve the integral's validity across the real line and extend analytically to the , excluding branch cuts along the imaginary axis where convergence issues arise. The inverse Gudermannian function \gd^{-1}(\phi) also possesses a dual integral representation: \gd^{-1}(\phi) = \int_0^\phi \sec u \, du = \int_0^\phi \frac{du}{\cos u}, for -\pi/2 < \phi < \pi/2, mirroring the structure of \gd(x) but in the trigonometric domain. This symmetry highlights the bidirectional mapping between hyperbolic and circular measures.

Advanced Identities

Argument-addition formulas

The argument-addition formula for the Gudermannian function \mathrm{gd}(x) expresses \mathrm{gd}(x + y) in terms of trigonometric functions of \mathrm{gd}(x) and \mathrm{gd}(y). Using the definition \mathrm{gd}(x) = \arctan(\sinh x), the formula follows from the hyperbolic addition theorem \sinh(x + y) = \sinh x \cosh y + \cosh x \sinh y and the identities \sinh z = \tan(\mathrm{gd}(z)), \cosh z = \sec(\mathrm{gd}(z)) for z \in \mathbb{R}. Substituting yields \sinh(x + y) = \frac{\sin(\mathrm{gd}(x))}{\cos(\mathrm{gd}(x)) \cos(\mathrm{gd}(y))} + \frac{\sin(\mathrm{gd}(y))}{\cos(\mathrm{gd}(x)) \cos(\mathrm{gd}(y))} = \frac{\sin(\mathrm{gd}(x)) + \sin(\mathrm{gd}(y))}{\cos(\mathrm{gd}(x)) \cos(\mathrm{gd}(y))}. Thus, \mathrm{gd}(x + y) = \arctan\left( \frac{\sin(\mathrm{gd}(x)) + \sin(\mathrm{gd}(y))}{\cos(\mathrm{gd}(x)) \cos(\mathrm{gd}(y))} \right). This identity bridges the additive structure of hyperbolic arguments with the corresponding circular angles via the Gudermannian mapping. An alternative form arises from the representation \mathrm{gd}(x) = 2 \arctan(\tanh(x/2)). The addition theorem for the tangent gives \tanh\left(\frac{x + y}{2}\right) = \frac{\tanh(x/2) + \tanh(y/2)}{1 + \tanh(x/2) \tanh(y/2)}, so \mathrm{gd}(x + y) = 2 \arctan\left( \frac{\tanh(x/2) + \tanh(y/2)}{1 + \tanh(x/2) \tanh(y/2)} \right). This expression highlights the connection to the addition formula, adjusted by the sign in the denominator due to the hyperbolic nature.

Complex extension

The Gudermannian function admits a natural extension to arguments via its defining : \gd(z) = \int_0^z \sech t \, \mathrm{d}t, where the path of integration lies within a simply connected avoiding the poles of \sech t at t = i\left(\frac{\pi}{2} + k\pi\right) for each k. This provides the from the real line to the , rendering \gd(z) holomorphic in regions excluding these singularities. The presence of poles in the integrand implies that \gd(z) is multi-valued, with branch points precisely at these locations; a principal can be defined by introducing branch cuts, typically extending vertically along the imaginary axis from each pole to i\infty or -\i\infty, ensuring continuity from the right across the cuts. Equivalent closed-form expressions valid on the real line extend by to the domain, inheriting structures from the constituent functions: \gd(z) = 2\arctan(e^z) - \frac{\pi}{2}, \gd(z) = \arctan(\sinh z), \gd(z) = \arcsin(\tanh z), \gd(z) = 2\arctan\left(\tanh\frac{z}{2}\right). The principal of \arctan w (with w \in \mathbb{C}) features branch cuts along the rays \{iy : y \geq 1\} and \{iy : y \leq -1\} on the imaginary , while \arcsin w has cuts along [-\infty, -1] and [1, \infty) on the real axis. Consequently, the branch cuts of \gd(z) arise where the inner functions (e.g., e^z, \sinh z) map to these cut loci, resulting in a network of cuts aligned with the poles of \sech t. For instance, in the expression involving e^z, discontinuities arise when e^z is purely imaginary with |\operatorname{Im}(e^z)| \geq 1, corresponding to half-lines in the z-plane parallel to the real axis. These extensions preserve key identities, such as \gd(\mathrm{i}w) = \mathrm{i} \gd^{-1}(w) for appropriate branches, linking the function to its inverse in the complex setting.

Historical Development

Origins and early uses

The Gudermannian function, which relates and trigonometric functions, has roots in 16th-century . In 1569, developed his conformal cylindrical , where the vertical coordinate for \phi is given by the \int_0^\phi \sec t \, dt, a form later recognized as the Gudermannian function \mathrm{gd}^{-1}(x) = 2 \arctan(\tanh(x/2)). This ensured constant scale along meridians, essential for , though Mercator computed it numerically without explicit hyperbolic connections. The function was formally introduced in the 1760s by , who coined the term "transcendent angle" for the relation between circular and angles in his work on . Lambert's formulation arose from stereographic projections linking circular sectors to ones, providing a bridge without complex numbers. This appeared in his logarithmic tables and treatises on , where it facilitated computations in and astronomy. Christoph Gudermann advanced the function's study in the 1830s through a series of papers on elliptic and hyperbolic integrals. In his 1830 memoir "Theorie der Potenzial- oder Cyklisch-Hyperbolischen Functionen" published in Crelle's Journal (Vol. 6, p. 165), he detailed the "longitude of u"—an early designation for the Gudermannian—as part of cyclic-hyperbolic functions, including extensive tables for practical evaluation. Gudermann's work emphasized analogies between trigonometric and hyperbolic systems, influencing later developments in special functions.

Naming and formalization

The Gudermannian function originated in the work of during the 1760s, as part of his foundational contributions to . Lambert introduced the concept as the "transcendent angle," defined implicitly through the integral relation connecting the arc length on a rectangular to angular measures, without relying on complex numbers. This served to bridge the geometric interpretations of circular and hyperbolic sectors via , laying the groundwork for its role in relating trigonometric and hyperbolic identities. In the 1830s, Christoph Gudermann significantly advanced the theory through a series of papers, most notably his 1830 memoir "Theorie der Potenzial- oder Cyklisch-Hyperbolischen Functionen," published in Crelle's Journal (Journal für die reine und angewandte Mathematik), volume 6. Gudermann formalized the function as the "longitude of u," denoted \lambda u or l u, and expressed it as the inverse of a hyperbolic operation: specifically, \lambda u = \int_0^u \sech t \, dt, which equals \arctan(\sinh u). He derived extensive identities, such as \tan(\lambda u) = \sinh u and \sec(\lambda u) = \cosh u, emphasizing its utility in expressing in terms of circular ones and vice versa. Gudermann also computed detailed tables and explored its differential properties, establishing it as a key tool in the analytic theory of . The modern name "Gudermannian function," often abbreviated as \gd(x), was proposed in 1862 by to honor Gudermann's pioneering efforts. In his article "On the Elliptic and " published in the , Cayley adopted the term for the function previously known by various designations, including Lambert's transcendent angle and Gudermann's longitude. This naming formalized its recognition as a distinct special function in English mathematical literature, with the standard definition \gd(x) = 2 \arctan(\tanh(x/2)) emerging from these historical developments to ensure consistency in subsequent works.

Generalizations and Extensions

Multidimensional generalizations

The Gudermannian function, originally defined for one-dimensional relations between hyperbolic and , extends to higher dimensions primarily through its role in coordinate systems for spaces and in modern approximation theory. In the context of , the function facilitates conformal representations of n-dimensional \mathbb{H}^n via stereographic projections from the to the sphere. Specifically, in polar-like coordinates, the hyperbolic radial distance r is related to an angular coordinate \theta by \theta = \gd(r), where \gd(r) = \int_0^r \sech t \, dt = \arctan(\sinh r). This yields the ds^2 = \sec^2 \theta \, (d\theta^2 + \sin^2 \theta \, d\Omega_{n-1}^2), where d\Omega_{n-1}^2 is the standard metric on the unit (n-1)-sphere. This form generalizes the 2D case, preserving conformality and enabling visualization and computation in higher dimensions, as seen in projections for solving the in curved spaces. In determination and parameterizations, the Gudermannian function generalizes to higher-dimensional map projections, particularly for 4D representations of . Here, the inverse Gudermannian \gd^{-1}(\phi) = \arsinh(\tan \phi) relates spherical latitudes \phi to hyperbolic parameters, extended via higher-order Mercator parameters m\mu defined as r(\phi) = 2 \arctanh(\tan(\phi/2)) for order m \geq 1. This construction maps angles to a parameter space with reduced (e.g., |\phi| < \pi/(2m)), linking to higher-order parameters through \tanh(2r_m) = r_{m(\mu)} and avoiding singularities in multi-dimensional groups SO(n). Such generalizations are crucial for modeling, where they provide periodic, unbounded representations analogous to Mercator projections on spheres. A distinct multidimensional extension arises in approximation theory, where the multivariate Gudermannian function serves as a density-inducing sigmoid for neural network operators in \mathbb{R}^n (n \geq 2). Defined via the standard Gudermannian \gd(x) = 2 \arctan(\tanh(x/2)) extended to a multidimensional kernel, it generates normalized, quasi-interpolation, Kantorovich, and quadrature operators for approximating continuous Banach space-valued functions on compact boxes or \mathbb{R}^n. These operators achieve pointwise and uniform convergence, with Jackson-type error estimates bounded by the multivariate modulus of continuity or Fréchet derivatives, e.g., \|L_{m,n}(f; X) - f(X)\| \leq C \omega_2(f; \delta/m^{1/n}) for second-order modulus \omega_2. This framework leverages the bounded, monotonic properties of \gd for multivariate sigmoid activation in one-hidden-layer networks, outperforming logistic sigmoids in high-dimensional function approximation. The Gudermannian function, denoted \operatorname{gd}(x), serves as a bridge between trigonometric and , expressing key identities that link their behaviors without invoking complex numbers. Specifically, \sin(\operatorname{gd} x) = \tanh x, \cos(\operatorname{gd} x) = \sech x, and \tan(\operatorname{gd} x) = \sinh x. These relations highlight its role in connecting elementary circular and , which, while not "special" in the non-elementary sense, form the foundational ties for the Gudermannian's applications in projections and . A primary connection exists with the Jacobi amplitude function \operatorname{am}(u, k), where the Gudermannian emerges as the limiting case when the elliptic modulus k = 1: \operatorname{am}(x, 1) = \operatorname{gd} x. In this degeneration, the Jacobi elliptic functions \operatorname{sn}(u, k), \operatorname{cn}(u, k), and \operatorname{dn}(u, k) reduce to hyperbolic functions: \operatorname{sn}(x, 1) = \tanh x, \operatorname{cn}(x, 1) = \sech x, and \operatorname{dn}(x, 1) = \sech x. This link underscores the Gudermannian's position as the "hyperbolic amplitude," facilitating transitions from elliptic to hyperbolic regimes in integrals and inverses. Approximations for small complementary modulus k' further relate \operatorname{am}(x, k) to \operatorname{gd} x via expansions involving hyperbolic terms, such as \operatorname{am}(x, k) = \operatorname{gd} x - \frac{1}{4} k'^2 (x - \sinh x \cosh x) \sech x + O(k'^4). The inverse Gudermannian function \operatorname{gd}^{-1}(\phi) is tied to the : \operatorname{gd}^{-1}(\phi) = \int_0^\phi \sec t \, dt, which represents the non-elementary secant integral. This connection appears in contexts like the , where latitude relates to hyperbolic distances. Additionally, the indefinite integral of the Gudermannian itself involves the \operatorname{Li}_2(z), providing a pathway to polylogarithmic in analytic continuations and series expansions.

Applications

Cartography and projections

The Gudermannian function plays a central role in the , a conformal cylindrical developed for navigation purposes. In this projection, geographic coordinates of \phi and \lambda are transformed to Cartesian coordinates (x, y) on a plane, where x = R \lambda (with R as the radius of the sphere) and y = R \cdot \mathrm{gd}^{-1}(\phi), with the inverse Gudermannian function \mathrm{gd}^{-1}(\phi) = \ln\left[\tan\left(\frac{\pi}{4} + \frac{\phi}{2}\right)\right]. This logarithmic scaling ensures that meridians and parallels are represented as equally spaced vertical and horizontal lines, respectively, preserving angles and making rhumb lines (lines of constant bearing) appear as straight lines, which is essential for maritime navigation. The , which converts the coordinates back to , directly employs the Gudermannian function: \phi = \mathrm{gd}(y/R), where \mathrm{gd}(x) = \int_0^x \sech(t) \, dt = 2 \arctan(\tanh(x/2)) = \arctan(\sinh x). This bridges of (circular on the ) with , reflecting the projection's conformal property derived from the stereographic mapping between spherical and hyperbolic geometries. The function's integral form arises naturally from the for preserving local shapes, as the scale factor along meridians must equal that along parallels. In modern cartography, the Gudermannian function extends to variants like the Web Mercator projection, used in digital mapping systems such as those employed by and . Here, the projection is adapted to the Pseudo-Mercator coordinate reference system (EPSG:3857), normalizing coordinates to a for web tiling: y(\phi) = \frac{\mathrm{gd}^{-1}(\phi)}{2\pi} + \frac{1}{2}, with latitude typically clipped to \pm 85.0511^\circ to avoid singularities at the poles. This application maintains the conformal distortion characteristics of the original Mercator while facilitating efficient rendering of large-scale interactive maps, though it introduces significant areal distortion in polar regions. The function's computational efficiency, via closed-form expressions like \mathrm{gd}(x) = \arctan(\sinh x), supports real-time geotransformations in geographic information systems (GIS). Beyond the standard Mercator, the Gudermannian function influences related projections, such as the transverse Mercator used in the grid system for large-scale military and maps. In these, hyperbolic elements analogous to the Gudermannian appear in series expansions for high-precision forward and inverse transformations, ensuring minimal distortion along a central . For instance, the latitude-to-y mapping incorporates terms derived from \sech integrals to achieve sub-meter accuracy over zones spanning 6° of . This underscores the function's broader utility in maintaining conformality across diverse cartographic frameworks.

Hyperbolic geometry

The Gudermannian function, denoted \gd(x), plays a pivotal role in by linking circular to their hyperbolic counterparts, enabling the translation of geometric concepts between and non-Euclidean settings. Defined as \gd(x) = \int_0^x \sech t \, dt = \arctan(\sinh x), it satisfies fundamental identities such as \sin(\gd(x)) = \tanh x, \cos(\gd(x)) = \sech x, and \tan(\gd(x)) = \sinh x. These relations stem from the geometric interpretation where \gd(x) represents the angle in a with opposite side \sinh x, adjacent side 1, and hypotenuse \cosh x, allowing hyperbolic lengths and angles to be expressed trigonometrically without complex arguments. A key application arises in the angle of parallelism, a characteristic feature of that describes the angle between a given line and an asymptotic parallel line at a a. This angle is given by \Pi(a) = \frac{\pi}{2} - \gd(a), where \sin(\Pi(a)) = \sech a. This formula underscores the function's utility in analyzing asymptotic behaviors and in models like the Poincaré half-plane, where vertical lines and semicircles serve as geodesics, and the Gudermannian aids in coordinate transformations for computing distances and areas. The inverse Gudermannian, \gd^{-1}(\phi) = \int_0^\phi \sec u \, du = \ln \left| \tan\left(\frac{\pi}{4} + \frac{\phi}{2}\right) \right|, extends these connections to metric computations in hyperbolic spaces. In the , it parametrizes hyperbolic distances along , such as the distance from the apex of a semicircular geodesic to a point, expressed via \gd^{-1} of the . More broadly, it derives area formulas for horocyclic figures, like the area of an elliptic horocyclic segment S = 2\rho^2 \gd^{-1}(a/(2\rho)) - \sin(a/(2\rho)), and volumes in three-dimensional , including the finite volume V = \pi \rho^3 \gd^{-1}(r/\rho) - \sin(r/\rho), where \rho is the curvature radius. These applications facilitate precise calculations of plane figures and spatial volumes in Lobachevskian geometry. Structurally, the Gudermannian relates a circular \phi to a \psi such that \cos \phi = 1 / \cosh \psi, providing a reciprocal cosine linkage that supports derivations in hyperbolic and extensions to multidimensional settings. This reciprocity is leveraged in inequalities for hyperbolic , where expressions like (\gd a)^2 + (\gd b)^2 \leq 2(\cosh^2(a/2) + \sinh^2(b/2)) bound side relations using series expansions of \gd. Such tools enhance conceptual understanding of triangle properties and defect in hyperbolic planes.

Other mathematical and physical contexts

The Gudermannian function appears in the study of Euler's elasticae, which describe the equilibrium shapes of thin rods under forces, a fundamental problem in . For borderline elasticae with monotone curvature that approach a straight line at , the tangential \theta_b(s) as a function of s is given by \theta_b(s) = 2 \gd(s), where \gd denotes the Gudermannian ; equivalent expressions include $2 \arcsin(\tanh s) or $4 \arctan(e^s) - \pi. In , the Gudermannian function emerges in soliton solutions of the sine-Gordon equation, a modeling phenomena such as dislocations in crystals, in superconductors, and self-induced transparency in . When the elliptic modulus k = 1, the Jacobi amplitude function reduces to the Gudermannian, yielding the single-soliton solution \phi_S(x, t) = \pm 4\beta \gd\left(\beta \sqrt{2\alpha} \frac{x - x_0 - vt}{\sqrt{1 - v^2/c^2}}\right), where \beta and \alpha are parameters, v is the (v < c), and the total energy is finite at E(\phi_S) = 8\sqrt{2} \alpha \beta. This form connects to Lorentz-invariant wave propagation in (1+1)-dimensional field theories. The function also features in quantum field theory applications to the dynamical Casimir effect, where accelerating mirrors in vacuum produce real photons from virtual pairs. In a specific model of time-dependent particle creation from such mirrors, hyperbolic functions analogous to the Gudermannian provide an exact trajectory for the mirror position that ensures finite energy and controlled particle flux, facilitating analytical computation of the Bogoliubov coefficients and radiation spectrum.

References

  1. [1]
    Gudermannian -- from Wolfram MathWorld
    The Gudermannian function is the odd function denoted either gamma(x) or gd(x) which arises in the inverse equations for the Mercator projection.
  2. [2]
  3. [3]
    Christoph Gudermann (1798 - 1852) - Biography - MacTutor
    Gudermann, at this time, was particularly interested in the theory of elliptic functions and in the expansion of functions by power series. In particular his ...Missing: Gudermannian | Show results with:Gudermannian
  4. [4]
    [PDF] Gudermann function
    A definition for the Gudermann function is: gdx = / x. 0 secht dt. The function is also called the Gudermannian function or the hyperbolic amplitude function.Missing: mathematical sources
  5. [5]
    A Structural Approach to Gudermannian Functions - ResearchGate
    Nov 9, 2023 · The Gudermannian function relates the circular angle to the hyperbolic one when their cosines are reciprocal. Whereas both such angles are ...
  6. [6]
    inverse Gudermannian function - PlanetMath
    Mar 22, 2013 · gd−1(x)=arsinh(tanx)=12ln1+sinx1−sinx=∫x0dtcost. gd - 1 ⁢ ( x ) = arsinh ⁢ ( tan ⁡ x ) = 1 2 ⁢ ln ⁡ 1 + sin ⁡ x 1 - sin ⁡ x = ∫ 0 x d ⁢ t ...Missing: mathematics | Show results with:mathematics
  7. [7]
    Inverse Gudermannian -- from Wolfram MathWorld
    The inverse function of the Gudermannian y=gd^(-1)phi gives the vertical position y in the Mercator projection in terms of the latitude phi and may be defined ...
  8. [8]
    The Gudermannian Growth Model: Theory, Application and ... - SciELO
    This work introduces the Gudermannian function as an option for modeling sigmoid growth curves. The original function was transformed.
  9. [9]
    A028296 - OEIS
    Gudermannian(x) = x - (1/6)*x^3 + (1/24)*x^5 - (61/5040)*x^7 + (277/72576)*x^9 + .... Absolute values are the Euler numbers A000364. Cf. A019669.
  10. [10]
    Hyperbolic secant: Series representations - Wolfram Functions
    Elementary Functions > Sech[z]. Series representations. Generalized power series. Expansions at z==z0. For the function itself. > Expansions at z==0.
  11. [11]
    [PDF] arXiv:2311.10526v3 [astro-ph.CO] 21 Mar 2024
    Mar 21, 2024 · are related through the Gudermannian function gd(x) = arctan(sinh(x)) = 2 arctan(tanh(x/2)) (3.1) so that for this blueshifting component we ...
  12. [12]
    Abramowitz and Stegun (AMS55), p. 77 - ConvertIt.com
    Online Edition of AMS55: Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables. ... Gudermannian Function ?r gd z=2 arctan ez-- 2 4.3.
  13. [13]
    Gudermannian Function—Wolfram Documentation
    Mathematical function, suitable for both symbolic and numerical manipulation. The Gudermannian function is generically defined by . Gudermannian[z] has ...Missing: sources | Show results with:sources
  14. [14]
    gud.py
    ... Gudermannian function. # # Discussion: # # The Gudermannian function relates the hyperbolic and trigonometric # functions. For any argument X, there is a ...
  15. [15]
    4.34 Derivatives and Differential Equations
    §4.34 Derivatives and Differential Equations ; 4.34.5, d d z ⁡ sech ⁡ z ; ⓘ. Symbols: d f d x : derivative of f with respect to x , sech ⁡ z : hyperbolic secant ...
  16. [16]
    [PDF] FORMULAS FOR nth ORDER DERIVATIVES OF HYPERBOLIC ...
    Jul 19, 1971 · Formulas for the higher derivatives of tanh, sech, and csch, which may be derived in a similar way, are tabulated in the next section. A method ...
  17. [17]
    The Complex Gudermannian Function
    ### Integral Representations of the Gudermannian Function
  18. [18]
    A Structural Approach to Gudermannian Functions
    Nov 9, 2023 · The Gudermannian function relates the circular angle to the hyperbolic one when their cosines are reciprocal.
  19. [19]
    [PDF] History of the Integral of the Secant
    To see how this happened, we must first discuss sailing and early maps so that we can explain why Mercator invented his famous map projection. From the time of ...Missing: Gudermannian | Show results with:Gudermannian
  20. [20]
    A Theory and Generalization of the Circular and Hyperbolic Functions
    " The name was given in honor of Gudermann who called the function the "longitude of u" (Crelle's Journal, Vol. 6, 1830, p. 165).-EDITOR. 2 The Smithsonian ...
  21. [21]
    Johann Heinrich Lambert (1728 - 1777) - Biography - MacTutor
    Aug 29, 2025 · ... transcendental and Lindemann proved that π is transcendental. Lambert also made the first systematic development of hyperbolic functions. A ...
  22. [22]
    (PDF) On Elliptic and Hyperbolic Modular Functions and the ...
    Oct 16, 2025 · In this article, we move back almost 200 years to Christoph Gudermann, the great expert on elliptic functions, who successfully put the ...
  23. [23]
    Earliest Known Uses of Some of the Words of Mathematics (G)
    May 28, 2018 · The object of this paper is to discuss a general type of correspondence between structures which I have called Galois connexions.<|control11|><|separator|>
  24. [24]
    [PDF] Slide Rules with Hyperbolic Functions
    Then Cayley, in 1862, in recognition of Gudermann's contributions, proposed the name guder- mannian for the angle that Lambert called transcendent. The name, ...<|control11|><|separator|>
  25. [25]
    [PDF] A description of several coordinate systems for hyperbolic spaces
    Dec 19, 2001 · For flat, Euclidean geometry, it is easy to find in textbooks all the sys- tems of coordinates where the Helmholtz equation is solvablea: 4 ...<|control11|><|separator|>
  26. [26]
    None
    ### Sections Discussing Gudermannian Function and Its Extension to Higher Dimensions
  27. [27]
  28. [28]
    DLMF: §22.16 Related Functions ‣ Properties ‣ Chapter 22 Jacobian Elliptic Functions
    ### Summary of Gudermannian Function and Related Connections in §22.16
  29. [29]
    Mercator Projection -- from Wolfram MathWorld
    ### Role of the Gudermannian Function in the Mercator Projection
  30. [30]
    [PDF] Some Principles of Web Mercator Maps and their Computation
    Mar 5, 2020 · As a short excursion, we consider the Gudermann or Gudermannian function gd (Weisstein 2020) also known as Mercator function (Walz 2017) ...Missing: cartography | Show results with:cartography
  31. [31]
  32. [32]
    (PDF) Transverse Mercator with an accuracy of a few nanometers
    Implementations of two algorithms for the transverse Mercator projection are described; these achieve accuracies close to machine precision.<|control11|><|separator|>
  33. [33]
    [PDF] Hyperbolic geometry - Brice Loustau
    Oct 15, 2021 · Hyperbolic geometry is a very special subject: it is the star of geometries, and geometry is the star of mathematics!<|control11|><|separator|>
  34. [34]
    (PDF) The inverse Gudermannian in the hyperbolic geometry
    Oct 20, 2017 · ... inverse Gudermannian is the function of. an angle of quasiparallelism. 4. The identity law in a rectangular trihedral of the ehp(I) type on c.<|control11|><|separator|>
  35. [35]
    None
    ### Summary of Gudermannian Function in Hyperbolic Triangle Inequalities
  36. [36]