function P =NevilleDeriv(X,Y,x)

Q=zeros(10);
Qprime=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));
        Qprime(i,j)= ( Q(i,j-1) + (x-X(i-j+1))*Qprime(i,j-1) - Q(i-1,j-1)  - (x-X(i))*Qprime(i-1,j-1) )/(X(i)-X(i-j+1));
    end
end
P=Qprime(10,10);