An Initial Solution Algorithm for Globally Convergent Homotopy Methods