function P = Neville(X,Y,x)

Q=zeros(10);
Q(:,1)=Y';
for i= 2:10
    for j=2:i
        Q(i,j)=( (x-X(i-j+1))*Q(i,j-1)-(x-X(i))*Q(i-1,j-1) )/(X(i)-X(i-j+1));
    end
end
P=Q(10,10);