［Octave］計算機イプシロンを求める

計算機イプシロンとは、不等式「1 + epsm > 1」が成り立つ最小の正数epsmのこと。これを手計算で求めてみる。

>> epsm = 1.0;
>> while (1.0 + epsm > 1.0)
> epsm /= 2.0;
> end
>> epsm *= 2.0;
>> disp(epsm)
2.2204e-16

この環境下では約2.2×10^-16であることがわかる。