Preprint
Article

This version is not peer-reviewed.

The Quadratic Equation for the Quaternions: The Closed Form Solution

Submitted:

29 April 2023

Posted:

30 April 2023

Read the latest preprint version here

Abstract
This article derives the closed form solution to finding the roots of a quadratic equation in SU3.
Keywords: 
;  ;  
Theorem 1. 
The Quadratic Equation , The General Depressed Case
Quadratic Quaternionic Calculator, Closed Form Solution
Let c = x + a   x + b   = x 2 + x   b + a   x + a   b , which can also be written: c = x 2 +   b   4 R +   a   4 L x + a   b We know that there exists a two dimensional basis in the four dimensional space of quaternions that describes the given roots. Namely the bisector of the roots (mistakenly known as the Axis of Symmetry), and the straight line that is between the two roots and the bisector.
We shall define the bisector as u = 1 2 a + b   as the first vector of the basis.
We now define the locator as v = 1 2 a b  
Preprints 72263 i001
EQ.1 u = 1 2 a + b EQ.2 v = 1 2 a b   , compelling a = u + v and b = u v EQ.3 Let t = x + u , therefore x = t u EQ.4a c = x 2 + x   b + a   x + a b EQ.5 c = x 2 + x u v + u + v x + u + v u v EQ.6 c = t u 2 + t u u v + u + v t u + u + v u v EQ.7 c = t 2 t u u t + u 2 + t u t v u 2 + u v + u t u 2 + v t v u + u 2 u v + v u v 2 EQ.8a c = t 2 t v + v t v 2 = t + v t v EQ.8b If c is upon the same great circle as v , then the traditional quadratic equation suffices to resolve the roots:
Let d = c + v 2   ;         0 = t 2 t v + v t + d , then l i m λ c λ v t = ± d = ± v 2 + c , where λ is the Great Circle Function (see page 5).
However, in the general case, it would be exceptionally rare for λ c = λ v , so we proceed.
EQ.9a c = t 2 t v + v t v 2 EQ.9b d = t 2 t v + v t                                           ;           d = c + v 2 .
EQ.9c d = t 2   v   4 L +   v   4 R t       . This is the fundamental middle handed form, from which we derive the solution.
EQ.10a Let w = t + v , such that t = w v EQ.11a c =                       w v 2                     w v v + v w v v 2 EQ.12a c = w 2 w v v w + v 2 w v v 2 + v w v 2 v 2 EQ.13a c = w 2 w v v w + v 2 w v + v 2 + v w v 2 v 2 EQ.14a c = w 2 2 w v EQ.15a c = w 2 2 w v = w w 2 v EQ.16a Let y = 2 v EQ.17a c = w 2 + w y     = w   w + y   . This is the second fundamental form, the left-handed form.
EQ.9b c = t 2 t v + v t v 2 (a restatement of EQ 9a)
EQ.10b Let z = t v , such that t = z + v , thence:
w = z + 2 v = z y z = w 2 v = w + y EQ.11b c =                   z + v 2                           z + v v     +   v z + v   v 2 EQ.12b c = z 2 + z v + v z + v 2 + z v v 2 + + v z + v 2 v 2 EQ.13b c = z 2 + 2 v z = z 2 y z =   z y   z . This is the third fundamental form, the right-handed form.
However, before we can proceed to resolve the roots of d = t 2 t v + v t     , some general definitions and lemmas are in order.
Definition 1. 
Orthogonal Imaginary Bases
EQ.7d i   j = k and λ μ = ν ; furthermore that, j   i = + k and μ λ = + ν .
i   k = + j and λ ν = + μ ; furthermore that, k   i = j and ν λ = μ .
j   k = i and μ ν = λ ; furthermore that, k   j = + i and ν μ = + λ .
If, and only if:
1.
Let α = A T A N 2 λ 2 λ 1 ; let β = A T A N 2 λ 3 λ 1 2 + λ 2 2 ; λ = λ 1 i + λ 2 j + λ 3 k = a i + b j + c k a 2 + b 2 + c 2
2.
λ = + i   c o s α c o s β + j s i n α c o s β + k s i n β
3.
μ =               i   s i n α + j c o s α                   + 0 k
4.
ν = i   c o s α s i n β + j s i n α s i n β + k c o s β
5.
This rigid body rotation of …
a.
i = 1 i + 0 j + 0 k to λ
b.
j = 0 i + 1 j + 0 k to μ
c.
k = 0 i + 0 j + 1 k to ν
which maintains the relative spatial property that ν   μ = λ = μ   ν     ;     μ   4 L ν = + μ   4 R ν .
Definition 2. 
A Two-Dimensional Plane in Hypercomplex Space, Absolute Shape.
In Cartesian space, we can define a plane as 0 = a x + b y + c z , however, we cannot do so over the quaternions, there are no real numbers for 72 a , b , c , d such that 0 = d q + a i + b j + c k . Rather, we define a plane as a continuum of points that satisfy the below property.
Space and Time Newton and Leibniz https://youtu.be/ZuDkAQvbWj8
EQ.8d Let P be the set of all points in the form p = q c o s θ + λ s i n θ + Ψ , where λ is a purely imaginary unit vector, and Ψ is a sum of
all remaining pairwise orthogonal imaginary unit vectors that are orthogonal to λ , such that there exists an chiral isomorphism from
λ , μ 1 ,   μ 2 . . . u n 2   onto d 1 ,   d 2 , d 3 . . . d n where n is the number of hypercomplex dimensions and that all d z is the set of orthogonal imaginary unit
vectors that this orthogonal basis is Cayley Geodesic, Orientable, Chiral and furthermore Cayley Algebraic (see Appendix A).
Thence Ψ = z = 1 z = n 2 α z μ z = z = 2 z = n β z d z ; such that, z , α z , β z R .
For the quaternions Ψ = α 1 μ + α 2 ν .
Regardless of the position and facing of different observers, all observers will agree that a collection of objects in P (from our perspective) defines a plane from all perspectives. Thus, position and orientation may be relative, but ShapeShape is Absolute, and Quadratic Equations have a well-defined Shape. In the below image there are three parabolas…but they are all the same, just from three different perspectives.
Preprints 72263 i002
Definition 3. 
Affine Two Dimensional Planes Hypercomplex Space.
EQ.9d Let P 1 be the hypercomplex plane such that p 1 P 1 , p 1 = q c o s θ + λ 1 s i n θ + Ψ 1   ,   θ   ;   Ψ 1 = α 1 μ 1 + β 1 ν 1   ,     α 1 , β 1 R EQ.10d Let P 2 be the hypercomplex plane such that p 2 P 2 , p 2 = q c o s ϕ + λ 2 s i n ϕ + Ψ 2   ,   ϕ ;   Ψ 2 = α 2 μ 2 + β 2 ν 2   ,     α 2 , β 2 R
Corollary 4. 
Parallel Two Dimensional Planes Hypercomplex Space.
If λ 1 = λ 2 and Ψ 1 Ψ 2 , then P 1 and P 2 define two parallel two-dimensional subspaces, separated by both  Ψ 3 , Ψ 4 = ± 1 2 1 4 Ψ 1 Ψ 2 2 , as preference can be given to neither Ψ 1 Ψ 2 nor Ψ 2 Ψ 1 and that Ψ 1 Ψ 2 Ψ 2 Ψ 1 , such that 1 4 Ψ 1 Ψ 2 2 define the unique line of intersection between two congruent cylinders about λ 1 and λ 2 with radii equal to one-half the real number magnitude of Ψ 1 Ψ 2 in λ , μ , ν space (which is rotated i ,   j ,   k space, they are the same space, with different bases).
Preprints 72263 i003
Lemma 5. 
Cayley Multiplication over the General Quaternions, The Square of a Quaternion, Polar and Cartesian Forms.
The rules of x 2 are well defined over any imaginary orthogonal basis. Let first separate x as x = x 1 + x 2 .
EQ.1e x = x 0 q + x 1 λ + x 2 μ + x 3 v = x 1 + x 2 EQ.2e x 1 = M 1 q c o s θ + λ s i n θ ; M 1 = x 0 2 + x 1 2 ; θ = A T A N 2 x 1 x 0 ; x 1 2 = M 1 2 q c o s 2 θ + λ s i n 2 θ EQ.3e x 2 = M 2 μ c o s ϕ + ν s i n ϕ ; M 2 = x 2 2 + x 3 2 ; ϕ = A T A N 2 x 3 x 2 ; x 2 2 = M 2 2 q EQ.4e x 1 x 2 = M 1 M 2 μ c o s θ c o s ϕ + ν c o s θ s i n ϕ v s i n θ c o s ϕ + u s i n θ s i n ϕ EQ.5e x 2 x 1 = M 1 M 2 μ c o s ϕ c o s θ + ν c o s ϕ s i n θ + v s i n ϕ c o s θ u s i n ϕ s i n θ EQ.6e x 1 x 2 + x 2 x 1 = 2 M 1 M 2 μ c o s ϕ + ν s i n ϕ c o s θ EQ.7e x 2 = M 2 2 q + M 1 2 q c o s 2 θ + λ s i n 2 θ + 2 M 1 M 2 μ c o s ϕ + ν s i n ϕ c o s θ EQ.8e x 2 = x 0 2 x 1 2 x 1 2 x 1 2 q + 2 x 0 x 1 λ + x 2 μ + x 3 ν (this is well known Cartesian form of a quaternionic square)
EQ.9e                     x 0 2 x 1 2 x 1 2 x 1 2 q = M 2 2 q + M 1 2 q c o s 2 θ EQ.10e             2 x 0 x 2 μ + x 3 ν = 2 M 1 M 2 μ c o s ϕ + ν s i n ϕ c o s θ EQ.11e                     2 x 0 x 1 λ = M 1 2 λ s i n 2 θ
Preprints 72263 i004
Definition 6. 
Hypercomplex Chiral Orthogonal Basis Conversion, For the Quaternions
EQ.1f x = c 0 q + c 1 i + c 2 j + c 3 k = d 0 q + d 1 λ + d 2 μ + d 3 ν ; c 0 = d 0 Since     c 0 = d 0 , we are only concerned with the imaginary bases, thus:
EQ.2f y = x c 0 q = c 1 i + c 2 j + c 3 k = d 1 λ + d 2 μ + d 3 ν We already know that:
1.
α = A T A N 2 c 2 c 1 ; β = A T A N 2 3 c 1 2 + c 2 2
2.
λ = + i   c o s α c o s β + j s i n α c o s β + k s i n β
3.
μ =                 i   s i n α + j c o s β                       + 0 k
4.
ν = i   c o s α s i n β + j s i n α s i n β + k c o s β
Let us rename the above as:
EQ.3f λ = γ 1,1 i + γ 1,2 j + γ 1,3 k EQ.4f μ = γ 2,1 i + γ 2,2 j + γ 2,3 k EQ.5f ν = γ 3,1 i + γ 3,2 j + γ 3,3 k Which yields the system of three linear equations:
EQ.5f d 1 γ 1,1 i + d 2 γ 2,1 i + d 3 γ 3,1 i = c 1 i EQ.6f d 1 γ 1,2 j + d 2 γ 2,2 j + d 3 γ 3,2 j = c 2 j EQ.7f d 1 γ 1,3 k + d 2 γ 2,3 k + d 3 γ 3,2 k = c 3 k EQ8f Let 𝚪 be a 3x3 real matrix whose pairwise entries are equal to γ m , n .
EQ9f Let C be a 1x3 real column matrix whose entries are c 1 ,   c 2 and c 3 respectively.
EQ10f Let D =𝚪-1 C , which is also a 1x3 real column matrix
Then d 1 , d 2 , d 3 are the respective entries of D from top to bottom.
Theorem 7. 
The Lambda Choice Quaternion Eraser
Statement One  0 = t v + v t when both t and v are on the same Great Circle of λ .
Statement Two  α 1 μ + α 2 v = t v + v t when both t and v are not on the same Great Circle and α 1 μ + α 2 v is orthogonal to both t and v .
Proof. 
EQ.1 Let v = v 0 q + v 1 λ , establishing v as the reference frame, compelling the orthogonal basis λ ,   μ ,   ν , then:
EQ.2 Let t = t 0 q + t 1 λ + t 2 μ + t 3 ν EQ.3 t v + v t = t 0 q + t 1 λ + t 2 μ + t 3 ν v 0 q + v 1 λ + v 0 q + v 1 λ t 0 q + t 1 λ + t 2 μ + t 3 ν EQ.4 0 = t 0 q t 1 λ v 0 q + v 1 λ + v 0 q + v 1 λ + t 0 q + t 1 λ , since they commute upon the same Great Circle.
EQ.5 t v + v t = t 2 μ t 3 ν v 0 q + v 1 λ + v 0 q + v 1 λ + t 2 μ + t 3 ν EQ.6 t v + v t = t 2 v 0 μ t 2 v 1 ν t 3 v 0 ν + t 3 v 1 μ + t 2 v 0 μ + t 3 v 0 ν t 2 v 1 ν + t 3 v 1 μ EQ.7 t v + v t = 2 t 3 v 1 μ t 2 v 1 ν , such that t v + v t is orthogonal to v , vanishing t 0 and t 1 and v 0 .
Q.E.D.
Likewise we could establish t as the reference frame via: t = t 2,0 q + t 2,1 λ 2 , compelling the orthogonal basis λ 2 , μ 2 , ν 2 , and the expression t v + v t will result in vector that is also orthogonal to t , vanishing the real part and λ 2 part of y , leaving only μ 2 and ν 2 as the remaining dimensions.
Regardless of which reference frame we choose, we know that t v + v t   is orthogonal to both t and v , and, by definition, anything in the form of M ( μ v c o s θ + ν v s i n θ ) is strictly orthogonal to v ; however, not everything in form of M ( μ v c o s θ + ν v s i n θ ) is strictly orthogonal to t .
Also observe that the matrix X = x   4 L + + x   4 R is non-invertible, since any vector with the same imaginary components, but any real component, will evaluate to X . This occurs not only for the hypercomplex numbers, but for the real and complex numbers as well. Consider the following equation:
EQ.8 0 = x y + y x . Obviously, any real value of x will suffice, and any complex value shall also suffice.
EQ.9 w = x y + y x . Even though w has a non-zero magnitude and is orthogonal to both x and y , any real or λ (in respect to y ) shall suffice to resolve x , hence the inherent invertibility of X = x   4 L + + x   4 R .
However, this does not imply that the appearance of some X = x   4 L + + x   4 R within an equation renders the solution impossible, rather it cannot appear alone, for instance z = x x y + y x most certainly has a solution for a given z and y , we would rewrite this as:
EQ.10 z = q x x y + y x = q 4 L   y   4 L +   y   4 R x   ;           x = q 4 L   y   4 L +   y   4 R 1 z .
Hence, the statement z = x 2 x y + y x also has a solution…two in fact.
However, the most important takeaway is that t v + v t = 2 t 3 v 1 μ t 2 v 1 ν , meaning the real part of v , which is v 0 , has no effect; thus, the equation d = t 2 t v + v t remains equal to d , no matter the real part of v . Allowing us to reduce the equation to (let M and N be positive reals):
EQ.11 d = t 2 t v + t v = t 2 t M q c o s θ + λ s i n θ + M q c o s θ + λ s i n θ t = t 2 t M λ s i n θ + M λ s i n θ t = t 2 N t λ + N λ t
Theorem 8. 
The Right Triangle Theorem
Thence, the expression: d = t 2 t v + v t geometrically compels d to be the hypotenuse of a right triangle, since t v + v t is orthogonal to both t and its square, as both t and t 2 lay upon the same Great Circle.
Although there exists an entire family of right triangles that share d as the hypotenuse, there are only two congruent right triangles within this family that satisfy t .
EQ.1 d = t 2 t v + v t   ;             v = v 0 q + v 1 λ v         λ v , μ v , ν v , which is the orthogonal basis in respect to v .
EQ.2 t v + v t = α μ v + β ν v EQ.3 μ v c o s θ + ν v s i n θ is orthogonal to μ v s i n θ + ν v c o s θ by definition. We shall choose this orthogonality to generate the family of solutions.
EQ.4 μ v c o s θ + ν v s i n θ is orthogonal to + μ v s i n θ ν v c o s θ by definition. We discard this orthogonality in favor of the former.
EQ.5 d = d 0 q + d 1 λ v + d 2 μ v + d 3 ν v EQ.6 t = t 0 q + t 1 λ v + t 2 μ v + t 3 ν v EQ.7 t 2 =       t 0 2 t 1 2 t 2 2 t 3 2 q + 2 t 0 t 1 λ v + 2 t 0 t 2 μ v + 2 t 0 t 3 ν v     ;                 d 0 =   t 0 2 t 1 2 t 2 2 t 3 2 ;       d 1 = 2 t 0 t 1 EQ.8 Let f = d 2 μ v + d 3 ν v   EQ.9 Let g = G   + μ v c o s β + ν v s i n β = t v + t v   = + 2 t 3 v 1 μ v   2 t 2 v 1 ν v = x μ v + y ν v EQ.10 Let t 2 = h = H μ v s i n β + ν v c o s β = + 2 t 0 t 2 μ v + 2 t 0 t 3 ν v = w μ v + z ν v EQ.11 f = G g + H h   EQ.12 d 2 μ v = + G μ v c o s θ H μ v s i n θ = + 2 t 3 v 1 μ v + 2 t 0 t 2 μ v = x μ v + w μ v EQ.13 d 3 ν v = + G ν v s i n θ + H ν v c o s θ = 2 t 2 v 1 ν v + 2 t 0 t 3 ν v = y ν v + z ν v EQ.14 t 2 = d g                 t = ± d g The above relationship provides us with enough information to brute the roots of the quadratic equation by simply comparing every value of the angular argument of β against the real number magnitude of error from the return on d in a preliminary search, and then converge rapidly upon the roots via bisection.
In fact, it was by empirical observation of the roots (using the above rapidly convergent algorithm) that I was able to resolve the closed form solution. We shall first simplify the quadratic equation further in lieu of those empirical results.
Theorem 9. 
The Orthogonal Basis Rotation Theorem; The Relative Frame Theorem.
EQ.1 Let d = d 0 q + d 1 λ v + d 2 μ v + d 3 ν v = A q c o s α + λ v s i n α + B μ v c o s ϕ + ν v s i n ϕ EQ.2a μ 2 = + μ v c o s ϕ + ν v s i n ϕ   EQ.2b ν 2 = μ v c o s ϕ + ν v c o s ϕ   EQ.3 d = d 0 q + d 1 λ v + d 2 μ v + d 3 ν v = A q c o s α + λ v s i n α + B μ 2 + 0 ν 2 We are able to perform this basis conversion because all we did was rotate μ v , ν v about λ v ; hence, λ v , μ 2 , ν 2 preserves the multiplicative relationships expected in the original basis. In fact, there is no preferred frame of reference for the μ and ν axes for an Observer on the Great Circle of λ , only λ is absolute from the Observer’s perspective. The Observer is free to rotate the μ v , ν v axes in any manner that simplifies the existing problem.
Thus, in the equation d = t 2 t v + v t , the v variable establishes λ , and the d variables establishes μ 2 and ν 2 . However, we have yet one more simplification to make:
Empirical observation revealed that the equation 1 B d = A B q c o s α + λ v s i n α + μ 2 = t 2 t v + v t , always yielded two roots, such that:
EQ.5 t = C ν 2 ± w 2 , where w = x q + y λ + z μ 2 . That is, when both roots were added together, they left no q ,   λ nor μ 2 part, and both roots always had the same ν 2 part. This was true for in excess of 100,000 random trials.
Although empirical observation cannot be used to claim a proof; most certainly assisted in realizing how to derive the closed form solution. Let us briefly examine why the form of t = C ν 2 ± w 2 would be the solution to our roots.
In the image below it shows ± w 2 as the purple vector, located at μ c o s θ for some arbitrary chosen theta. It follows that + w C ν 2 and w C ν 2 will form an isosceles triangle, since + w and w will have an equal but opposite reaction with C v . The green vectors are forced via t v + v t , since they must be orthogonal to both λ and t 2 . The question of course is which value of θ in the below diagram produces the correct pair of magnitudes for the blue and green legs that have a sum square of one when they interact with the magnitude of v .
Preprints 72263 i005
Theorem 10. 
The Fully Depressed Case of the Quadratic Equation
We now combine Theorems 14 and 15 to yield the fully depressed case of the quadratic equation. We shall remove the real part of v and rotate the orthogonal axes in respect to d and then divide by the orthogonal magnitude of d … and we shall simply write μ 2 as μ EQ.1a d = t 2 t v + v t = t 2 t v 0 q + N λ + v 0 q + N λ t , where λ v is in respect to v .
EQ.1b d = t 2 t N λ + N λ t EQ.1c d = ξ q c o s α + λ v s i n α + ω 1 μ v + ω 2 ν v , where μ v , ν v is the initial orthogonal basis in respect to v .
EQ.2a Let Ω = ω 1 2 + ω 2 2 EQ.2b Let ϕ = A T A N 2 ω 2 ω 1 EQ.2c Let μ 2 = + μ v c o s ϕ + ν v s i n ϕ EQ.2d Let ν 2 = μ v s i n ϕ + ν v c o s ϕ EQ.2e Let d = d 0 q + d 1 λ + Ω μ 2 + 0 ν 2 EQ.3a t = t 0 q + t 1 λ + t 2 μ 2 + t 3 ν 2 EQ.3b t 2 = t 0 2 t 1 2 t 2 2 t 3 2 q + 2 t 0 t 1 λ + 2 t 0 t 2 μ 2 + 2 t 0 t 3 ν 2 EQ.3c t N λ + N λ t = 2 N t 3 μ 2 2 N t 2 ν 2 (Lambda Choice Quaternion Eraser)
EQ.4a d 0 = t 0 2 t 1 2 t 2 2 t 3 2 EQ.4b d 1 = 2 t 0 t 1 EQ.4c Ω = 2 t 0 t 2 + 2 N t 3 EQ.4d 0 = 2 t 0 t 3 2 N t 2
Lemma 11. 
The Mu Part and Nu Part Equivalence.
EQ.5a 0 = 2 t 0 t 3 2 N t 2 EQ.5b 0 = t 0 t 3 N t 2 EQ.5c t 0 = N t 2 t 3 EQ.6a Ω = 2 t 0 t 2 + 2 N t 3 EQ.6b Ω 2 N t 3 = 2 t 0 t 2 EQ.6c t 0 = Ω 2 N t 3 2 t 2 EQ.7a N t 2 t 3 = Ω 2 N t 3 2 t 2 EQ.7b 2 N t 2 2 = Ω t 3 2 N t 3 2 EQ.7c t 2 2 = Ω t 3 2 N t 3 2 2 N .
Lemma 12. 
The Real Part and Nu Part Equivalence.
EQ.1a 0 = 2 t 0 t 3 2 N t 2 EQ.1b 0 = t 0 t 3 N t 2 EQ.2a t 2 = t 0 t 3 N EQ.2b Ω = 2 t 0 t 2 + 2 N t 3 EQ.2c Ω 2 N t 3 = 2 t 0 t 2 EQ.2d t 2 = Ω 2 N t 3 2 t 0 EQ.3a t 0 t 3 N = Ω 2 N t 3 2 t 0 EQ.3b 2 t 0 2 t 3 = Ω N 2 N 2 t 3 EQ.3c t 0 2 = Ω N 2 N 2 t 3 2 t 3 EQ.3d 1 t 0 2 = 2 t 3 Ω N 2 N 2 t 3
Lemma 13. 
The Lambda Part Identity
EQ.1 d 1 = 2 t 0 t 1 EQ.2 t 1 = d 1 2 t 0 EQ.3 t 1 2 = d 1 2 4 t 0 2 = d 1 2 4 2 t 3 Ω N 2 N 2 t 3 = 2 d 1 2 t 3 4 Ω N 8 N 2 t 3
Lemma 14. 
The Real Part Identity
EQ.1 d 0 = t 0 2 t 1 2 t 2 2 t 3 2 EQ.2 d 0 = t 0 2 d 1 2 4 t 0 2 t 2 2 t 3 2 EQ.3 d 0 = Ω N 2 N 2 t 3 2 t 3 2 d 1 2 t 3 4 Ω N 8 N 2 t 3 Ω t 3 2 N t 3 2 2 N t 3 2 EQ.4 d 0 = Ω N 2 N 2 t 3 2 t 3 2 d 1 2 t 3 4 Ω N 8 N 2 t 3 Ω t 3 2 N t 3 2 2 N + t 3 2 EQ.5 d 0 = Ω N 2 N 2 t 3 2 t 3 2 d 1 2 t 3 4 Ω N 8 N 2 t 3 Ω t 3 2 N EQ.6 d 0 = Ω N 2 2 N 3 t 3 Ω t 3 2 2 N t 3 2 d 1 2 t 3 4 Ω N 8 N 2 t 3 EQ.7 d 0 = Ω N 2 2 N 3 t 3 Ω t 3 2 2 N t 3 2 N t 3 2 N t 3 2 d 1 2 t 3 4 Ω N 8 N 2 t 3 EQ.8 d 0 = Ω N 2 2 N 3 t 3 Ω t 3 2 2 N t 3 4 N d 1 2 t 3 2 8 Ω N 2 t 3 16 N 3 t 3 2 EQ.9 d 0 = 4 Ω N 8 N 2 t 3 4 Ω N 8 N 2 t 3 Ω N 2 2 N 3 t 3 Ω t 3 2 2 N t 3 4 N d 1 2 t 3 2 8 Ω N 2 t 3 16 N 3 t 3 2 EQ.10 d 0 = 4 Ω 2 N 3 8 Ω N 4 t 3 4 Ω 2 N t 3 2 8 Ω N 4 t 3 + 16 N 5 t 3 2 + 8 Ω N 2 t 3 3 4 N d 1 2 t 3 2 8 Ω N 2 t 3 16 N 3 t 3 2 EQ.11 d 0 = + 8 Ω N 2 t 3 3 + 16 N 5 4 N d 1 2 4 Ω 2 N t 3 2 16 Ω N 4 t 3 + 4 Ω 2 N 3 8 Ω N 2 t 3 16 N 3 t 3 2 EQ.12 8 Ω N 2 d 0 t 3 16 N 3 d 0 t 3 2 = + 8 Ω N 2 t 3 3 + 16 N 5 4 N d 1 2 4 Ω 2 N t 3 2 16 Ω N 4 t 3 + 4 Ω 2 N 3 EQ.13 0 = 8 Ω N 2 t 3 3 + 16 N 5 4 N d 1 2 4 Ω 2 N t 3 2 16 Ω N 4 t 3 + 4 Ω 2 N 3 8 Ω N 2 d 0 t 3 + 16 N 3 d 0 t 3 2 EQ.14 0 = 8 Ω N 2 t 3 3 + 16 N 5 + 16 N 3 d 0 4 N d 1 2 4 Ω 2 N t 3 2 16 Ω N 4 + 8 Ω N 2 d 0 t 3 + 4 Ω 2 N 3 EQ.15 0 = 2 Ω N 2 t 3 3 + 4 N 5 + 4 N 3 d 0 N d 1 2 Ω 2 N t 3 2 4 Ω N 4 + 2 Ω N 2 d 0 t 3 + Ω 2 N 3 EQ.16 0 = 2 Ω N t 3 3 + 4 N 4 + 4 N 2 d 0 d 1 2 Ω 2 t 3 2 4 Ω N 3 + 2 Ω N d 0 t 3 + Ω 2 N 2 EQ.17 0 = A t 3 3 + B t 3 2 + C t 3 + D EQ.17a A = + 2 Ω N EQ.17b B = + 4 N 4 + 4 N 2 d 0 d 1 2 Ω 2 EQ.17c C = 4 Ω N 3 2 Ω N d 0 EQ.17d D = + Ω 2 N 2
Definition 15. 
Cardano’s Theorem: The Real Cubic Identity of the Nu Part
We now use the Cardano Method to depress the Cubic of the Nu Part.
EQ.1 0 = A t 3 3 + B t 3 2 + C t 3 + D EQ.2 0 = ζ 3 + p ζ + q EQ.3 ζ = t 3 + B 3 A EQ.4 p = 3 A C B 2 3 A 2 EQ.5 q = 2 B 3 9 A B C + 27 A 2 D 27 A 3 Completing Cardano’s depression of the Cubic. We now implement Vieta’s Substitution:
Definition 16. 
Vieta’s Theorem: The Resolution of the Nu Part.
EQ.1 ζ = w p 3 w EQ.2 0 = w 3 + q p 3 27 w 3 EQ.3 0 = w 6 + q w 3 p 3 27 EQ.4 y = w 3 EQ.5 0 = y 2 + q y p 3 27 EQ.6 y = q 2 ± q 2 4 + p 3 27 EQ.7 w = q 2 ± q 2 4 + p 3 27 3 , either sign of the square root shall suffice.
EQ.8 t 3 + B 3 A = w p 3 w EQ.9   t 3 = B 3 A + w p 3 w We now use the identities from t 3 to yield t 0 ,   t 1 and t 2 .
EQ.10 t 0 = ± Ω N 2 N 2 t 3 2 t 3 EQ.11 t 1 = d 1 2 t 0 EQ.12 t 2 = ± Ω t 3 2 N t 3 2 2 N Of course, we have a serious dilemma. Which of the three real roots do we accept for t 3 ? Which sign of the above squares do we choose in unison? Only the polar form of solution will elucidate which root of t 3 to accept, and then how to produce t 0 , t 1 and t 2 .
Theorem 17. 
The Offset Circle Theorem
For the moment, let us suppose we know which cubic root to select as t 3 , then we now examine the following relationship: N t 2 t 3 = Ω 2 N t 3 2 t 2 . This equation informs us that the coordinate t 2 μ 2 + t 3 ν 2 lays upon a circle, with a radius of Ω 4 N , offset from the origin by + Ω 4 N .
Let t 2 μ 2 + t 3 ν 2 = M μ 2 c o s θ + ν 2 s i n θ , then the Law of Cosines reveals that:
EQ.1 M 2 = Ω 4 N 2 + Ω 4 N 2 2 Ω 4 N 2 c o s 2 θ EQ.2 M 2 = 2 Ω 4 N 2 1 c o s 2 θ EQ.3 M = Ω 2 N s i n θ , upholding the Law of Sines.
We now examine the relationship t 0 = N t 2 t 3 .
EQ.4 t 0 = N M c o s θ M s i n θ EQ.5 t 0 = N c o t θ EQ.6 t 1 = d 1 2 t 0 = d 1 2 N t a n θ EQ.7 d 0 = t 0 2 t 1 2 t 2 2 t 3 2 EQ.8 d 0 = N 2 c o t 2 θ d 1 2 4 N 2 t a n 2 θ M 2 c o s 2 θ M 2 s i n 2 θ EQ.9 d 0 = N 2 c o t 2 θ d 1 2 4 N 2 t a n 2 θ M 2 c o s 2 θ + s i n 2 θ ; M 2 = Ω 2 4 N 2 s i n 2 θ ; c o s 2 θ + s i n 2 θ = 1 EQ.10 d 0 = N 2 c o t 2 θ 1 4 N 2 d 1 2 t a n 2 θ + Ω 2 s i n 2 θ , which leads to a nasty degree six equation with 6 pairs of conjugate solutions for θ .
Before we proceed, the below image is the geometric appearance of the question at hand in μ 2 , ν 2 space.
In the following url link, q is Omega, N is N, and t is theta: https://www.desmos.com/calculator/q2bfcbs7wq
Preprints 72263 i006
However, when we yield the roots of the cubic to produce t 3 we can solve for theta without any of the hassle that the polar form introduces.
EQ.11 t 2 μ 2 + t 3 ν 2 = M μ 2 c o s θ + ν 2 s i n θ EQ.12 t 3 = M s i n θ EQ.13 t 3 = Ω 2 N s i n 2 θ EQ.14 2 N t 3 Ω = s i n 2 θ ; 2 N t 3 Ω = s i n θ EQ.15 θ = A r c s i n + 2 N t 3 Ω . We know to take the positive root, since the t 2 , t 3 coordinate resides in the first quadrant, since both magnitude variables, N and Ω , are positive by definition, forcing the red circle (in the above image) in the upper two quadrants. We also both angles for the Arcsine function. This is not because we cannot resolve the ambiguity; rather, both θ solutions fulfill d = t 2 t v + v t simultaneously. Hence:
EQ.16 θ 1 = A r c s i n + 2 N t 3 Ω ;     θ 2 = π θ 1 ; yielding the empirically observed form: t = t 3 ν 2 ± Ψ 2 , where Ψ has no ν 2 part.
With both values of theta known, we simply use the identities above to yield t 0 ,   t 1 , t 2 EQ.17a t 0,1 = N c o t θ 1 ; t 1,1 = d 1 2 t 0,1 ; t 2,1 = Ω 2 N s i n θ 1 c o s θ 1 EQ.17b t 0,2 = N c o t θ 1 ; t 1,2 = d 1 2 t 0,1 ; t 2,2 = Ω 2 N s i n θ 2 c o s θ 2 ; t 0,1 = t 0,2 ; t 1,1 = t 1,2 ; t 2,1 = t 2,2 ; t 3,1 = t 3,2 . Q.E.D.
Theorem 18. 
The Six Root Equivalence Theorem
At this moment in the paper, there is no known method to determine which of the three cubic roots to choose for t 3 . However, a simple guess and check on the resulting t 2 , t 3 coordinates suffices, since the coordinate t 2 μ 2 + t 3 Ω 4 N ν 2 must have a magnitude of Ω 4 N . That is, given:
EQ.1 t 2,1 = Ω 2 N s i n θ 1 c o s θ 1 ; t 3 = Ω 2 N s i n 2 θ 1 , when does:
EQ.2 Ω 4 N = Ω 2 N s i n θ 1 c o s θ 1 2 + Ω 2 N s i n 2 θ 1 Ω 4 N 2 EQ.3 Ω 2 N s i n θ 1 c o s θ 1 2 = Ω 2 s i n 2 2 θ 1 16 N 2 EQ.4 Ω 2 N s i n 2 θ 1 Ω 4 N 2 = Ω 2 c o s 2 ( 2 θ 1 ) 16 N 2 EQ.5 Ω 4 N = Ω 2 s i n 2 2 θ 1 16 N 2 + Ω 2 c o s 2 2 θ 1 16 N 2 EQ.5 Ω 4 N = Ω 2 s i n 2 2 θ 1 + Ω 2 c o s 2 2 θ 1 16 N 2 EQ.6 Ω 4 N = 1 4 N Ω 2 s i n 2 2 θ 1 + c o s 2 2 θ 1 EQ.7 Ω 4 N = Ω 4 N c o s 2 2 θ 1 + s i n 2 2 θ 1 EQ.8 Ω 4 N = Ω 4 N 1 EQ.9 Ω 4 N = Ω 4 N Which informs us that all three cubic roots suffice, and thus a Quadratic Equation over the Quaternions has three pairs of roots that conjugate over the ν 2 307 axis, since all resultant Arcsine values yield the reflexive property Ω 4 N = Ω 4 N . Thus there exists six solutions to a Quaternionic Quadratic.
However, only one of the three real roots with return a Real Value of Theta for A r c s i n + 2 N t 3 Ω , while the other two roots shall return imaginary values for theta.
That does not mean however that we discard these solutions, since any resulting imaginary unit is equal to λ , but rather that it introduces non-real coefficients against the vector components, and culmination of these coefficients in the form of α q + β λ recombine when multiplied against the q ,   λ   ,   μ 2 and ν 2 vector components to yield the exact same value of t as the real coefficients that resulted from the t 3 root whose θ was also real.
That is, all three Arcsine arguments of t 3 produce the same t vector after recombination.
Q.E.D.
Theorem 19. 
The Closed Form Solution for a General Quadratic Equation for the Quaternions.
We now combine all of the steps to solve original query:
EQ.1 f = x 2 + x   b + a   x EQ.2 Let c = f + a   b EQ.3 c = x 2 + x   b + a   x + a   b EQ.4 u = 1 2 a + b EQ.5 v = 1 2 a b   EQ.6 Let t = x + u , therefore x = t u EQ.7 c = t 2 t   v + v   t v   2 EQ.8 Let       d = c + v 2 = d 0,0 q + d 1,0 i + d 2,0 j + d 3,0 k EQ.9 d = t 2 t   v + v   t EQ.10 v = v 0,0 q + v 1,0 i + v 2,0 j + v 3,0 k EQ.11 α = A T A N 2 v 2 v 1 EQ.12 β = A T A N 2 v 3 v 1 1 + v 2 2 EQ.13 λ = + i c o s α c o s β + j s i n α c o s β + k s i n β = γ 1,1 i + γ 1,2 j + γ 1,3 k EQ.14 μ 1 = i s i n α + j c o s α + 0 k = γ 2,1 i + γ 2,2 j + γ 2,3 k 329 EQ.15 ν 1 = i c o s α s i n β + j s i n α s i n β + k c o s β = γ 3,1 i + γ 3,2 j + γ 3,3 k EQ.16 Let 𝚪 be a 3x3 real matrix whose pairwise entries are equal to γ m , n .
EQ.17 Let A be a 1x3 real column matrix whose entries are v 1,0 ,   v 2,0 and v 3,0 respectively.
EQ.18 Let B be a 1x3 real column matrix whose entries are d 1,0 ,   d 2,0 and d 3,0 respectively.
EQ.19 Let 𝐕 =𝚪 A , which is also a 1x3 real column matrix, let it the results be named N , 0,0 . N is our first primary variable.
EQ.20 Let D =𝚪 B , which is also a 1x3 real column matrix, let it the results be named d 1,1 , d 2,1 , d 3,1 We do not require the inverse Gamma Matrix for this process.
𝚪 Gamma Matrix A Matrix B Matrix 𝚪 A = V Matrix 𝚪 B = D Matrix
Preprints 72263 i012
EQ.21 ϕ = A T A N 2 d 3,1 d 2,1 ; Ω = d 2,1 2 + d 3,1 2 EQ.22 μ 2 = + μ 1 c o s ϕ + ν 1 s i n ϕ = + i s i n α c o s ϕ c o s α s i n β s i n ϕ + j + c o s α c o s ϕ + s i n α s i n β s i n ϕ + k c o s β s i n ϕ EQ.23 μ 2 = μ 1 i + μ 2 j + μ 3 k EQ.24 ν 2 = μ 1 s i n ϕ + ν 1 c o s ϕ = + i + s i n α s i n ϕ c o s α s i n β c o s ϕ + j c o s α s i n ϕ + s i n α s i n β c o s ϕ + k c o s β c o s ϕ EQ.25 ν 2 = ν 1 i + ν 2 j + ν 3 k EQ.26 v = 1 2 a b = v 0,0 q + N λ + 0 μ 1 + 0 ν 1 = v 0,0 q + N λ + 0 μ 2 + 0 ν 2 EQ.27 d = d 0,0 q + d 1,1 λ + d 2,1 μ 1 + d 3,1 ν 1 = d 0,0 q + d 1,1 λ + Ω μ 2 + 0 ν 2 EQ.28 A = + 2 Ω N EQ.29 B = + 4 N 4 + 4 N 2 d 0,0 d 1,1 2 Ω 2 EQ.30 C = 4 Ω N 3 2 Ω N d 0,0 EQ.31 D = + Ω 2 N 2 EQ.32 0 = A t 3 3 + B t 3 2 + C t 3 + D EQ.33 0 = ζ 3 + p ζ + q EQ.34 ζ = t 3 + B 3 A EQ.35 p = 3 A C B 2 3 A 2 EQ.36 q = 2 B 3 9 A B C + 27 A 2 D 27 A 3 EQ.37 ζ = w p 3 w EQ.38 0 = w 3 + q p 3 27 w 3 EQ.39 0 = w 6 + q w 3 p 3 27 EQ.40 y = w 3 EQ.41 0 = y 2 + q y p 3 27 EQ.42 y = q 2 ± q 2 4 + p 3 27 EQ.43 w = q 2 ± q 2 4 + p 3 27 3 , either sign of the square root shall suffice, and any cube root will suffice.
EQ.44 t 3 + B 3 A = w p 3 w EQ.45   t 3 = B 3 A + w p 3 w . All roots will be real.
EQ.46 θ = A r c s i n + 2 N t 3 Ω . If θ is a complex number, then λ is the imaginary unit.
The complex arguments of θ can be used, substituting i for λ and then wielding these complex coefficients against the q , λ , μ 2 and ν 2 vector components of t . However, to avoid confusion, it is most practical check all three roots of t 3 and just select the real argument of θ . EQ.47 t 1 = + q N c o t θ + λ d 1 , , 1 2 N t a n θ + μ 2 Ω 2 N s i n θ c o s θ + ν 2 Ω 2 N s i n 2 θ = t 0,1 q + t 1,1 λ + t 2,1 μ 2 + t 3,1 ν 2 EQ.48 t 2 = q N c o t θ λ d 1 , , 1 2 N t a n θ μ 2 Ω 2 N s i n θ c o s θ + ν 2 Ω 2 N s i n 2 θ = t 0,2 q + t 1,2 λ + t 2,2 μ 2 + t 3,2 ν 2 The above two equations are the roots in the proper orthogonal basis of λ ,   μ 2 , ν 2 , however we must now convert back to i ,   j ,   k .
EQ.49 t 0,1 q + t 1,1 λ + t 2,1 μ 2 + t 3,1 ν 2 = t 0,3 q + t 1,3 i + t 2,3 j + t 3,3 k EQ.50 t 0,3 = t 0,1 EQ.51 t 1,3 = t 1,1 γ 1,1 + t 2,1 μ 1 + t 3,1 ν 1 EQ.52 t 2,3 = t 1,1 γ 1,2 + t 2,1 μ 2 + t 3,1 ν 2 EQ.53 t 2,3 = t 1,1 γ 1,3 + t 2,1 μ 3 + t 3,1 ν 3 EQ.54 t 0,2 q + t 1,2 λ + t 2,2 μ 2 + t 3,2 ν 2 = t 0,4 q + t 1,4 i + t 2,4 j + t 3,4 k EQ.55 t 0,4 = t 0,2 EQ.56 t 1,4 = t 1,2 γ 1,1 + t 2,2 μ 1 + t 3,2 ν 1 EQ.57 t 2,4 = t 1,2 γ 1,2 + t 2,2 μ 2 + t 3,2 ν 2 EQ.58 t 2,4 = t 1,2 γ 1,3 + t 2,2 μ 3 + t 3,2 ν 3 Recall that t = x + u and therefore x = t u and that u = 1 2 a + b .
EQ.59 t 1 = x 1 + u = t 0,3 q + t 1,3 i + t 2,3 j + t 3,3 k EQ.60 t 2 = x 2 + u = t 0,4 q + t 1,4 i + t 2,4 j + t 3,4 k EQ.61 u = u 0 q + u 1 i + u 2 j + u 3 k EQ.62 x 1 = t 0,3 u 0 q + t 1,3 u 1 i + t 2,3 u 2 j + t 3,3 u 3 k EQ.63 x 2 = t 0,4 u 0 q + t 1,4 u 1 i + t 2,4 u 2 j + t 3,4 u 3 k The above two equations satisfy the original query f = x 2 + x   b + a   x , proving that all Quadratic Equations over the Quaternions adhere to the same closed form solution.
Q.E.D.

Statements and Declarations

I have nothing to declare.

Appendix A. Overview of Cayley Geodesicity, Orientability, Chirality, Algebraicity

In Volume Two, when we move beyond the Quaternions, we shall give this topic proper treatment. For now, a summary of the definitions of Geodesicity, Orientability, Chirality and Algebracity shall suffice. Below are three Cayley Tables of four, six and twenty-two dimensional hypercomplex spaces that satisfy these four properties. We shall not address odd dimensional spaces until Volume Three.
1 . A Cayley Table is said to be Geodesic if each orthogonal imaginary unit vector appears once, and only once, in each row and column.
2 . A Cayley Table is said to be Orietnable if the square of each imaginary unit vector is equal to q and the square of ± q = + q .
3 . A Cayley Table is said to be Chiral if it is both Geodesic and Chiral, and the resulting sign of each Cayley product is the opposite of its transposed product for the Even Dimensional Spaces.
4 . A Cayley Table is said to be Algebraic if it is Geodesic, Oreitnable and Chiral, and that the resulting product signs adhere to uniform pattern over the product grid that allow one to calculate all m t h roots of n -unity for all hypercomplex vectors in that space.
We shall cover the specific pattern of sign placement that admits Algebracity in Volume Two. For now, see if you can discern the pattern of positive 91 and negative sign placement in the 22D table. The first dimension, " 1 " is q in the below table.
Algebraic Icosidgions 
Preprints 72263 i007
In all of the below tables, the dimension 1 is equal to q .
Algebraic Hexonions
Preprints 72263 i008
Algebraic Octonions
Preprints 72263 i009
Algebraic Sedenions
Preprints 72263 i010
I do not feign to know why the below paradigm works, I only know that it works (via copious brute force experiments). Only the below layout of dimension and sign placement allows one to take the Mth root of N unity (see Appendix B) for the Even Hypercomplex Complex Numbers and locate the roots of a quadratic equation.
Preprints 72263 i011

