function [C,L]=lagran(X,Y)
%variables de entrada 
%          - X es un vector que contiene las 
%            abcisas de los puntos
%          - Y es un vector que contiene las 
%            ordenadas de los puntos 
%variables de salida 
%          - C es una matriz cuyas entradas son 
%            los coeficientes del polinomio 
%            de interpolación de Lagrange (VECTOR COLUMNA)
%          - L es una matriz cuyas entradas son 
%            los coeficientes del polinomio 
%            de interpolación de Lagrange

w=length(X);
n=w-1; 
L=zeros(w,w); 
%Formación de los coeficientes del polinomio

for k=1:n+1
   V=1; 
   for j=1:n+1 
      if k~=j
         V=conv(V,poly(X(j)))/(X(k)-X(j)); 
      end 
   end 
   L(:,k)=V; 
end
 
%Cálculo de los coeficientes del polinomio 
%interpolador de Lagrange 

L
C=L*Y'