This code is quick and dirty.  Feel free to modify it as you
wish.

If you download the precompiled Windows binary from the GP-PARI website,
you will need CYGWIN1.DLL zipped up in this directory. Place all three
text files in a directory GP-PARI can read.  Change the seed for the
random number generator at the top of BATCH.TXT for each run!  Run GP-PARI
giving it at least 40 MB memory by using the '-s 40M' option.  At the GP
prompt, type '\r batch.txt'.  (If you are using the Windows precompiled
binary, you may need to give the path to batch.txt, as in '\r
/pari/batch.txt', and you may need to similarly change the \r lines in
batch.txt and rslll.txt.)  This will load the other two files and start a
search with a basis of 50 terms in random order, outputting any normal
solutions found to outsol and any prime solutions found to poutsol.  If
any of these solutions better the records at http://euler.free.fr, mail
them in!  If you want to change the basis used, modify the files as
necessary.  Note that changing the number of terms used by the LLL
algorithm requires changing BOTH the value of workterms AND the "1..50" in
the first vecextract command of each function.

Greg