Appendix B. The M-th Root of N-Unity for Class of Algebraic Hypercomplex Numbers of Even Dimensions, The Great Circle Theorem

Assuming that we are in a Hypercomplex Space that is Cayley Algebraic, then let x n ,   m be the function that returns the m t h principal root unity for x .
EQ.11d Let q be the observation vector.
EQ.12d Let D be the set of pairwise orthogonal imaginary unit vectors, D = n and n must be odd, such that D       q   is even.
EQ.13d Let x = α 0 q + z = 1 z = n α z d z such that z ,   α z R EQ.14d Let β = + z = 1 z = n α z 2 , which is the real number magnitude of the imaginary part of x .
EQ.15d Let γ = + α 0 2 + z = 1 z = n α z 2 , which is the real number magnitude of x .
EQ.16d Let λ = 1 β x α 0 q   , compelling λ to be a unit vector.
EQ.17d Let τ = 1 β α 0 , giving us the ratio between the magnitudes of the real part and the imaginary part.
EQ.18d Let θ = A T A N 2 1 τ = A C O T A N 2 τ , that is, the four-quadrant arccotagent of τ .
EQ.19d x = γ q c o s θ + λ s i n θ EQ.20d x n , m = γ n q c o s θ + 2 π m n + λ s i n θ + 2 π m n , m , n Z , m n

