#include <stdio.h>
#include <conio.h>
#define maxint 32767

float table[10][10];
int n;

void initable()
{
	int i,j;
	for(i=0;i<10;i++)
		for(j=0;j<10;j++) table[i][j]=maxint;
}

void getvals()
{
	int i;
	printf("Enter the number of terms: ");
	scanf("%d",&n);
	for(i=0;i<n;i++) scanf("%f %f",&table[i][0],&table[i][1]);
}

void getdiff(int d)
{
	int i;
	for(i=0;i<=n-d-1;i++)
	{
		table[i][d+1]=(table[i+1][d]-table[i][d])/(table[d+i][0]-table[i][0]);
	}
}

void showtable()
{
	int i,j;
	printf("\nThe table after finding differences is \n\n");
	for(i=0;i<n;i++)
	{
		for(j=0;j<n+1;j++)
		{
			if (table[i][j]!=maxint) printf("%10.2f",table[i][j]);
		}
		printf("\n");
	}
}

float interpolate(float x)
{
	float ans=table[0][1],xval=1;
	int i;
	for(i=0;i<n-1;i++)
	{
		xval=xval*(x-table[i][0]);
		ans=ans+table[0][i+2]*xval;
	}
	return ans;
}


void main()
{
	int j;
	float iplt=0,ans=0;
	clrscr();
	initable();
	getvals();
	for(j=1;j<n+1;j++) getdiff(j);
	showtable();
	printf("\nEnter the value of x to interpolate at: ");
	scanf("%f",&iplt);
	ans=interpolate(iplt);
	printf("\n\nThe value of f(%.2f) is %.2f",iplt,ans);
	getch();
}