% Session p. 124/125

f = inline('ReturnProbability(epsilon,n)-0.5','epsilon','n');
for n=10*2.^(0:3)
  out = sprintf('n = %3i\t N = %6i\t\t epsilon* = %18.16f',n,...
     (2*n+1)^2-1,fzero(f,[0.06,0.07],optimset('TolX',1e-16),n)); 
  disp(out);
end