Appendix B.2. Corollary: The Square Root of a Quaternion, The Well Defined Positive and Negative Square Roots

For a quaternion x = α 0 q + α 1 i + α 2 j + α 3 k   , the square root is given by:
EQ.21d
+ x = + α 0 2 + α 1 2 + α 2 2 + α 3 2 q c o s 0 + 1 2 A T A N 2 + α 1 2 + α 2 2 + α 3 2 α 0 + 1 + α 1 2 + α 2 2 + α 3 2 x α 0 q   s i n 0 + 1 2 A T A N 2 + α 1 2 + α 2 2 + α 3 2 α 0 x = + α 0 2 + α 1 2 + α 2 2 + α 3 2 q c o s π + 1 2 A T A N 2 + α 1 2 + α 2 2 + α 3 2 α 0 + 1 + α 1 2 + α 2 2 + α 3 2 x α 0 q   s i n π + 1 2 A T A N 2 + α 1 2 + α 2 2 + α 3 2 α 0

Appendix C. The Quadratic Equation , The Trivial Case of Symmetric Roots, For All Hypercomplex Dimensions

EQ.1e C = x 2 + x y + y x can be readily solved for x even if C and y are not on the same Great Circle.
EQ.2e Let t = x + y , such that x = t y EQ.3e C = t y 2 + t y y + y t y EQ.4e C = t 2 t y y t + y 2 + t y y   2 + y t y 2   EQ.5e C = t 2 y 2 EQ.6e t 2 = y 2 C EQ.7e t = ± y 2 C           x + y = ± y 2 C     EQ.8e x = y ± y 2 C   EQ.9e 0 = y ± y 2 C 2 + y ± y 2 C y + y y ± y 2 C + C Q.E.D
It is our goal to transform the earlier equation, c = w 2 + w y , into the Symmetric Case via a series of additional substitutions.
The Symmetric Case occurs when C = x 2 +   y   4 L +   y   4 R x .

References

  1. https://mathworld.wolfram.com/VietasSubstitution.html.
  2. https://www.math.ucdavis.edu/~kkreith/tutorials/sample.lesson/cardano.html.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2026 MDPI (Basel, Switzerland) unless otherwise stated