function y=base_lagrange(x)
    n=length(x)
    for i=1:n
        a=x
        a(i)=[]
        y(1)=poly(a,'x')
        y(1)=y(i)/horner(y(i),x(i))
    end
end function