MATLABで簡易逆ガンマ関数を実装する

ガンマ関数を使えば階乗を一般化でき便利である。

\begin{align}
\Gamma(n+1)=n!
\end{align}

今回は\(n!\)となる\(x\)の値をmatlabを使って求める。

n=8256;

f=@(x) abs(n-gamma(x));
x0=3;

x=fminsearch(f,x0);

gamma(x)

プログラムではfminsearchを使っているので、いずれ使用しない形に書き改める。