// Mat01169, script para Metodo de Hutton para Pi 
function d=digse(ra,r)
 if (r==ra) | (abs(r)<%eps)
    d=16.0;   // maximo que pode ser obtido em IEEE 754 //
   // com precisao dupla (por enquanto...) //
 else
    d=max(0,-log10( 2*abs( (ra-r)/r ))); // digse //
 end
endfunction; // funcao que calcula os digses //

// Calculo de Pi pelo Metodo de Hutton //
nmax=20;
sn=4; n=0;  // inicializacao de s0 //
printf("%i %20.16f %10.4f \n",0,sn,0); 
for n=1:nmax
  sa=sn;
  sn=sn + 4*(-1)**n*( (1/2)**(2*n+1) + (1/3)**(2*n+1) )/(2*n+1);  
  // atualizacao do termo geral //
  printf("%i %20.16f %10.4f \n",n,sn,digse(sa,sn)); 
end
// criado Ago 2011 //

