/* * Calculate Basel sum \sum_{i=1}^{\infty} \frac{1}{i*i} */ #include #include #define N 10000 int main(void) { long i; double s, ans; ans = M_PI*M_PI/6; s = 0.; for (i = 1; i <= N; i++) { s += 1./((double) i*i); } printf("%20.14f %20.14f %20.14f\n", s, ans, ans-s); return(0); }