pi2.cc~

CFD First time code using gnuplot C++ and STL

//Inclue bibliotecas de C
#include <stdlib.h>
#include <stdio.h>
//Fun玢o main
int main(int argc , char * argv [])
{
register double width , sum;
register int intervals , i; //N鷐ero de intervalos
intervals = atoi(argv [1]);
width = 1.0 / intervals ; //largura
sum = 0;
for (i=0; i<intervals ; ++i)
{
register double x = (i + 0.5) * width;
sum += 4.0 / (1.0 + x * x);
}
sum *= width;
//Mostra o resultado
printf ( "Estimation of pi is %f\n" , sum);

return(0);
}