Hi All,

I'm trying to solve function that contains a call to psi (digamma) but
my call to fsolve fails every time, ie:

> f=psi(x/2)==0;
psi(1/2*x)==0
> fsolve(f,x,2.7,3.1);
fsolve(): function does not evaluate numerically

but I'd expect the result to be x~=2.923...

Am I missing something here?