#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
double series(double x, double eps)
{
double t = 1, s = 1;
for(int k = 1; abs(t) > eps; ++k)
s += t *= x*(3-2*k)/(2*k);
return s;
}
int main(int argc, char * argv[])
{
for(double x = 0; x < 1; x += 0.1)
{
cout << setprecision(3) << setw(5) << x
<< setprecision(8) << setw(10) << sqrt(1+x) << " "
<< setprecision(8) << setw(10) << series(x,1e-8) << endl;
}
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8aW9tYW5pcD4KI2luY2x1ZGUgPGNtYXRoPgoKdXNpbmcgbmFtZXNwYWNlIHN0ZDsKCmRvdWJsZSBzZXJpZXMoZG91YmxlIHgsIGRvdWJsZSBlcHMpCnsKICAgIGRvdWJsZSB0ID0gMSwgcyA9IDE7CiAgICBmb3IoaW50IGsgPSAxOyBhYnModCkgPiBlcHM7ICsraykKICAgICAgICBzICs9IHQgKj0geCooMy0yKmspLygyKmspOwogICAgcmV0dXJuIHM7Cn0KCmludCBtYWluKGludCBhcmdjLCBjaGFyICogYXJndltdKQp7CiAgICBmb3IoZG91YmxlIHggPSAwOyB4IDwgMTsgeCArPSAwLjEpCiAgICB7CiAgICAgICAgY291dCA8PCBzZXRwcmVjaXNpb24oMykgPDwgc2V0dyg1KSA8PCB4CiAgICAgICAgICAgIDw8IHNldHByZWNpc2lvbig4KSA8PCBzZXR3KDEwKSA8PCBzcXJ0KDEreCkgPDwgIiAgICIKICAgICAgICAgICAgPDwgc2V0cHJlY2lzaW9uKDgpIDw8IHNldHcoMTApIDw8IHNlcmllcyh4LDFlLTgpIDw8IGVuZGw7CgogICAgfQp9Cg==