--- test_math.py.original	2008-07-24 16:38:11.000000000 +0300
+++ test_math.py	2008-07-24 21:57:19.000000000 +0300
@@ -261,4 +261,18 @@
         self.ftest('degrees(-pi/4)', math.degrees(-math.pi/4), -45.0)
 
+    def testErf(self):
+        self.assertRaises(TypeError, math.erf)
+        self.ftest('erf(0)', math.erf(0), 0.0)
+        self.assertEquals(math.erf(INF), 1.0)
+        self.assertEquals(math.erf(NINF), -1.0)
+        self.assert_(math.isnan(math.erf(NAN)))
+
+    def testErfc(self):
+        self.assertRaises(TypeError, math.erfc)
+        self.ftest('erfc(0)', math.erfc(0), 1.0)
+        self.assertEquals(math.erfc(INF), 0.0)
+        self.assertEquals(math.erfc(NINF), 2.0)
+        self.assert_(math.isnan(math.erfc(NAN)))
+
     def testExp(self):
         self.assertRaises(TypeError, math.exp)
@@ -403,4 +417,25 @@
             self.assert_(math.isnan(math.ldexp(NAN, n)))
 
+    def testLgamma(self):
+        self.assertRaises(TypeError, math.lgamma)
+        self.ftest('lgamma(1.0)', math.lgamma(1.0), 0.0)
+        self.ftest('lgamma(2.0)', math.lgamma(2.0), 0.0)
+        self.ftest('lgamma(3)', math.lgamma(3), math.log(2))
+        self.assertEquals(math.lgamma(INF), INF)
+        self.assertEquals(math.lgamma(NINF), INF)
+        self.assertRaises(ValueError, math.lgamma, -1)
+        self.assertRaises(ValueError, math.lgamma, 0)
+        self.assert_(math.isnan(math.lgamma(NAN)))
+
+    def testTgamma(self):
+        self.assertRaises(TypeError, math.tgamma)
+        self.ftest('tgamma(0.5)', math.tgamma(0.5), math.sqrt(math.pi))
+        self.ftest('tgamma(-0.5)', math.tgamma(-0.5), -2*math.sqrt(math.pi))
+        self.ftest('tgamma(3.5)', math.tgamma(3.5), math.sqrt(math.pi)*15/8)
+        self.ftest('tgamma(-1.5)', math.tgamma(-1.5), math.sqrt(math.pi)*4/3)
+        self.assertEquals(math.tgamma(INF), INF)
+        self.assertRaises(ValueError, math.tgamma, 0)
+        self.assert_(math.isnan(math.tgamma(NAN)))
+
     def testLog(self):
         self.assertRaises(TypeError, math.log)