--- mathmodule.c.original	2008-07-24 16:37:31.000000000 +0300
+++ mathmodule.c	2008-07-24 22:51:29.000000000 +0300
@@ -328,4 +328,8 @@
 FUNC1(cosh, cosh, 1,
       "cosh(x)\n\nReturn the hyperbolic cosine of x.")
+FUNC1(erf, erf, 0,
+      "erf(x)\n\nReturn the error function of x.")
+FUNC1(erfc, erfc, 0,
+      "erfc(x)\n\nReturn the complementary error of x.")
 FUNC1(exp, exp, 1,
       "exp(x)\n\nReturn e raised to the power of x.")
@@ -354,4 +358,6 @@
 	     "This is the largest integral value <= x.");
 
+FUNC1(lgamma, lgamma, 0,
+      "lgamma(x)\n\nReturn the logarithm of gamma function of x.")
 FUNC1(log1p, log1p, 1,
       "log1p(x)\n\nReturn the natural logarithm of 1+x (base e).\n\
@@ -367,4 +373,6 @@
 FUNC1(tanh, tanh, 0,
       "tanh(x)\n\nReturn the hyperbolic tangent of x.")
+FUNC1(tgamma, tgamma, 0,
+      "tgamma(x)\n\nReturn the gamma function of x.")
 
 /* Precision summation function as msum() by Raymond Hettinger in
@@ -1073,4 +1081,6 @@
 	{"cosh",	math_cosh,	METH_O,		math_cosh_doc},
 	{"degrees",	math_degrees,	METH_O,		math_degrees_doc},
+	{"erf",	math_erf,	METH_O,		math_erf_doc},
+	{"erfc",	math_erfc,	METH_O,		math_erfc_doc},
 	{"exp",		math_exp,	METH_O,		math_exp_doc},
 	{"fabs",	math_fabs,	METH_O,		math_fabs_doc},
@@ -1083,4 +1093,5 @@
 	{"isnan",	math_isnan,	METH_O,		math_isnan_doc},
 	{"ldexp",	math_ldexp,	METH_VARARGS,	math_ldexp_doc},
+	{"lgamma",	math_lgamma,	METH_O,		math_lgamma_doc},
 	{"log",		math_log,	METH_VARARGS,	math_log_doc},
 	{"log1p",	math_log1p,	METH_O,		math_log1p_doc},
@@ -1095,4 +1106,5 @@
 	{"tan",		math_tan,	METH_O,		math_tan_doc},
 	{"tanh",	math_tanh,	METH_O,		math_tanh_doc},
+	{"tgamma",	math_tgamma,	METH_O,		math_tgamma_doc},
  	{"trunc",	math_trunc,	METH_O,		math_trunc_doc},
 	{NULL,		NULL}		/* sentinel */