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?
------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and
threat landscape has changed and how IT managers can respond. Discussions
will include endpoint security, mobile security and the latest in malware
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-users
