In their precedent work, the authors constructed closed oriented hyperbolic surfaces with pseudo-Anosov homeomorphisms from certain class of integral matrices. In this paper, we present a very simple algorithm to compute the Teichmüller polynomial corresponding to those surface homeomorphisms by first constructing an invariant track whose first homology group can be naturally identified with the first homology group of the surface, and computing its Alexander polynomial.