我试图解决在Matlab中的函数：digamma（X）= C

Question

例如，C=0，

solve('psi(x)=0') 

ans = 
-226.83295306016122662496413158295 

psi(ans) 

???Error using ==> psi 
Input must be single or double. 

我不能得到正确的答案

Answer 1

有趣......这似乎是一个在solve给我的错误...无论我试图解决什么值，我总是得到一个奇怪的值大约-227。甚至当我尝试通过给digamma的近似欺骗MATLAB，我得到了相同的结果或者更糟：

>> solve('(gamma(x+0.01)-gamma(x))/0.01/gamma(x)=0') 
ans = 
    matrix([[-226.83790783643886637282996154237]]) 

>> solve('(gammaln(x+0.01)-gammaln(x-0.01))/0.02 = 0') 
??? Error using ==> mupadmex 
Error in MuPAD command: cannot differentiate equation [numeric::fsolve] 

以下数值方法的工作原理：

%// value of the digamma to solve for 
Y = -10; 

%// Solve using numerical scheme 
X = fsolve(@(x)psi(max(0,x)) - Y, exp(Y)) 

%// Check solution: psi(X) ≈ Y 
psi(X)