Fast solution of toeplitz systems of equations and computation of Padé approximants