# HG changeset patch
# User Burcin Erocal <burcin@erocal.org>
# Date 1337093678 -7200
# Node ID 5637fa8f6b36c2c33b786b2a70493db7dbddd557
# Parent  7d14009cc1693eaeabeb27ca4e955a512fe8b8e6
Add doctests to check Pynac evaluates psi1 numerically if input is not exact.

diff --git a/sage/functions/other.py b/sage/functions/other.py
--- a/sage/functions/other.py
+++ b/sage/functions/other.py
@@ -1002,6 +1002,11 @@
             \psi\left(x\right)
             sage: loads(dumps(psi1(x)+1))
             psi(x) + 1
+
+            sage: t = psi1(x); t
+            psi(x)
+            sage: t.subs(x=.2)
+            -5.28903989659219
         """
         GinacFunction.__init__(self, "psi", nargs=1, latex_name='\psi',
                 conversions=dict(maxima='psi[0]', mathematica='PolyGamma'))