This version has been tested to give essentially the java results up to log t = -5, but diverges after that.
Authored by: guidrymw 2013-02-22
Parent: [r120]
Child: [r122]