Polinomal İnterpolasyon

Numerik analiz dersinde kökleri girdi alarak fonksiyonu, polinomal interpolasyon yöntemi ile bir program yapalım.
#include "stdio.h"
#include<math.h>
#include<stdlib.h>
int main()
{
printf("Bu proje Polinomal Interpolasyon yöntemiyle girilen kökleri ven" 
" o kökteki değerleri girdi şeklinde alarak denklemi ortaya cikariyornnn"); 

int i,j,k,bilinmeyen_sayisi,denklem_sayisi
,a,b;
float katsayi,c,dizi[99][99]={0};
char cevap,e,E;
int nokta;

printf("Lütfen Nokta Sayisini Giriniz:");
scanf("%d",&nokta);
int x[nokta],y[nokta];
printf("nLutfen Noktalarin Sirasiyla Degerlerini Giriniz:n");
for(i=0;i<nokta;i++){
printf("X%d=",i+1);
scanf("%d",&x[i]);
}
printf("nLutfen Noktadaki Degerleri Giriniz:n");
for(i=0;i<nokta;i++){
printf("Y%d=",i+1);
scanf("%d",&y[i]);
} 

denklem_sayisi=nokta;
bilinmeyen_sayisi=nokta;

/*burda girdiğimiz degerleri matrise yerleştiriyoruz*/
for(i=0;i<denklem_sayisi;i++)
{
for(j=0;j<=bilinmeyen_sayisi;j++ )
{
if(j!=bilinmeyen_sayisi)
{dizi[i][j]=pow(x[i],j);}
if(j==bilinmeyen_sayisi)
{dizi[i][j]=y[i];}

}
}



/*burada satırlardaki işlemleri yok ediyoruz. Matristeki - ile carpıp eklemeyi düşünün*/

for(i=0;i<denklem_sayisi;i++)
{
for(j=0;j<denklem_sayisi;j++)
{
if(i!=j)
{
katsayi=(-1)*(dizi[j][i]/dizi[i][i]);
for(k=0;k<=bilinmeyen_sayisi;k++ )
{
dizi[j][k]=dizi[j][k]+katsayi*dizi[i][k];
}
}
}
}

/*değerler son olarak ekrana basılıyor?*/
printf("nbu fonksiyonun katsayilarin degeri = n");
for(i=0;i<denklem_sayisi;i++)
{
dizi[i][i]=dizi[i][bilinmeyen_sayisi]/dizi[i][i];
printf("ndenklemdeki a(%d) nin katsayisi : %fn",i,dizi[i][i]);
}
system ("pause");
return 0;

}