The algorithm was implemented using the LiDIA library for computational number
theory and it will be part of the multiprecision floating-point arithmetic of
the package in release 1.4. LiDIA is available from
ftp://crypt1.cs.uni-sb.de/pub/systems/LiDIA/LiDIA-1.2.1.tgz
http://www-jb.cs.uni-sb.de/LiDIA/linkhtml/lidia/lidia.html
Here is the output of the program: