see http://nswt.tuwien.ac.at:8000/htdocs/internet/unix/perl/math-pari.html