#include <stdio.h>
#include <math.h>
int main()
{
double a, b, h,n;
puts("Input a, b,h ,n"); scanf("%lf%lf%lf%lf", &a
, &b
, &h
, &n
);
for(double x = a; x <= b; x += h)
{
double cs
[8] = { 1, 1/sqrt(2), 0, -1/sqrt(2) }; double S = 1, t = 1;
for(int k = 1; k <= n; ++k) S += cs[k%4]*(1 - 2*((k/4)%2))*(t *= x/k);
printf("%.3lf %15.7lf %15.7lf\n",x
,Y
,S
); }
}
I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxtYXRoLmg+CgppbnQgbWFpbigpCnsKCiAgICBkb3VibGUgYSwgYiwgaCxuOwogICAgcHV0cygiSW5wdXQgYSwgYixoICxuIik7IHNjYW5mKCIlbGYlbGYlbGYlbGYiLCAmYSwgJmIsICZoLCAmbik7CgogICAgZm9yKGRvdWJsZSB4ID0gYTsgeCA8PSBiOyB4ICs9IGgpCiAgICB7CiAgICAgICAgZG91YmxlIHhxID0geC9zcXJ0KDIpOwogICAgICAgIGRvdWJsZSBZID0gZXhwKHhxKSpjb3MoeHEpOwoKICAgICAgICBkb3VibGUgY3NbOF0gPSB7IDEsIDEvc3FydCgyKSwgMCwgLTEvc3FydCgyKSB9OwogICAgICAgIGRvdWJsZSBTID0gMSwgdCA9IDE7CiAgICAgICAgZm9yKGludCBrID0gMTsgayA8PSBuOyArK2spIFMgKz0gY3NbayU0XSooMSAtIDIqKChrLzQpJTIpKSoodCAqPSB4L2spOwoKICAgICAgICBwcmludGYoIiUuM2xmICAlMTUuN2xmICAlMTUuN2xmXG4iLHgsWSxTKTsKICAgIH0KfQoK