Reach Us
+44-1522-440391

Engineering Mathematics and Physics Department, Faculty of Engineering, Alexandria University, Egypt

- *Corresponding Author:
- Shebl SA

Engineering Mathematics and Physics Department

Faculty of Engineering

Alexandria University, Alexandria

2438, Egypt

Tel: 2035910052

E-mail: [email protected]

**Received** September 01, 2015; **Accepted** September 07, 2015; **Published** September 09, 2015

**Citation:** Shebl SA (2015) Semi Analytical Approach for Establishing the Two Real Roots of Sextic Polynomial. J Appl Computat Math 4:254. doi:10.4172/2168-9679.1000254

**Copyright:** © 2015 Shebl SA. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Applied & Computational Mathematics

Zeros of polynomial equations are analytically hard to be determined beyond the special cases of the quartic equations. Under some particular conditions, quintic and sextic polynomial equations may be solved iteratively. This paper presents analytical-graphical solution for accomplishing zeros of particular sextic polynomial equation of two real zeros. The concerned polynomial has been modeled geometrically as the radical problem in geodesy which is the geodetic height of a point on the terrain surface of the earth. The earth’s model to be adopted is the triaxial ellipsoidal surface. The achieved solution may be utilized as initial values for a convenient and convergent iterative process.

Sextic; Triaxial ellipsoid; Mongean projection; Geodetic height; Umbilic points

