655a656,686> static PyObject *> math_expm1(PyObject *self, PyObject *arg)> {> double x = PyFloat_AsDouble(arg);> if (PyErr_Occurred())> return NULL;> double k = 2.;> double s = x;> double r = x;> double u0;> double u1;> if (fabs(x) < log(2.)) {> do {> u0 = x/k*s;> u1 = x/(k+1)*u0;> s = u0;> r += u0;> k++;> } while (fabs(u0)>fabs(u1));> }> else {> r = exp(x) -1;> }> return PyFloat_FromDouble(r);> > }> > PyDoc_STRVAR(math_expm1_doc,> "expm1(x)\n\n\> Return exp(x) minus 1 for small x");> 1366a1398> {"expm1", math_expm1, METH_O, math_expm1_doc}, 1408a1441> PyModule_AddObject(m, "C", PyFloat_FromDouble(Py_MATH_C));