#include <stdio.h>
#include <math.h>
double S(double x, double eps)
{
double sum = x, term = x;
for(int n
= 0;fabs(term
) > eps
;++n
) sum += term *= -x*x/(2*n+2)/(2*n+3);
return sum;
}
int main()
{
const double eps = 1e-8;
printf(" x S(x) sin(x)\n---------------------------\n"); for(double x = 0.0; x < 1.6; x += 0.1)
printf("%3.1lf %.7lf %.7lf\n",x
,S
(x
,eps
),sin(x
)); }
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxtYXRoLmg+Cgpkb3VibGUgUyhkb3VibGUgeCwgZG91YmxlIGVwcykKewogICAgZG91YmxlIHN1bSA9IHgsIHRlcm0gPSB4OwogICAgZm9yKGludCBuID0gMDtmYWJzKHRlcm0pID4gZXBzOysrbikKICAgICAgICBzdW0gKz0gdGVybSAqPSAteCp4LygyKm4rMikvKDIqbiszKTsKICAgIHJldHVybiBzdW07Cn0KCmludCBtYWluKCkKewogICAgY29uc3QgZG91YmxlIGVwcyA9IDFlLTg7CiAgICBwcmludGYoIiB4ICAgICAgIFMoeCkgICAgICAgc2luKHgpXG4tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS1cbiIpOwogICAgZm9yKGRvdWJsZSB4ID0gMC4wOyB4IDwgMS42OyB4ICs9IDAuMSkKICAgICAgICBwcmludGYoIiUzLjFsZiAgICUuN2xmICAgJS43bGZcbiIseCxTKHgsZXBzKSxzaW4oeCkpOwp9Cg==