#include "stdafx.h"
#include <iostream>
#include <math.h>
using namespace std;
double y(double x);
int main(int argc, char* argv[])
{
double a = 0.1, b = 1.0, h = 0.1;
double sum, prev;
int k;
cout << "k = "; cin >> k;
for (double x = a; x < b + h / 2; x += h) {
//start values
sum = prev = pow(x,3) / 3;
for (int n = 2; n <= k; ++n) {
prev *= -1 * (x * x);
sum += prev / (4 * n * n - 1);
}
cout << "y(x) = " << y(x) << " AND s(x) = " << sum << endl;
}
return 0;
}
double y(double x)
{
return (1 + x * x) / 2 * atan(x) - x / 2;
}
I2luY2x1ZGUgInN0ZGFmeC5oIgojaW5jbHVkZSA8aW9zdHJlYW0+CiNpbmNsdWRlIDxtYXRoLmg+CiAKdXNpbmcgbmFtZXNwYWNlIHN0ZDsKCmRvdWJsZSB5KGRvdWJsZSB4KTsKIAppbnQgbWFpbihpbnQgYXJnYywgY2hhciogYXJndltdKQp7Cglkb3VibGUgYSA9IDAuMSwgYiA9IDEuMCwgaCA9IDAuMTsKCWRvdWJsZSBzdW0sIHByZXY7CglpbnQgazsKCgljb3V0IDw8ICJrID0gIjsgY2luID4+IGs7CiAKCWZvciAoZG91YmxlIHggPSBhOyB4IDwgYiArIGggLyAyOyB4ICs9IGgpIHsKCQkvL3N0YXJ0IHZhbHVlcwoJCXN1bSA9IHByZXYgPSBwb3coeCwzKSAvIDM7CgkJZm9yIChpbnQgbiA9IDI7IG4gPD0gazsgKytuKSB7CgkJCXByZXYgKj0gLTEgKiAoeCAqIHgpOwoJCQlzdW0gKz0gcHJldiAvICg0ICogbiAqIG4gLSAxKTsKCQl9CiAKCQljb3V0IDw8ICJ5KHgpID0gIiA8PCB5KHgpIDw8ICIgICBBTkQgcyh4KSA9ICIgPDwgc3VtIDw8IGVuZGw7Cgl9CiAKCXJldHVybiAwOwp9Cgpkb3VibGUgeShkb3VibGUgeCkgCnsKCXJldHVybiAoMSArIHggKiB4KSAvIDIgKiBhdGFuKHgpIC0geCAvIDI7Cn0=