/* error5.c
積み残し現象
1/(1^2) + 1/(2^2) + 1/(3^2) + ... の理論値は (pi^2)/6 だが ...
*/

#include<stdio.h>
#define pi 3.14159265358979324 
main()
{
    int i,j,n;
    float x;
    printf("理論値は %f\n\n",pi*pi/6);
    printf("昇順に加えると\n",n);
    for(j=2; j<=7; j++){
        n=(int)pow(10,j);
        x=0.0;
        for(i=1;i<=n;i++) x+=1.0/i/i;
        printf("%8d項目までの和 = %f\n",n,x);
    }
    printf("\n");
    printf("降順に加えると\n",n);
    for(j=2; j<=7; j++){
        n=(int)pow(10,j);
        x=0.0;
        for(i=n;i>=1;i--) x+=1.0/i/i;
        printf("%8d項目までの和 = %f\n",n,x);
    }
}