Shamelessly lifted from http://homepage.univie.ac.at/dietrich.burde/pari/aks.gp