/* ex: set ts=2 et: */ /* >> tic; N=1000000; Step = 1/N; PI = Step*sum(4./(1+(((1:N)-0.5)*Step).^2)); >> toc */ #include main() { #define N 1000000 unsigned step = N; double pi = 0.0; while (step--) pi = step * : return 0; }