#include #include double const pi = 3.1415926535897932385; static double f0(double x){ if (x<0) return f0(-x); if (x>pi) return -f0(x-pi); return 1-2*x/pi;} static double fs(double x, int n){ double s = 0; while (n--) if(n&1) s += cos(n*x)/(n*n); return 8./(pi*pi)*s;} static void rp(double x){ printf("%e %20.15f %20.15f\n", x, f0(x), fs(x,10000));} int main(){ rp(pi/2); rp(0); rp(1); rp(-4); rp(7); rp(11.53); return 0;}