Worked Examples in the Dual-Number Algebra
Introduction
This article collects explicit computations in the dual-number algebra $\mathbb{D}'_R = R[\varepsilon]/(\varepsilon^2)$, in the spirit of the worked computation that accompanies the algebra of a number system. It follows Dual-Numbers Algebra for the conventions, Dual-Numbers Norm and Invertibility for the criterion for invertibility, Dual-Numbers Ideals and the Maximal Ideal for the ideal structure, Dual-Numbers Zero Divisors for the zero-divisor classification. Every computation was recomputed in exact rational arithmetic; where a floating-point value appears it is flagged.
The treatment is mathematically honest: every claim is either proved or stated as a definition. No physics is invoked. Unless a statement is labelled otherwise, the base is $R = \mathbb{R}$ and the algebra is $\mathbb{D}' = \mathbb{D}'_{\mathbb{R}}$. A general dual number is
$$ A = a + \varepsilon a', \qquad a, a' \in \mathbb{R}, $$
with $a = \operatorname{Re} A$, $a' = \operatorname{Inf} A$, dual conjugation $\bar A = a - \varepsilon a'$, unit criterion $a \neq 0$, and maximal ideal $\mathrm{M} = (\varepsilon)$. The two distinguished submodules are $R_{\mathbb{D}'} = \mathbb{R}\cdot 1$ and $\varepsilon R_{\mathbb{D}'} = \varepsilon\mathbb{R}$.
Each computation below is a worked instance of a general result proved in a companion article; the instance is stated with its reference, so that the article remains a table of examples rather than a second proof of the theory.
The Algebra in Coordinates
The Basis and the Product
The algebra has the basis $\{1, \varepsilon\}$ with the multiplication table
| $\cdot$ | $1$ | $\varepsilon$ |
|---|---|---|
| $1$ | $1$ | $\varepsilon$ |
| $\varepsilon$ | $\varepsilon$ | $0$ |
Example (product, worked). For $A = 2 + 3\varepsilon$ and $B = 4 + 5\varepsilon$,
$$ AB = (2\cdot 4) + (2\cdot 5 + 3\cdot 4)\varepsilon = 8 + 22\varepsilon. $$
The real parts multiply, and the infinitesimal part is $a\,b' + a'\,b$; the cross terms do not interact, because $\varepsilon^2 = 0$.
Example (squaring, worked). For $A = a + \varepsilon a'$, the square is
$$ A^2 = a^2 + 2aa'\,\varepsilon. $$
In particular $(2 + 3\varepsilon)^2 = 4 + 12\varepsilon$ and $(\varepsilon)^2 = 0$.
The Conjugate
Example (worked). $\overline{2 + 3\varepsilon} = 2 - 3\varepsilon$, and
$$ (2 + 3\varepsilon)(2 - 3\varepsilon) = 4 + (-6 + 6)\varepsilon = 4. $$
The Two Involutions
The Identity Involution
The identity map $\operatorname{id}(A) = A$ is an involutive algebra automorphism. Its fixed-point submodule is the whole algebra $\mathbb{D}'$, and the associated $\pm1$-eigenspace decomposition is the trivial one, with eigenvalue $+1$ everywhere.
The Dual Conjugation
Dual conjugation is
$$ \bar{\cdot} : a + \varepsilon a' \longmapsto a - \varepsilon a'. $$
Example (worked). $\overline{2 + 3\varepsilon} = 2 - 3\varepsilon$; applying it twice returns $2 + 3\varepsilon$, so $\bar{\bar A} = A$. It is an algebra automorphism: $\overline{AB} = \bar A\bar B$, checked on the worked product by $\overline{8 + 22\varepsilon} = 8 - 22\varepsilon$ and $(2 - 3\varepsilon)(4 - 5\varepsilon) = 8 - 22\varepsilon$.
The Fixed-Point Submodules
Proposition. The $+1$-eigenspace of dual conjugation is the real submodule
$$ R_{\mathbb{D}'} = \{A : \bar A = A\} = \{a : a \in \mathbb{R}\}, $$
and the $-1$-eigenspace is the infinitesimal submodule
$$ \varepsilon R_{\mathbb{D}'} = \{A : \bar A = -A\} = \{\varepsilon a' : a' \in \mathbb{R}\}. $$
Proof. $\bar{a + \varepsilon a'} = a - \varepsilon a'$ equals $a + \varepsilon a'$ exactly when $a' = 0$, and equals $-(a + \varepsilon a')$ exactly when $a = 0$; the two sets are as displayed.
Example (worked projections). For $A = 2 + 3\varepsilon$,
$$ \tfrac{1}{2}(A + \bar A) = 2, \qquad \tfrac{1}{2}(A - \bar A) = 3\varepsilon, $$
and indeed $2 + 3\varepsilon = A$.
The Intersection
Proposition. $R_{\mathbb{D}'} \cap \varepsilon R_{\mathbb{D}'} = \{0\}$, and $\mathbb{D}' = R_{\mathbb{D}'} \oplus \varepsilon R_{\mathbb{D}'}$.
Proof. An element of the intersection satisfies $a' = 0$ and $a = 0$, so it is zero. The sum decomposition is the worked projection identity above.
So the algebra has exactly one nontrivial involution, and the two submodules are the eigenspaces. The biquaternion algebra, by contrast, carries four conjugations and six eigenspace-type subspaces; here the lattice degenerates to two submodules and the single relation $R_{\mathbb{D}'} \cap \varepsilon R_{\mathbb{D}'} = 0$.
The Nilpotents and the Maximal Ideal
Squaring the Maximal Ideal
Example (worked). Every element of $\mathrm{M} = \varepsilon\mathbb{R}$ is nilpotent of index two:
$$ (3\varepsilon)^2 = 9\varepsilon^2 = 0, \qquad (-5\varepsilon)^2 = 25\varepsilon^2 = 0, \qquad (7\varepsilon)^2 = 0. $$
So $\mathrm{M}^2 = 0$ and $\mathrm{M} \neq 0$: the maximal ideal is nilpotent of index exactly two.
Example (worked). A mixed element is never nilpotent: $(2 + 3\varepsilon)^2 = 4 + 12\varepsilon \neq 0$, and more generally $(a + \varepsilon a')^n = a^n + na^{n-1}\varepsilon a'$, which vanishes only when $a = 0$.
The Powers of an Element
Proposition. For every integer $n \geq 1$, $(a + \varepsilon a')^n = a^n + n\,a^{n-1}\varepsilon a'$, interpreted for $n = 1$ as the element itself.
Proof. Induction on $n$ using $(a + \varepsilon a')^n(a + \varepsilon a') = a^{n+1} + (n a^{n-1}a'\,a + a^n a')\varepsilon = a^{n+1} + (n+1)a^n \varepsilon a'$.
Example (worked). $(1 + \varepsilon)^{5} = 1 + 5\varepsilon$, and $(2 + \varepsilon)^3 = 8 + 12\varepsilon$.
Explicit Zero-Divisor Pairs
The General Pair
By Dual-Numbers Zero Divisors the zero divisors are exactly the nonzero elements of $\mathrm{M}$, that is, the elements $\varepsilon a'$ with $a' \neq 0$, and any two of them multiply to zero:
$$ (\varepsilon a')(\varepsilon b') = a' b'\,\varepsilon^2 = 0 \qquad \text{for all } a', b'. $$
A Table
Each row is a pair of nonzero zero divisors whose product is zero; the product is computed in the last column.
| $A$ | $B$ | $AB$ |
|---|---|---|
| $\varepsilon$ | $\varepsilon$ | $0$ |
| $2\varepsilon$ | $5\varepsilon$ | $0$ |
| $-3\varepsilon$ | $\varepsilon$ | $0$ |
| $7\varepsilon$ | $-4\varepsilon$ | $0$ |
The point is that the annihilation is not special to any pair: it holds for every pair of nonzero elements of the line $\mathrm{M}$. Equivalently, $\operatorname{Ann}(\varepsilon a') = \mathrm{M}$ for every $a' \neq 0$, so every nonzero zero divisor annihilates the entire punctured line.
Example (worked, mixed pair). The element $2 + 3\varepsilon$ is not a zero divisor and annihilates nothing nonzero:
$$ (2 + 3\varepsilon)(b + \varepsilon b') = 2 b + (2 b' + 3 b)\varepsilon, $$
which vanishes only when $b = 0$ and $2b' + 3b = 0$, that is $b = b' = 0$.
Units and Inverses, Worked
Units
Example (worked). $A = 2 + 3\varepsilon$ has real part $2 \neq 0$, so it is a unit. Its inverse is
$$ A^{-1} = \frac{1}{2} - \frac{3}{4}\,\varepsilon, $$
and the check is
$$ (2 + 3\varepsilon)\Bigl(\tfrac{1}{2} - \tfrac{3}{4}\varepsilon\Bigr) = 1 + \Bigl(-\tfrac{3}{2} + \tfrac{3}{2}\Bigr)\varepsilon = 1. $$
Non-Units
Example (worked). $A = 0 + 3\varepsilon$ has real part $0$, so it is not a unit. Any attempted inverse fails: if $(3\varepsilon)(b + \varepsilon b') = 3b\,\varepsilon = 1$, then $3b = 1$ and $0 = 1$, a contradiction.
Example (worked). $A = -4 + \varepsilon$ has real part $-4 \neq 0$, so it is a unit, since the criterion is $a \neq 0$ and not $a > 0$:
$$ (-4 + \varepsilon)^{-1} = -\frac{1}{4} - \frac{1}{16}\varepsilon. $$
The Inverse Worked, in General
Example (worked). For $A = a + \varepsilon a'$ with $a \neq 0$, the inverse formula gives
$$ A^{-1} = \frac{1}{a} - \frac{a'}{a^2}\,\varepsilon, $$
and the general check is
$$ (a + \varepsilon a')\Bigl(\frac{1}{a} - \frac{a'}{a^2}\varepsilon\Bigr) = 1 + \Bigl(-\frac{a'}{a} + \frac{a'}{a}\Bigr)\varepsilon = 1. $$
The Matrix Model and the Shear
The Matrix of a Worked Product
By Dual-Numbers Matrix Element Representation the algebra embeds faithfully in $M_2(\mathbb{R})$ through
$$ \Phi(a + \varepsilon a') = \begin{pmatrix} a & a' \\ 0 & a \end{pmatrix}. $$
Worked, $\Phi(2 + 3\varepsilon) = \begin{pmatrix} 2 & 3 \\ 0 & 2 \end{pmatrix}$ and $\Phi(4 + 5\varepsilon) = \begin{pmatrix} 4 & 5 \\ 0 & 4 \end{pmatrix}$, and
$$ \Phi(2 + 3\varepsilon)\,\Phi(4 + 5\varepsilon) = \begin{pmatrix} 8 & 2\cdot 5 + 3\cdot 4 \\ 0 & 8 \end{pmatrix} = \begin{pmatrix} 8 & 22 \\ 0 & 8 \end{pmatrix} = \Phi(8 + 22\varepsilon), $$
matching the coordinate product worked above. The matrix is upper triangular with equal diagonal entries, and $\Phi(1 + \varepsilon) = \begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix}$.
Multiplication by a Unit Is a Shear
Let $A = a + \varepsilon a'$ be a unit, so $a \neq 0$, and write $s = a'/a$, giving $A = a(1 + s\varepsilon)$. Acting on $B = b + \varepsilon b'$,
$$ A B = a(1 + s\varepsilon)(b + \varepsilon b') = a b + (a b' + a s b)\varepsilon = a\bigl(b + (b' + s b)\varepsilon\bigr). $$
In the coordinates $(b, b')$ of $B$ this is the linear map
$$ (b, b') \longmapsto \bigl(a b,\ a(b' + s b)\bigr), $$
a dilation by $a$ followed by the shear $(b, b') \mapsto (b, b' + s b)$; the shear alone is multiplication by $1 + s\varepsilon$, which fixes the real axis pointwise.
Example (worked). For $A = 1 + \varepsilon$ and $B = 2 + 5\varepsilon$,
$$ (1 + \varepsilon)(2 + 5\varepsilon) = 2 + (5 + 2)\varepsilon = 2 + 7\varepsilon, $$
so the point $(2, 5)$ goes to $(2, 7)$, raised by the shear amount $s b = 1 \cdot 2 = 2$.
Example (worked). For $A = 2 + 3\varepsilon$ one has $a = 2$, $s = 3/2$, and for $B = 2 + 5\varepsilon$,
$$ (2 + 3\varepsilon)(2 + 5\varepsilon) = 4 + (10 + 6)\varepsilon = 4 + 16\varepsilon = 2\Bigl(2 + \bigl(5 + \tfrac{3}{2}\cdot 2\bigr)\varepsilon\Bigr), $$
a dilation by $2$ of the sheared point $(2, 8)$.
Worked Hyperbolic Elements
The Two Factors of a Unit
Every unit factors as $A = a(1 + s\varepsilon)$ with $a \neq 0$ and $s = a'/a$: the real factor $a$ is the scale and $1 + s\varepsilon$ the parabolic factor.
Example (worked). $2 + 3\varepsilon = 2\bigl(1 + \tfrac{3}{2}\varepsilon\bigr)$ and $-4 + \varepsilon = -4\bigl(1 - \tfrac{1}{4}\varepsilon\bigr)$.
Powers of a Worked Element
Example (worked). For $A = 2 + 3\varepsilon$,
$$ (2 + 3\varepsilon)^3 = 2^3 + 3\cdot 2^2 \cdot 3\,\varepsilon = 8 + 36\varepsilon, $$
and in the factorised form $2^3\bigl(1 + \tfrac{3}{2}\varepsilon\bigr)^3 = 8\bigl(1 + 3\cdot\tfrac{3}{2}\varepsilon\bigr) = 8 + 36\varepsilon$, the two computations agreeing.
The Hyperbolic Case Degenerates
In the split-complex algebra the units $a + jb$ satisfy $a^2 - b^2 = 1$ and form a hyperbola, a one-parameter multiplicative family with a genuine hyperbolic parameter; in the complex algebra the units form a circle. In $\mathbb{D}'$ the corresponding family of units is the pair of parallel lines $a = \pm 1$, and each line is carried by the parabolic family $1 + s\varepsilon$ with the additive parameter $s$. The hyperbola and the circle have both degenerated to lines, the hyperbolic parameter being replaced by the shear parameter; this is the arithmetic content of $\varepsilon^2 = 0$.
Relation Between the Decompositions
The two decompositions of the dual-number algebra coincide, and the worked examples make the coincidence explicit.
- The conjugate decomposition is $\mathbb{D}' = R_{\mathbb{D}'} \oplus \varepsilon R_{\mathbb{D}'}$, from the eigenspaces of $\bar{\cdot}$; worked for $A = 2 + 3\varepsilon$ as $A = 2 + 3\varepsilon$.
- The real–infinitesimal decomposition is the same identity written as $A = \operatorname{Re}(A) + \operatorname{Inf}(A)\varepsilon$.
There is no third decomposition, and no nontrivial Peirce decomposition: the only idempotents are $0$ and $1$, so the corners of Dual-Numbers Ideals and the Maximal Ideal are trivial. The single relation $R_{\mathbb{D}'} \cap \varepsilon R_{\mathbb{D}'} = \{0\}$ is the only intersection of the two submodules, and the worked example $2 + 3\varepsilon = A$ is its arithmetic content.
Summary
The worked computations in the dual-number algebra are these. The two involutions are the identity and dual conjugation; the identity has fixed set all of $\mathbb{D}'$, while dual conjugation has $+1$-eigenspace $R_{\mathbb{D}'} = \mathbb{R}$ and $-1$-eigenspace $\varepsilon R_{\mathbb{D}'} = \varepsilon\mathbb{R}$, meeting only in $0$. Every element of the maximal ideal is nilpotent of index two, $(\varepsilon a')^2 = 0$, while a mixed element has $(a + \varepsilon a')^n = a^n + n a^{n-1}\varepsilon a'$ and is never nilpotent unless $a = 0$. The zero divisors are the nonzero elements $\varepsilon a'$ of $\mathrm{M}$, and any two of them multiply to zero; for instance $(2\varepsilon)(5\varepsilon) = 0$. The criterion for a unit is $a \neq 0$, and the inverse is $a^{-1} - a^{-2}\varepsilon a'$; the check $(2 + 3\varepsilon)(\tfrac{1}{2} - \tfrac{3}{4}\varepsilon) = 1$ is the worked instance. Under the matrix model $\Phi(a + \varepsilon a') = \begin{pmatrix} a & a' \\ 0 & a \end{pmatrix}$ the worked product $\Phi(2 + 3\varepsilon)\Phi(4 + 5\varepsilon) = \Phi(8 + 22\varepsilon)$ reproduces the coordinate product, and multiplication by the unit $a(1 + s\varepsilon)$ acts as the dilation by $a$ followed by the shear $(b, b') \mapsto (b, b' + s b)$. The two decompositions of the algebra coincide, and no nontrivial Peirce decomposition exists.
Summary of Notation
| Symbol | Meaning |
|---|---|
| $\mathbb{D}' = \mathbb{R}[\varepsilon]/(\varepsilon^2)$ | Dual-number algebra; $\varepsilon^2 = 0$ |
| $A = a + \varepsilon a'$ | General dual number |
| $a = \operatorname{Re} A$, $a' = \operatorname{Inf} A$ | Real and infinitesimal parts |
| $\bar{A} = a - \varepsilon a'$ | Dual conjugation, the nontrivial involution |
| $R_{\mathbb{D}'}$, $\varepsilon R_{\mathbb{D}'}$ | Real and infinitesimal submodules |
| $\mathrm{M} = (\varepsilon)$ | Maximal ideal, $\mathrm{M}^2 = 0$ |
| $\Phi(a + \varepsilon a') = \begin{pmatrix} a & a' \\ 0 & a \end{pmatrix}$ | The faithful matrix model |
| $a(1 + s\varepsilon)$ | Factorisation of a unit, $s = a'/a$ |
| $(b, b') \mapsto (b, b' + s b)$ | The shear of multiplication by $1 + s\varepsilon$ |
Further Reading
- Eduard Study, Geometrie der Dynamen (Teubner, Leipzig, 1903), for the classical coordinate computations in the dual numbers and their geometry.
- I. M. Yaglom, Complex Numbers in Geometry (Academic Press, New York, 1968), for worked computations in the dual, complex and split-complex algebras side by side.
- John H. Conway and Derek A. Smith, On Quaternions and Octonions (A K Peters, Natick, 2003), for worked examples in the two-dimensional real algebras.
- Andreas Griewank and Andrea Walther, Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation (SIAM, Philadelphia, 2008), for the computed algebra of nilpotent extensions.