массивы в C++

[Monkey]

Вроде бы все просто... Есть такие переменные:
double func(double x,double y)
{return 100.*(y-x*x)*(y-x*x)+(1-x)*(1-x);}
double f[7];
double h;
/**/
пытаюсь сделать с ними вот что:

for (k=0;k<n+1;k++){f[k]=func(x[k],y[k]);}
t=0;
for (k=0;k<n+1;k++){
if(k!=n){h=fabs(f[k]-f[k+1]);}
if(k==n){h=fabs(f[k]-f[0]);}
if(h<eps) t++;
cout<<"h="<<h<<endl;
fprintf(u,"\nh=%d",h) ;
}

а переменная h почему-то не хочет принимать значения разности f, получается вот что:
x[0]= 1.0000
y[0]= 2.0000
x[1]= 1.9659
y[1]= 2.2588
x[2]= 1.2588
y[2]= 2.9659
f[0]=100.0000
f[1]=258.8712
f[2]=190.8661
h=11276
h=-29110
h=-13874

Подскажите, с чем это может быть связано?