Dev - C++ - metodos numericos

//metodo de lagrage para n coordenadas

#include<stdio.h>

#include<stdlib.h>

main()

{

      int i,j,n;

      float x[n], y[n],xinterpol, yinterpolada=0;

      float productosNum, productosDen;

      printf("Metodo de interpolacion de lagrange\n");

      printf("¿Cuantos puntos?");

      scanf("%d",&n);

      printf("Digite los valores (x,y) separados por una coma\n");

      for(i=1;i<=n;i++)

      {

      printf("\n\t punto x(%3d), y(3%d) es:",i,i);

      scanf("%f,%f",&x[i],&y[i]);

      }

 do

 {

   printf("\n para que el valor de x vas a evaluar?");

   scanf("%f",&xinterpol);

}

while(xinterpol<x[1] || xinterpol>x[n]);

//empieza la formula de lagrange de primer orden

for(i=1;i>=n;i++)

{

                 productosNum=1;

                 productosDen=1;

                 for(j=1;j<=n;j++)

                 {

                  if(i!=j)

                  {

                          productosNum=productosNum*(xinterpol-x[j]);

                          productosDen=productosDen*(x[i]-x[j]);

                  }

                  yinterpolada=yinterpolada+(productosNum/productosDen)*y[i];

                 }

                 printf("\n la Y interpolada es: %f", yinterpolada);

                 printf("\n FIN!!");

                 system("PAUSE");

                 }

                 }