function coeff = Divdiffcoeff(X,Y)

F = zeros(10);
F(:,1)=Y';
for i=2:10
    for j=2:i
        F(i,j)=(F(i,j-1)-F(i-1,j-1))/(X(i)-X(i-j+1));
    end
end
coeff=diag(F);