format long

x=[0.6 0.7 0.8 0.9];
fx=[0.6-cos(0.6) 0.7-cos(0.7) 0.8-cos(0.8) 0.9-cos(0.9)];
n=4;
 for k=1:n
Q(1:n,1)=fx;
for i=2:n
for j=2:i
Q(i,j)=((-x(i-j+1))*Q(i,j-1)-(-x(i))*Q(i-1,j-1))/(x(i)-x(i-j+1));
end
end
end
Q(5,5)