#include <stdio.h>
#include <stdlib.h>
typedef double T;
T f( T x )
{
return 1/(x*x);
}
T F( T x )
{
return -1/x;
}
T Integral0( T left, T right, long N, T (*pf)(T) )
{
T S = 0.0;
T x, dx;
dx = (right - left) / N;
for ( x = left; x <= right; x += dx ) {
S += dx * ( pf( x ) + pf( x+dx ) ) / 2.0;
}
return S;
}
int main()
{
T L = 0.000001;
T R = 1000;
T F0 = F(R)-F(L);
long n;
for ( n=10000; n<100000000; n*=1.2 ) {
printf( "%ld %g\n",n, (Integral0(L,R,n,f)-F0)/F0);
}
return 0;
}