Implementation of Bestvina-Handel algorithm to find train tracks for surface homeomorphisms