Load Doc/factor-a.sm1, then factorization can be computed by invoking asir, of which ftp cite is at http://www.math.kobe-u.ac.jp