In previous papers we have outlined a program for deriving the infrared behavior of the axial gauge gluon propagator in a pure Yang-Mills theory. The program is based on an integral equation for the gluon propagator derived from the Dyson equation and the Ward identities. Here we present a solution to this equation, obtained numerically. The solution exhibits a Singularity in the infrared, and therefore presumably predicts confinement of color. The method is supposed to be exact in the infrared. Away from the infrared, therefore, our solution is only approximate. Nevertheless, even in the ultraviolet, our solution for the propagator is not very different from the known asymptotic freedom result, so it may be that it is a reasonable approximation over the entire range of momentum.