An algorithm producing hyperbolicity equations for a link complement in S^3