SoloCodigo

/*Figueroa Bernal jaime eduardo
UNAM FES ACATLAN METODOS NUMERICOS II*/
#include<stdio.h>
#include<conio.h>
#include<stdlib.h>
#include<ctype.h>
double *cons(int);
void leer(double []int);
void dest(double []);
void solution(double[]double[]doubleint);
void main(void)
{
   int ing;
   char dt;
   double *X*Yein;
   clrscr();
   textcolor(10);
   cprintf("\t\t\t Interpolacion de Lagrange.\n");
   puts("\n¨No de elementos de la tabla?:");
   scanf("%d"&n);
   X=cons(n);
   Y=cons(n);
   do
   {
   puts("\n¨Grado del polinomio.?");
   scanf("\t%d"&g);
   if(g<n)
   {
      puts("¨Punto a interpolar.?");
      scanf("%lf"&in);
      do
      {
         puts("\nIngresa los elementos de x:");
         leer(Xn);
         puts("Ingresa los elemtos de F(x).");
         leer(Yn);
         puts("¨Datos correctos.?");
         puts("(s/n)");
         }while(d=tolower(getch())!='s');
         solution(XYing);
         dest(X);
         dest(Y);
   }
    else
   {
      puts("No hay sufucientes puntos.");
   }
   getch();
   puts("\n¨Deseas continuar.?\n");
   puts("(s/n)");
   }while(t=tolower(getch())!='n');
 }
double *cons(int N)
{
   double *T;
   T=(double*)calloc(Nsizeof(double));
   if(T==NULL)
   {
      puts("No hay memoria disponible");
   }
   return(T);

}
void leer(double A[]int N)
{
   int i;
   for(i=0;i<N;i++)
   {
      printf("\t%d="i);
      scanf("%lf"&A);
   }
}
void dest(double A[])
{
   free(A);
}
void solution(double X[]double Y[]double xint gr)
{
   int ij;
   double num=1den=1s=0;
   for(i=0;i<=gr;i++)
   {
      for(j=0;j<=gr;j++)
      {
         if(i!=j)
         {
            num*=(x-X[j]);
            den*=(X-X[j]);
         }
      }

      num*=Y;
      printf("\nL(%d)=%lf"inum/den);
      s+=num/den;
      num=den=1;
   }
   printf("\nEl resultado es=%lf"s);
}











Autor: jaime