For a long time, solutions of polynomial equations of higher order have been dealt with either graphically or analytically [1]. The works of Abel and Galois have shown that the general polynomial equations of degree higher than the fourth cannot be solved in radicals [2]. While Abel published the proof of impossibility of solving these equations (Abel's impossibility theorem), Galois gave a more rigorous proof using the group theory. This does not mean that there is no algebraic solution to the general polynomial equations of degree five and above [3]. In fact these equations are solved algebraically by employing symbolic coefficients: the general quintic is solved by using the Bring radicals, while the general sextic can be solved in terms of Kampe de Feriet functions [4].

Under certain limited conditions concerning the coefficients and roots, Kulkarni described a method to decompose the given sextic equation into two cubic polynomials as factors. The cubic polynomials are then equated to zero and solved to obtain the six roots of the sextic equation in radicals [5].

This paper introduces a straightforward and accurate solution for sextic polynomial equation of two real roots involving an analyticalgraphical analysis. The related polynomial characterizes the geometrical model of a common problem in geodesy which is the geodetic height [6]. Therefore, we pursue the solution of sextic polynomial equation of two real roots originated accordingly from the geodetic height problem.

Such a model symbolizes the coordinates of two foot points of the two normals from a point on the irregular surface of the earth (terrain) to the triaxial ellipsoidal surface as the earth’s model [7]. This has been yielded two curves characterize the locii of the traces of the normals from the terrain point to the surface circular sections whose intersections are the desired two foot points.

The graphical solution was granted and considered to identify such locii and, consequently, their intersections whose coordinates may be exploited as a preliminary estimate of initial values when iterative process is considered.

The sextic polynomial of interest takes the form:

(1)

and, according to **Figure 1**, considering a triaxial ellipsoid surface’s center as the coordinate system’s origin, aligning z-axis to the polar axis direction and x- axis to the east direction, then equation of the surface may be expressed as:

(X^{2}/a^{2}) + (Y^{2}/b^{2}) + (Z^{2}/c^{2}) = 1 (2)

Such reference surface has dimensions (a) along the semi-major axis (x-axis), (b) along the semi-minor axis (y-axis) and (c) along the polar axis (z-axis), and a>b>c (**Figure 1**).

There are six normals from E to the reference surface [6-8], four of which are imaginary and two are real. The two normals from the terrain point E(X_{E}, Y_{E}, Z_{E}), **Figure 2**, to the reference surface has two foot points P of coordinates (X_{P}, Y_{P}, Z_{P}), one of them gives the minimum distance from E to the surface, i. e. the geodetic height [6].

The evaluation of the model relies absolutely on the location of the terrain point E and the size of the reference surface as well (**Figure 2**).

The foot points P may be located graphically as the points of intersection of two curves μ and η, each of which is the locus of the normals from E to a set of circular sections of the surface. The technique traced for constructing μ, η. Coordinates of P is utilized as follows.

In **Figure 2**, PE is the normal from E to the surface, and then the parametric coordinates of P are:

X_{P} = X_{E}/t

Y_{P} = Y_{E}(1−e_{1}^{2})/(t –e_{1}^{2})

Z_{P} = Z_{E}(1−e^{2}_{2})/(t –e^{2}_{2}) (3)

Where, when applying "3" to "2", t stands for the roots of "1", e_{1} is the equatorial eccentricity, and e^{2}_{2} is the polar eccentricity, where:

e_{1}^{2} = 1−(b^{2}/a^{2}) (4)

and:

e^{2}_{2} = 1− (c^{2}/a^{2}) (5)

Substituting these coordinates into the equation of the surface (Equation "2") and comparing the coefficients of the resulting equation with the coefficients of "1", then:

a_{6} = 1 (6)

a_{5} = −2e_{1}^{2} −2e^{2}_{2} (7)

a_{4} = − F^{2} + e^{2}_{4} + e^{1}_{4} + 4e_{1}^{2}e^{2}_{2}– G^{2} − e_{1}^{2}G^{2} − H^{2} + e^{2}_{2}H^{2} (8)

a_{3} = 2e^{2}_{2}F^{2} + 2e_{1}^{2}F^{2} − 2e^{1}e_{2}^{4} − 2e_{1}4e^{2}_{2} + 2e^{2}_{2}G^{2} − 2e_{1}^{2}e^{2}_{2}G^{2} + 2e_{1}^{2}H^{2}− 2e_{1}^{2}e^{2}_{2}H^{2} (9)

a_{2}=−e_{2}^{4}F^{2} – e_{1}^{4}F^{2} − 4e_{1}^{2}e^{2}_{2}F^{2} + e_{1}^{4}e_{2}^{4} − e_{2}^{4}G^{2} + e_{1}^{2}e_{2}^{4}G^{2} − e_{1}^{4} H^{2}+ e1^{4}e_{2}^{2}H^{2} (10)

a_{1} = 2e_{1}^{2}e_{2}^{4}F^{2} + 2e_{1}^{4}e_{2}^{4} F^{2} (11)

a_{0} = − e_{1}^{4}e_{2}^{4}F^{2} (12)

where:

F^{2} = X_{E}^{2}/a^{2} (13)

G^{2} = YE^{2}/b^{2 } (14)

H^{2}= ZE^{2}/c^{2} (15)

Solving "6" to "15" in terms of a_{0}, a_{1}, a_{3}, a_{4} and a_{5} yields:

e_{1}^{2}= −(a_{5}/4) – ε (16)

e_{2}^{2} = −(a_{5}/4) + ε (17)

F^{2} = − a_{1}^{2}/(a_{0}a_{5}^{2}) (18)

G^{2} = [α / (4(1 − e_{1}^{2}))] + [β / (4ε (1− e_{1}^{2})) (19)

H^{2} = [α / (4(1 – e_{2}^{2}))] − [β / (4ε (1− e_{2}^{2}))] (20)

where:

ε = [(a_{5}^{2} /16) – (a_{0}a_{5}/a_{1})]1/2 (21)

α = (2a_{1}^{2} /(a_{0} a_{5}^{2})) + (a_{5}^{2} / 2) + (4 a_{0} a_{5}/a_{1}) – 2 a_{4} (22)

β = − (a_{4}a_{5}/2) + a_{3} – (a_{1}^{2} / (2a_{0}a_{5})) + a_{5}^{3}/8 (23)

Assuming a reasonable value for a, then values of b and c can be determined in terms of the polynomial coefficients from "16" and "17", as:

b = a (1+ ε + a_{5}/4)1/2 (24)

c = a (1 − ε + a_{5}/4)1/2 (25)

Also, values of X_{E}, Y_{E}, and Z_{E} can be obtained in terms of the polynomial coefficients from "13" through "15".

Now, to obtain the desired roots t (from Equation "3"), we seek the foot point P(X_{P},Y_{P} ,Z_{P}) which, as previously mentioned, may be located graphically as the points of intersection of two curves μ and η, each of which is the locus of points of intersection between normals from E to a set of circular sections, with these sections. The technique traced for constructing μ, η and coordinates of P is utilized as follows.

**Figure 3 **displays a significant aspect associated with the triaxial ellipsoidal surface that is the circular sections. There exist two sets of such sections [8]. The first set of these sections is generated as the intersections of the surface with a set of parallel planes whose equations are (**Figure 3**):

Z = (X – δ )[c^{2} (a^{2} – b^{2}) / (a^{2}(b^{2} – c^{2}))]1/2 (26)

where δ is the length along x -axis traced by the planes (**Figure 4**). These planes are parallel to planes tangent to the surface at the umbilic points U_{1} and U_{2} whose coordinates are (**Figure 4**):

(U_{1}, U_{2}) = (± a [(a^{2} – b^{2}) / (a^{2}– c^{2})]1/2, 0,c [(b^{2} – c^{2}) / (a^{2}– c^{2})]1/2) (27)

The locus of these centers is a line joining U_{1} and U_{2} and its equation is:

Z = − X [c^{2} (b^{2} – c^{2}) / (a^{2}(a^{2} – b^{2}))]1/2 (28)

and the centers of the circles are:

C_{1} = (± δ (a^{2} – b^{2}) / (a^{2} – c^{2}), 0,δ c (a^{2} – b^{2})1/2(b^{2} – c^{2})1/2/ (a(a^{2} – c^{2}))) (29)

Similarly, the other set of the circular sections are generated by intersecting the surface with the parallel planes of equations:

Z = − (X – δ)[c^{2} (a^{2} – b^{2}) / (a^{2}(b^{2} – c^{2}))]1/2 (30)

whose corresponding umbilics U_{3} and U_{4} are of coordinates:

(U_{3}, U_{4}) = (± a [(a_{2} – b^{2}) / (a_{2}– c^{2})]1/2, 0, ± c [(b^{2} – c^{2})/ (a_{2}– c^{2})]1/2) (31)

The locus of these centers is a line joining U_{3} and U_{4} whose equation is:

Z = X [c^{2}(b^{2} – c^{2})/(a_{2}(a_{2} – b^{2}))]1/2 (32)

and the centers of the circles are:

c^{2} = (± δ (a_{2} – b^{2})/(a_{2} – c^{2}), 0, ± δ c(a_{2} – b^{2})1/2(b^{2} – c^{2})1/2/(a(a_{2} – c^{2}))) (33)

Radii of the circles (r), in terms of (δ), are:

r = b [1 – (δ^{2}(a^{2} – b^{2})/a^{2}(a^{2} – c^{2}))]1/2 (34)

Employing these circular sections, the graphical representation of the model produces the foot points P_{1} and P_{2}. This may be practiced using the auxiliary orthogonal projection (Mongean Projection) of the circular sections onto an auxiliary plane parallel to the circles' planes (**Figure 5**). In such projection, the circles are projected as circles and, consequently, the projection of any line normal to any circle is projected as the line joining the projection of E (E_{1} or E_{2}) and the projection of the center of the circle. The line joins points of intersection of the normals to the circles are a curve μ and η (**Figure 5**).

Obviously, intersection of the two locii of all normals to the two sets of circular sections, i. e., μ and η, yields P_{1} and P_{2} and, consequently, EP_{1} and EP_{2} are the normals from E to the surface. Substituting coordinates of P_{1} and P_{2} in "3" gives the two real roots of the polynomial. Either P_{1} or P_{2} is related to the minimum distance from the terrain point E to the surface, i.e. the geodetic height.

The above illustrated analysis is applicable only for the desired conditions of "1". These conditions arose according to the mutual relation between size of the surface and the polynomial coefficients.

These conditions may be formulated from "6" to "12" as:

a_{6} = 1

−4 < a_{5} < 0

a_{1} > 0

a_{0} < 0

a_{0}/a_{1} > a_{5} /16

−1 −a_{5}/2 < a_{0}a_{5}/a_{1} < a_{5}^{2}/16

a_{4} < a_{1}^{2}/(a_{0}a_{5}^{2}) + (a_{5}^{2}/4) + (2a_{0}a_{5}/a_{1})

a^{2} = (3a_{1}/a_{5}) + (3a_{0}^{2}a_{5}2/a_{1}^{2}) – (a_{0}a_{4}a_{5}/a_{1}) + (a_{3}a_{5}/4) (35)

Now, we conclude the process as:

1) Verifying conditions (Equation "35").

2) Supposing a reasonable value for a.

