(1,N);
x(1) = sqrt(2);
for i = 1:(N-1)
  x(i+1) = sqrt(2 + x(i));
  wynik = 2 * 2^i/prod (x(1:i));
  disp ([i, wynik, pi - wynik]);
endfor

