In this article we present the methodology, according to which it is possible to derive approximate solutions for the roots of the general sextic polynomial equation as well as some other forms of sextic polynomial equations that normally cannot be solved by radicals; the approximate roots can be expressed in terms of polynomial coefficients. This methodology is a combination of two methods. The first part of the procedure pertains to the reduction of a general sextic equation H(x) to a depressed equation G(y), followed by the determination of solutions by radicals of G(y) which does not include a quintic term, provided that the fixed term of the equation depends on its other coefficients. The second method is a continuation of the first and pertains to the numerical correlation of the roots and the fixed term of a given sextic polynomial P(x) with the radicals and the fixed term of the sextic polynomial Q(x), where the two polynomials P(x) and Q(x) have the same coefficients except for the fixed term which might be different. From the application of the methodology presented above, the following formulation is derived; For any given general sextic polynomial equation P with coefficients within the interval [a, b], a defined polynomial equation Q corresponds which has equal coefficients to P except for its fixed term which might be different and dependent on the other coefficients so that Q has radical solutions. If we assume a pair of equations P, Q with coefficients within a predetermined interval [a, b], the numerical correlation through regression analysis of the radicals of Q, the roots of P and the fixed terms of P, Q, leads to the derivation of a mathematical model for the approximate estimation of the roots of sextic equations whose coefficients belong to the interval [a, b].