3) Calculating ε; b; c; F; G and H from "21"; "16"; "17"; "18"; "19" and "20" respectively.

4) Calculating X_{E}; Y_{E} and Z_{E} from "13"; "14" and "15" respectively.

5) Constructing Mongean Projection and getting the coordinates of the trace P_{1} and P_{2}.

6) Substituting coordinates of P_{1} and P_{2} in "3" to get the two real roots of "1".

We introduce a numerical example for the purpose of assessing the applicability of the procedure and verifying its accuracy. Consider the polynomial equation:

t^{6} − 1.74 t^{5} – 1.2695 t^{4} + 3.064728 t^{3} – 1.7921016 t^{2} + 0.46002816 t – 0.048540902 = 0

The first step is checking the coefficients whether they satisfy conditions in "35" or not, where:

a_{6} = 1

−4 < (a_{5}=- 1.74) < 0

(a_{1} = 0.46002818) > 0

(a_{0} = −0.0485409024) < 0

(a_{0}/a_{1} = −0.1055172405) < (a_{5} /16 = −0.10875)

(−1 − a_{5}/2 = −0.13)<(a_{0}a_{5}/a_{1} = 0.1835973508)< (a_{5}^{2}/16 = 0.189225)

(a_{4} = – 1.2695) < ( a_{1}^{2}/( a_{0}a_{5}^{2}) + (a_{5}^{2}/4) + (2a_{0}a_{5}/a_{1})= −0.4995000133)

