Неверное решение, ошибка при компиляции.

[Janliv]

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

#include <stdio.h>						
#include <conio.h>
#include <cmath>
int main()

{

double a,b,h,x,s,r,y,z;
int n, zn = -1, k;
puts("Input a,b,h,n");
scanf("%lf %lf %lf %d", &a, &b, &h, &n);
for (x=0, x<b, x+=h)
	{
		y=((1+x*x)/2)*(atan(x))-x/2;
		r=s=1;
			for (k=1, k<=n, k++)
				{	
				r=pow(zn,k+1)*((pow(x,2*k+1)/(4*k*k-1)));
				s+=r;
				}
	printf ("\n Y(x)= %lf \n",y );
	printf ("\n S(x)= %lf \n",s );
	z=fabs(y-s);
	printf ("\n |Y(x)-S(x)|=%lf \n",z);
	} 
getch();
return 0;
}