Для каждого x, изменяющегося от a до b с шагом h, найти.

[Mat Bourn]

Всем доброго. Проверьте пожалуйста задачку. Математику я не знаю. Делал по аналогии чужих примеров.

Для каждого x, изменяющегося от a до b с шагом h, найти значения функции Y(x), суммы S(x) и |Y(x)–S(x)| и вывести в виде таблицы. Значения a, b, h и n вводятся с клавиатуры.
Работу программы проверить для a = 0,1; b = 1,0; h = 0,1; значение параметра n выбрать в зависимости от задания.

#include <stdio.h>
#include <iostream.h>
#include <windows.h>


 int _tmain(int argc, _TCHAR* argv[])
{
 SetConsoleOutputCP(1251);
 SetConsoleCP(1251);

 const double pi = 3.14;
 double a, b, h, n, x;
 double Yx, Sx, k;


 cout << "  a  =  ";
 cin >> a;

 cout << "  b  =  ";
 cin >> b;

 cout << "  h  =  ";
 cin >> h;

 cout << "  n  =  ";
 cin >> n;

 for ( x = a; x <= b; x += h )
{
 Yx = -0.5 * log(1 - 2 * x * cos(pi / 4) + pow(x,2));


 for( k = 1; k <= n; k++ )
{
 Sx = (pow(x,k) * cos(k * pi / 3)) / k;
}

 //cout << " |   " << Yx << "   |     " << Sx << "     |       " << Yx - Sx << "        |";
 printf (" |     %.6f",Yx);
 printf ("     |     %.6f",Sx);
 printf ("     |          %.6f",Yx - Sx);
 printf ("         |");
 cout << endl;
}

 cout << endl;
 cout << "  ";
 system("pause");
 return 0;
}

Использовал C++ Builder 10. Нахождение суммы S(x) я думаю что не правильно.

[ViktorR]

Есть:

Yx = -0.5 * log(1 - 2 * x * cos(pi / 4) + pow(x,2));

По условию задачи:

Yx = -0.5 * log(1 - 2 * x * cos(pi / 3) + pow(x,2));

А тут следует добавить:

Sx = 0;   // Тут накапливаем сумму
for( k = 1; k <= n; k++ )
{
 Sx =  Sx + pow(x,k) * cos(k * pi / 3) / k;
}

PS: Это так, взгляд с боку ...

[Mat Bourn]

Спасибо.