a^{2} = (3a_{1}/a_{5}) + (3a_{0}^{2}a_{5}^{2}/a_{1}^{2}) – (a_{0}a_{4}a_{5}/a_{1}) + (a_{3}a_{5} /4) = – 1.7921016

Obviously, the coefficients satisfy the conditions, so we move to the second and third step, as we adopt 10 units for the value of a, and using "21"; "16"; "17"; "18"; "19" and "20", then:

ε = 0.075

b = 8

c = 7

Then the fourth step is calculating the terrain point E(X_{E}, Y_{E}, Z_{E}) using "13"; "14" and "15" respectively, as:

X_{E} = 12

Y_{E} = 10

Z_{E} = 8

The next step is the construction of Mongean graphical technique (**Figure 5**), in order to obtain the coordinates of the trace P_{1} and P_{2}. Using the values calculated above, the elements could be plotted using a triaxial ellipsoid model of dimensions a, b and c equal 10, 8 and 7 units respectively. The technique produces two locii μ and η which intersect at P_{1} and P_{2}, whose coordinates are:

P_{1} (6.84, 4.59, 3.51)

P_{2} (−8.45,−3.60, −2.03)

The final step is the determination of the desired real roots(t_{1}, t_{2}) of the proposed polynomial by substituting coordinates of P_{1} and P_{2} into "3", then we obtain: t_{1}=1.7540406851 and

t_{2}= −1.4202492291

Many geodetic applications engage high accuracy as in analyzing the earth's crustal movements, militant activities, satellite observations, etc. [6,7].

Thus, the significance of the proposed analysis is evaluated versus accuracy, as we apply the values of both t_{1} and t_{2} to the proposed sextic polynomial. The outputs are errors of values 4.1229×10^{-8} and 3.3855×10^{-7}which may be accepted in the favor of the graphical solution.

Finally, the solution enables getting the shortest distance EP_{1} from E to the surface (geodetic height), (**Figure 5**) which equals 8.7208829828 units.

Polynomial equations of degree higher than the fourth cannot be solved in radicals [2] except Kulkarni decomposition approach, furthermore, sextic polynomials are barely solvable analytically [5].

As the numerical methods are considered as the regular approach for handling the solvable sextic polynomials, they lack the initial guess required for setting up the solution. So, algorithm has to be adapted in order to isolate intervals for the real roots, roughly approximated values of the real roots are calculated then iterations are dealt with using appropriate convergent iteration algorithm, such as Newton method, in order to calculate such roots [9,10]. Furthermore, Newton's method for solving higher order polynomials will always converge if the initial point is sufficiently close to the root and if this root is not singular [11].

While Kulkarni introduced an analytical approach, yet the method is limited for the reducible sextic over the real field by the mean of decomposition [12].

Although the offered analysis conveys graphical involvement, it may be considered simple and straightforward method for computing the two real roots of the sextic polynomial. The results show an acceptable accuracy since the errors yielded are of values 4.1229×10^{-8} and 3.3855×10^{-7} for each root.

