前回のは[-π/2,π/2]で成り立っていたが、 実数全域でも成り立つようにしよう。 素直に書くと、 double sin2(double x) { int i = (int)(x * (1 / M_PI) + 0.5); if(x i--; if(i & 1) x = i * M_PI - x; else x = x - i * M_PI; double x2 = x * x; return x …