Ряд разложения функции. В чем ошибся?

Question

[yaroslavklimuk]

Набросал, но не уверен что верно работает, может кто-нибудь помочь?

Для каждого 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 = 15;

for(x = a; x<=b; x+=h)
{
                p1 = s = p = 1;
                for(k = 1; k<=n; k++)
                {
                        f=pow(-1,k+1)*(pow(x,2*k))/(2*k*(2*k-1);
                        s+=f;

                  }

y=-log(sqrt(1+pow(x,2)))+x*atan(x);
p1=fabs(y-s);