Thanks! The problem came from the 1/2 ! For convenience I've found the
function digamma on numpy
*http://docs.scipy.org/doc/scipy/reference/generated/scipy.special.psi.html
But it's quite hard to find it! Maybe we can ask to add digamma in the
title between parenthesis?
Fabien
*
2012/7/10 Damon
On Tue, Jul 10, 2012 at 08:57:24AM -0400, Benjamin Root wrote:
> On Tue, Jul 10, 2012 at 7:05 AM, Damon McDougall
> wrote:
>
> > On Tue, Jul 10, 2012 at 12:27:59PM +0200, Fabien Lafont wrote:
> >
> > > But It returns zero division error even when x is in ]0,1]
> >
> > I think it blows up at x = 0.
On Tue, Jul 10, 2012 at 7:05 AM, Damon McDougall
wrote:
> On Tue, Jul 10, 2012 at 12:27:59PM +0200, Fabien Lafont wrote:
> > 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.polygamm
On Tue, Jul 10, 2012 at 12:27:59PM +0200, Fabien Lafont wrote:
> 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))
You wan
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):
