Question

Неточно выводит значения. Не смогла найти, где ошиблась , и отмечу как лучший ответ кину спам, если не в тему напишите

using namespace std;
#include
#include
#include
#include
HANDLE c;
int main() {
c = GetStdHandle(STD_OUTPUT_HANDLE);
SetConsoleTextAttribute(c, 11);
setlocale(LC_ALL, "ru");
double x, y, s, a, b, h, next;
int k;

cout << " a= ";
cin >> a;
cout << endl << " b= ";
cin >> b;
cout << endl << " h= ";
cin >> h;
cout << endl << " k= ";
cin >> k;
cout << "\nЗначения x, y, s для каждого x+const*h из [a; b].\n\n";

x = a;
b = b + (h / 3);
do {
y = cosh(x);
s = 0;
next = 1;
for (double w = 0; w < k; w++) {
next = next * ((x * x) / (4 * w * w + 6 * w + 2));
s = s + next;
}
cout << setw(5) << x << endl;
cout << setw(14) << y << endl;
cout << setw(14) << setprecision(3) << s << endl;
cout << "\n\n";
x = x + h;
} while (x < b);
return 0;
}

Задание
Вывести на экран таблицу значений функции y(x) и ее разложения в ряд s(x) для x, изменяющихся от a до b с шагом h=(b-a)/10.
a=0.1
b=1
y(x)=ch(x)
s(x)=sum(от n=0 до k) (x^(2n))/((2n)!)