function A=St(x,f) [baris,kolom]=size(x); A=zeros(baris,baris); A(:,1)=f; for j = 1:baris-1 for i = 1:baris - j A(i,j+1)=(A(i+1,j)-A(i,j))/(x(i+j)-x(i)); end end endfunction