Accurate algebraic formula for the quintic & Solution by iteration of radicals Abdel Missa, Chrif Youssfi January 15, 2021 arXiv:2101.05610v1 [math.GM] 12 Jan 2021 1 Contents 1 Introduction 3 2 Reduced Forms 4 3 Iteration of radicals algorithm to solve equations in Form 3 5 3.1 Algorithm ............................. 5 3.2 Unique root near the positive real axis . 7 3.3 Absoluteerrortheorem. .. .. 8 3.3.1 Starting point properties . 9 3.3.2 Proofofabsoluteerrortheorem . 11 3.4 Relativeerrortheorem . 20 3.5 Speedofconvergence . 21 4 Solving equations in Form 2 24 5 Solving equations in Form 1 25 5.1 Challenge with iterative radicals . 25 5.2 Proposedalgorithm........................ 26 6 Trigonometric algorithm 30 7 Examples 31 7.1 Example 1: x5 + x + 0.01 = 0 .................. 32 7.1.1 Trigonometric algorithm . 32 7.1.2 Iteration of radicals algorithm . 33 7.2 Example 2: x5 + x +(3.08 + 1.68i)= 0 ............ 33 7.2.1 Trigonometric algorithm . 33 7.2.2 Iteration of radicals algorithm . 34 8 Future direction of work 35 2 Abstract According to the Abel-Ruffini theorem [1] and Galois theory [2], there is no solution in finite radicals to the general quintic equation. This article takes a different approach and proposes a new method to solve the quintic by it- eration of radicals. But, the most intriguing result is an accurate algebraic formula for absolute and relative root approximation: 3 2 formula- root < 4.32 10− and formula/root -1 < 2.51 10− . We then | | × | | × expand some of the geometric properties discussed to construct a trigono- metric algorithm that derives all roots. 1 Introduction Solving polynomial equations occupies a special place in the history of math- ematics. The first algebraic solution to the quadratic equation is attributed to Al-Khwarizmi [3] in the 9th century. In the 16th century, a solution to the cubic was found by Del Ferro, Tartaglia and Cardano [4] and was arguably the first significant discovery of the European Renaissance. Soon thereafter, Ferrari created a method to solve the general quartic equation [5]. The success achieved in finding formulas to lower degree equations stimulated a search to achieve the same ambition for the quintic. Despite efforts by Euler, Gauss, Lagrange and others, no general solutions were found. Finally, al- most 300 years after solving the cubic, Abel provided a complete proof of the impossibility of solving general quintic polynomials by finite combination of arithmetic operations and radicals. Subsequently, Galois went a step further and provided an exact criterion that characterizes solvable equations: their Galois group must be solvable. Naturally, the quintic equation required more tools. The first significant de- velopment was the discovery of the Bring-Jerrard normal form x5 +ux+v =0 by using respectively quadratic and quartic Tschirnhaus tranformations [6]. In 1858, Hermite [7] used a simplified normal form to provide the first known solution to general quintic equations in terms of elliptic modular functions. Shortly afterwards, Brioschi [8] and Kronecker [9] published equivalent solu- tions. In 1860, Cockle [10] and Harley [11] developed a method for solving the quintic using differential equations, leading to a solution in hypergeo- metric functions. And in 1884, Klein [12] provided an icosahedral solution of 3 the quintic. More recently, Doyle and McMullen [13] solved the quintic using an iterative method, while Glasser [14] provided a derivation to trinomial equations that leads to the same hypergeometric solution for the quintic as Cockle and Harley. In this article, we explore a new method that solves the quintic by itera- tion of radicals. Perhaps the most astonishing finding is the ability to find an accurate global approximation in radicals to one the five roots from the 3 first iteration: 1st iteration- root < 4.32 10− . We then proceed by proving the speed of convergence| of the proposed| iterative algorithm. Next, we discuss Bring radicals. Recall that the Bring radical of a com- plex number a is any of the five roots of: x5 + x + a = 0. More importantly, quintic equations can be solved in closed form using radicals and Bring rad- 5 icals. One might think that a simple iterative algorithm xk+1 = √a + xk − can solve the Bring equation. However, this is not the case as we show a counter example (when a is real) where the algorithm does not converge re- gardless of the starting point’s choice outside of the real root itself. On other hand, we prove that our proposed iterative algorithm converges globally for any complex number a. Furthermore, we provide a radical approximation to 2 one of the Bring radicals: BR(z)–approximation < 2.90 10− . | | Finally, we expand some of the properties used in the proposed algorithm to provide the location of the five roots in the complex plane. This leads to a trigonometric bisection algorithm that solves general quintic equations. 2 Reduced Forms The general quintic equation can be transformed to the Bring–Jerrard normal form: 5 v + d1v + d0 =0 From hereon, we define the n-th root of a complex number as the one with its argument in [ π/n, π/n[. − Excluding the trivial cases of d1 = 0 and d0 = 0, a change of variable 4 4 x = v/√d1 can further simplify the equation to Form 1: d0 x5 + x + a = 0 where a = (a =0) (1) 4 5 d1 6 p Another change of variable z = a/x leads to Form 2: z5 + z4 a4 = λ where λ = (2) 2 − 2 Let y = uz, (2) becomes: y5 + uy4 = u5λ 2 Choose u to satisfy: 1 λ λ 5 π π u5 = | | or u = | | = eiθ θ [ , [. λ λ ∈ − 5 5 This leads to Form 3: 5 4 y + uy + = ξ ξ = λ R ∗ (3) 2 | | ∈ Think of Form 3 as a simple rotation of Form 2 with the benefit of solving the equation close to the real axis. Notice that s is solution of (3) if and only if its conjugates ¯ is a solution of y5 +¯uy4 = ξ 2 Therefore we can assume from hereon that 0 θ π/5. ≤ ≤ 3 Iteration of radicals algorithm to solve equa- tions in Form 3 3.1 Algorithm Let us consider Equation (3) 5 4 y + uy + = ξ ξ R ∗ (3) 2 ∈ 5 π u = eiθ θ 0, ∈ 5 It can be exploited in two ways. First, by factorizing:h i 2ξ y = 4 (4) su + y Second, by completing the quintic: u 5 2u2 2u3 u4 u5 y + =2ξ + y3 + y2 + y + (5) 5 5 25 125 3125 By substituting y in the right-hand side of (5): y = G(ξ,y) (6) with: 2 3 4 5 5 2u 2u u u u G (ξ,y)= 2ξ + t3 + t2 + t + (7) r 5 25 125 3125 − 5 2ξ t = 4 (8) su + y This inspires the following fixed point algorithm with the following sequence: (yk)k N by: ∈ Starting point1 2 ξ 9 2+ √2 y0 = with α = (9) α s 4 Iterative process yk+1 = G (ξ,yk) (10) To prove the convergence and the accuracy of this algorithm, we proceed with the following steps: • Section 3.2 shows that there exists a unique root y∗ for Equation (3) near the positive real axis. 1The intuition behind α is discussed below. 6 + • Section 3.3 shows that ξ R ∗ and θ [0,π/5]: ∀ ∈ ∀ ∈ 3 y1 y∗ < 4.32 10− | − | In other words, the closed form formula: 2 3 4 5 5 2u 3 2u 2 u u u y1 = 2ξ + t + t + t + r 5 25 125 3125 − 5 with 2ξ 4 t = 2 v ξ 9 uu + α u t provides an accurate absolute approximation of y∗ • Section 3.4 shows that the relative error of the formula is also small: ξ R+ and θ [0,π/5] ∀ ∈ ∀ ∈ y1 2 1 < 2.58 10− y − ∗ • Finally, section 3.5 shows the sequence (yk)k N converges toward y∗ ∈ and: 1 k N yk+1 y∗ < yk y∗ K 15.44 ∀ ∈ | − | K | − | ≈ 3.2 Unique root near the positive real axis Consider equation (3). When θ = 0, u = 1 and it is obvious that there is a unique real solution that satisfies (3). Now assume θ > 0. We will show in this section that there exists a unique root y of Equation (3) near the real axis: θ y = reiσ with r> 0 and σ , 0 (11) ∈ −4 The real and the imaginary parts of the equation are: r5 cos(5σ)+ r4 cos(θ +4σ)=2ξ (12) and r5 sin (5σ)+ r4 sin (θ +4σ)=0 (13) 7 From (13): sin (θ +4σ) r = (14) − sin (5σ) Substituting r in (12): sin4 (θ +4σ) sin (σ θ) f(σ)= − =2ξ (15) sin5 (5σ) Let’s write: f(σ)= g1(σ)g2(σ) Where sin4 (θ +4σ) sin (σ θ) g1(σ)= and g2(σ)= − sin4 (5σ) sin (5σ) θ g1 is positive and stricly increasing in [ , 0[. Likewise g2 is positive and − 4 stricly increasing since: sin(4σ + θ) + 4 sin(θ σ) cos(5σ) g′ (σ)= − > 0 2 sin2 (5σ) Therefore f is continuous and strictly increasing. In addition: θ f = 0 and lim f (σ)=+ (16) −4 σ 0− ∞ → Consequently, there exists a unique σ [ θ , 0[ satisfying (15) or a unique ∈ − 4 root y∗ such as: θ Arg(y∗) < 0 (17) − 4 ≤ 3.3 Absolute error theorem + In this section, we will prove the absolute error theorem: ξ R ∗ and θ [0,π/5]: ∀ ∈ ∀ ∈ 3 y1 y∗ < 4.32 10− | − | We will start by discussing the unique properties of the selected starting point. 8 3.3.1 Starting point properties • P.1 y∗ y0 (18) | |≤| | • P.2 u + y∗ u + y0 (19) | |≤| | • P.3 Define 4 u + y0 h0 = u + y r ∗ π Arg (h0) (20) | | ≤ 20 Proof: • P.1: y∗ y0 | |≤| | Using (3) 8 2 2 y∗ u + y∗ =4ξ | | | | Therefore 8 2 2 y∗ 1+ y∗ +2 y∗ cos(θ σ) =4ξ | | | | | | − 2 Since 1 + y∗ 2 y∗ : | | ≥ | | 1 2 9 9 2 θ σ 2 ξ y∗ cos − ξ or y∗ 2 θ σ | | 2 ≤ | | ≤ cos − 2 ! Also: π θ 0 θ and σ 0 ≤ ≤ 5 − 4 ≤ ≤ Which leads to: θ σ π 0 − ≤ 2 ≤ 8 and θ σ π 2+ √2 cos − cos = = α 2 ≥ 8 s 4 It follows that: 2 ξ 9 y∗ = y0 | | ≤ α | | 9 • P.2: u + y∗ u + y0 | |≤| | Since u = 1 and the angle between u and y∗ is θ σ: | | − 2 2 u + y∗ =1+ y∗ +2 y∗ cos(θ σ) | | | | | | − Likewise, since the angle between u and y0 is θ: 2 2 u + y0 =1+ y0 +2 y0 cos(θ) | | | | | | From P.1: y∗ y0 .