Помогите разобраться с программой С++ Builder

[spectrum988]

Для каждого x, изменяющегося от a до b с шагом h, найти значения функции Y(x), суммы S(x) и |Y(x)–S(x)| и вывести в виде таблицы. Значения a, b, h и n вводятся с клавиатуры. Так как значение S(x) является рядом разложения функции Y(x), при правильном решении значения S и Y для заданного аргумента x (для тестовых значений исходных данных) должны совпадать в целой части и в первых двух-четырех позициях после десятичной точки.
Работу программы проверить для a = 0,1; b = 1,0; h = 0,1; значение параметра n=7.
Помогите пожалуйсто,не знаю как сделать,чтобы программа работала правильно.
Условие.PNG
#include <stdio.h>
#include <conio.h>
#include <math.h>
void main()
{
double a, b, x, h,s,y,r,rez;
const double PI = 3.1415926;
int n, k;
puts("Input a,b,h,n");
scanf("%lf%lf%lf%d", &a, &b, &h, &n);
printf("\n x | Y(x) | S(x) | |Y(x)-S(x)|");
printf("\n------------------------------------------");
for(x = a; x<=b; x+=h)
{
r=s=1;
y=-1/2*log(1-2*x*cos(PI/3)+pow(x,2));
for(k=1;k<=n;k++)
{
s=(pow(x,k)*cos(k*PI/3))/k;
}
rez=fabs(y-s);
printf("\n %1.1lf | %1.4lf | %1.4lf | %1.4lf |", x,y,s,rez);
}
printf("\n--------------------------------------------");
system "pause";
}

[Smitt&Wesson]

s=(pow(x,k)*cos(k*(PI/3)))/k;
Вроде так.

[spectrum988]

Ну да,согласен,но почему-то в самой программе не считает значение y(выводит на экран 000000000000)а функцию s cчитает нормально.

[BDA]

y=-log(1-2*x*cos(PI/3)+pow(x,2))/2;

1/2 = 0 (целая арифметика)

[spectrum988]

Ответы выражений у и s должны быть одинаковыми,но почему-то они не совпадают.

[BDA]

#include <stdio.h>
#include <conio.h>
#include <math.h>

int
main()
{
    double a, b, x, h, s, y, rez;
    //const double PI = 3.1415926;
    int n, k;
    puts("Input a, b, h, n");
    scanf("%lf%lf%lf%d", &a, &b, &h, &n);
    printf("\n x | Y(x) | S(x) | |Y(x)-S(x)|");
    printf("\n------------------------------------------");
    for(x = a; x <= b; x += h) {
        s = 0;
        y = -log(1 - 2 * x * cos(M_PI / 3) + x * x) / 2;
        for(k = 1; k <= n; ++k)
            s += (pow(x, k) * cos(k * M_PI / 3)) / k;
        rez = fabs(y - s);
        printf("\n %1.1lf | %1.4lf | %1.4lf | %1.4lf |", x, y, s, rez);
    }
    printf("\n--------------------------------------------");
    return 0;
}

[spectrum988]

Cпасибо за помощь,программа работает правильно,буду признателен если обьясните мою ошибку,что я не учел.

[BDA]

Вы не считали сумму (s += ...), а переписывали значение в s новым числом (s = ...). После цикла в s хранился последний вычисленный элемент.

[spectrum988]

Спасибо еще раз,очень выручили,скажите а запись M_PI это предопределенная константа,которую не нужно присваивать в самой программе(вы написали мне как в комментарии //const double PI = 3.1415926; ).Значит будет и так и так правильно ?!

[BDA]

Пожалуйста.
Да, предопределенная, поэтому можно самому не объявлять её, но если откажется компилироваться, то просто расскоментируете свое определение.
Точнее, в math.h есть такая запись: #define M_PI 3.14159265358979323846.