我試圖解決在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)