trigamma for very large arguments (PR#14020)

1 message · Jerry.Lewis at BiogenIdec.com

Wed, Oct 21, 2009 7:20 PM #

Full_Name: Jerry W. Lewis
Version: 2.9.2
OS: Windows XP Professional
Submission from: (NULL) (96.237.55.233)


trigamma(x) returns 0 for x>1e152, yet
  trigamma <- function(x) 1/x
gives machine accuracy for any x>1e16