По численным методам

[VovaZZ]

в общем проблема в написании интерполяции ньютона назад
вот мой код:
#include<iostream.h>
#include<fstream.h>
double Nx(double,double,int);
double d(int,int);
double *a;
double d(int n,int k)
{
double f;
if(n==1) f=a[k+1]-a[k];
else f=d(n-1,k+1)-d(n-1,k);
return f;
}

double Nx(double f0,double t,int n)
{
double s=f0,tk=t,fakt=1;
for(int k=0;k<n-1;k++)
{
s+=tk*d(k+1,0)/fakt;
fakt*=(k+2);
tk*=(t-k-1);
}
return s;
}

void main()
{ifstream in ("dannye.txt");

double x0,h;
in>>x0>>h;
a=new double [5];
for(int k=0;k<5;k++)
in>>a[k];
double x;
cin>>x;
double t=(x-x0)/h;
cout<<'\n'<<Nx(a[0],t,5);
delete []a;
}

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