#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
double s(double x, double eps = 1e-10)
{
double a = x*x*x;
x = x*x;
double sum = a/3.0;
for(int k = 2; fabs(a) >= 4*eps*k*k; ++k)
{
sum += (a*=-x)/(4.0*k*k-1);
}
return sum;
}
int main(int argc, const char * argv[])
{
for(double x = 0.1; x <= 1.0; x += 0.1)
{
cout << setw(4) << x << setw(12)
<< ((1+x*x)*atan(x)-x)/2.0
<< setw(12) << s(x) << endl;
}
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8aW9tYW5pcD4KI2luY2x1ZGUgPGNtYXRoPgoKdXNpbmcgbmFtZXNwYWNlIHN0ZDsKCmRvdWJsZSBzKGRvdWJsZSB4LCBkb3VibGUgZXBzID0gMWUtMTApCnsKICAgIGRvdWJsZSBhID0geCp4Kng7CiAgICB4ID0geCp4OwogICAgZG91YmxlIHN1bSA9IGEvMy4wOwogICAgZm9yKGludCBrID0gMjsgZmFicyhhKSA+PSA0KmVwcyprKms7ICsraykKICAgIHsKICAgICAgICBzdW0gKz0gKGEqPS14KS8oNC4wKmsqay0xKTsKICAgIH0KICAgIHJldHVybiBzdW07Cn0KCgoKaW50IG1haW4oaW50IGFyZ2MsIGNvbnN0IGNoYXIgKiBhcmd2W10pCnsKICAgIGZvcihkb3VibGUgeCA9IDAuMTsgeCA8PSAxLjA7IHggKz0gMC4xKQogICAgewogICAgICAgIGNvdXQgPDwgc2V0dyg0KSA8PCB4IDw8IHNldHcoMTIpCiAgICAgICAgICAgIDw8ICgoMSt4KngpKmF0YW4oeCkteCkvMi4wCiAgICAgICAgICAgIDw8IHNldHcoMTIpIDw8IHMoeCkgPDwgZW5kbDsKICAgIH0KfQoKCg==