In addition, the analysis attends the geodetic height or the shortest distance from a given point to the triaxial ellipsoid surface (as an earth's model) avoiding the engagement of solving system of nonlinear equations using iteration solution scheme produced by Bektas [13].

Under certain conditions, a sextic polynomial equation of two real roots can be modeled as a familiar geodetic problem which is the geodetic height determination. The triaxial ellipsoidal surface is considered as the surface represents the earth’s irregular surface. Surface’s parameters and position of the terrain point are evaluated in terms of the polynomial’s coefficients, and vice versa.

Graphical solution is carried out for determining two foot points of the normal from a given point to the surface. Such points involve the desired roots and one of them associates the geodetic height.

Numerical justification for examining the approach was introduced and showed a reasonable accuracy versus the graphical solution.

- Dickenstein A,Emiris LZ (2005) Solving Polynomial Equations-Foundations,Algorithms and Applications. Springer verlag,Germany.
- King BR (1996) Beyond the Quartic Equation. Birkhauser, Boston.
- Kulkarni RG (2006) A versatile technique for solving quintic equations. Math ComputEduc 40: 205-215.
- Weisstein EW (2015) Sexticequation, WolframMathWorld.
- Kulkarni RG(2008) Solving sextic equations. Atlantic Electr J Math 3: 56-60.
- Shebl SA (1995)Conformalmapping ofthe triaxial ellipsoid and its applications in geodesy. Department ofEngineering Mathematics and physics, Faculty of Engineering, AlexandriaUniversity, Egypt.
- Snyder JP (1985) Conformal mapping of the triaxialellipsoid. Survey Rev 28: 130-148.
- Dresden A (1946)Solid Analytical Geometry andDeterminations. John Wiley and Sons, NewYork.
- Wang L, WankK, Zhou Z, Cui C(2011) Numerical methodsfor solving high order polynomial equations. Proceedings of the 4thInternational Conference on Information and Computing, Phuket Island, Thailand 150-153.
- Kalban K, El-Hajj A, Khoury S,Yousuf F(2002) Root Computations of Real-coefficient Polynomials usingSpreadsheets. Int J Engng, Great Britain18: 89-97.
- Sebah P,Gourdon X (2001) Newton’s Methodand High Order Iterations.
- Kulkarni RG (2008) Sextics by Division Method. The Teaching of Mathematics 11: 93-96.
- Bekats S (2014) Shortest Distance from a Point to TriaxialEllipsoid. International Journal ofEngineering and Applied Sciences 4.

Select your language of interest to view the total content in your interested language

- Adomian Decomposition Method
- Algebraic Geometry
- Analytical Geometry
- Applied Mathematics
- Axioms
- Balance Law
- Behaviometrics
- Big Data Analytics
- Binary and Non-normal Continuous Data
- Binomial Regression
- Biometrics
- Biostatistics methods
- Clinical Trail
- Complex Analysis
- Computational Model
- Convection Diffusion Equations
- Cross-Covariance and Cross-Correlation
- Differential Equations
- Differential Transform Method
- Fourier Analysis
- Fuzzy Boundary Value
- Fuzzy Environments
- Fuzzy Quasi-Metric Space
- Genetic Linkage
- Hamilton Mechanics
- Hypothesis Testing
- Integrated Analysis
- Integration
- Large-scale Survey Data
- Matrix
- Microarray Studies
- Mixed Initial-boundary Value
- Molecular Modelling
- Multivariate-Normal Model
- Noether's theorem
- Non rigid Image Registration
- Nonlinear Differential Equations
- Number Theory
- Numerical Solutions
- Physical Mathematics
- Quantum Mechanics
- Quantum electrodynamics
- Quasilinear Hyperbolic Systems
- Regressions
- Relativity
- Riemannian Geometry
- Robust Method
- Semi Analytical-Solution
- Sensitivity Analysis
- Smooth Complexities
- Soft biometrics
- Spatial Gaussian Markov Random Fields
- Statistical Methods
- Theoretical Physics
- Theory of Mathematical Modeling
- Three Dimensional Steady State
- Topology
- mirror symmetry
- vector bundle

- Total views:
**12402** - [From(publication date):

October-2015 - Jun 19, 2019] - Breakdown by view type
- HTML page views :
**8589** - PDF downloads :
**3813**

**Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals**

International Conferences 2019-20