#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==