Интерполяционный многочлен Ньютона

#include <stdio.h>
#include <conio.h>
#include <math.h>
#define NUMB 10
#define X0  -1.5
#define Xn  -0.6
#define h   0.1
int fact(int n)
{
    if (n==1 || n==0) return 1;
   return n*fact(n-1);
}
double f(double x)
{
    return x*x-4*x*cos(x);
}
double ticalc(double t,int n)
{
    if (n==1) return 1;
    double mult;
   for (int k=0;k<n-1;k++)
    mult=t*ticalc(t-1.0,n-1);
   return mult;
}
double N(double x)
{
    int k,i=0;
    double dely0,t,Y[10],sum=0.0,p;
   t=(x-X0)/h;
   for (p=X0;p<=Xn;p=p+h);
   {
    Y[i]=f(p);
      i++;
   }
   dely0=Y[1]-Y[0];
   sum+=Y[0];
   for (k=1;k<=NUMB;k++)
   sum+=(ticalc(t,k+1)/fact(k))*pow(dely0,k);
   return sum;
}
void main()
{
    int i=1;
    double x;//=0.406,sum=2.6744;
   for (x=X0;i<=NUMB;x+=h)
    {
    printf ("V TOCHKE %0.2lf : ZNACHENIE FUNKCII = %0.4lf ZNACHENIE MNOGOCHLENA = %0.4lf\n",x,f(x),N(x));
      i++;
   }
getch();
}

dely0=Y[1]-Y[0];
...
...
sum+=(ticalc(t,k+1)/fact(k))*pow(dely0,k);

   for (int k=0;k<n-1;k++)
        mult=t*ticalc(t-1.0,n-1);