Lagrange's Interpolation

#include <stdio.h>

int main()
{
	int i, j, n;

	printf("Enter the number of data points: ");
	scanf("%d", &n);

	float x[n], y[n];

	printf("Enter the values of x and y:\n");
	for (i = 1; i <= n; i++)
	{
		scanf("%f %f", &x[i], &y[i]);
	}
	
	int xp, yp = 0;
	printf("Enter interpolation point: ");
	scanf("%f", &xp);
	printf("\n");
	
	for (i = 1; i <= n; i++)
	{
		int p = 1;
		for (j = 1; j <= n; j++)
		{
			if (i != j)
			{
				p = p * (xp - x[j]) / (x[i] - x[j]);
			}
		}
		yp = yp + p * y[i];
		printf("%.3f\n", yp);
	}

	printf("Interpolated value at %.3f is %.3f\n", xp, yp);

	return 0;
}