function LagrangianBasis(x,xcoor)
  
n  =  length(xcoor);
L  =  ones(n,1);


for i=1:n
  for j=1:n-1
      if abs(i-j)>0
        L(i)  =  L(i)*(x-xcoor(j))/(xcoor(i)-xcoor(j));
      end
  end              
end  
  
end
