#include
#include
using namespace std;
int main()
{
setlocale(LC_ALL, "Rus");
double a, b, s, h, y, t, q;
double x = a = -2, b = -0.1;
h = (b - a) / 10;
int n;
for (x = a; x < b; x += h)
{
y = x * cos(x);
s = x;
t = x;
q = -x * x;
for (n = 1; n < 141; n++)
{
t = t * q / (2 * n * (2 * n - 1));
s += t;
cout << "y= " << y << setw(12) << "s= " << s << endl;
}
}
}

Надеюсь, вопрос по этой задачи - последний, хотя кто его знает... А насчёт "обогрева вселенной" - это бы ещё чья корова мычала, ибо сумма по условию задачи должна браться от нулевого n до n, равного ста сорокà, хотя достаточно взять количество членов ряда на порядок меньше - результаты не изменятся. Вон же ссылка есть во втором ответе, кому надо - пусть проверят и сам во всём убедятся! Ой, умора, ха-ха-ха!.
(✷‿✷)