Added sage version of code by Mark van Hoeij to compute divisors of modular units