% Change format to show many decimals
format long

f=@(x)sum(1./gamma(1:x));

f(100)