<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Mar 6, 2018 at 4:38 AM, Matej Usaj <span dir="ltr"><<a href="mailto:m.usaj@utoronto.ca" target="_blank">m.usaj@utoronto.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi All,<br>
<br>
I'm trying to solve function that contains a call to psi (digamma) but<br>
my call to fsolve fails every time, ie:<br>
<br>
> f=psi(x/2)==0;<br>
psi(1/2*x)==0<br>
> fsolve(f,x,2.7,3.1);<br>
fsolve(): function does not evaluate numerically<br>
<br>
but I'd expect the result to be x~=2.923...<br>
<br>
Am I missing something here?<br></blockquote><div><br></div><div>Can it be that psi function returns values in terms of the Euler constant? Then you may need to make float evaluation wrapper for it.</div><div><br></div><div>Best wishes,</div><div>Vladimir</div><div><div>-- </div><div>Vladimir V. Kisil                 <a href="http://www.maths.leeds.ac.uk/~kisilv/">http://www.maths.leeds.ac.uk/~kisilv/</a></div><div>  Book:     Geometry of Mobius Transformations     <a href="http://goo.gl/EaG2Vu">http://goo.gl/EaG2Vu</a></div><div>  Software: Geometry of cycles          <a href="http://moebinv.sourceforge.net/">http://moebinv.sourceforge.net/</a></div></div><div><br></div></div></div></div>