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? Thanks for your help! Cheers, Matej