From: Fabien L. <laf...@gm...> - 2012-07-10 10:28:05
|
Hello everyone, I try to plot the digamma function of (1/2 + 1/x) but I'm not sure that I'm plotting the good one. I've tried: special.polygamma(0, (1/2 + 1/x)) and special.polygamma(1, (1/2 + 1/x)) but I don't have the same result as with mathcad. I've tried to code it like that: def F(x): return mpmath.diff(lambda x: gamma(1/2 + 1/x),1)/gamma(1/2 + 1/x) But It returns zero division error even when x is in ]0,1] Any idea? |