%% function [A]=s102021105(fx,xi,p) % teacher revised
p = 0; % teacher revised
fx=[0.3 0.4 0.5 0.6];xi=[0.5488-0.3 .4493-.4 .3679-.5 .3012-.6];
  %run s102021105(fx,xi,0);then A(4,4)is the ans
  format long;fx=fx';xi=xi';n=size(fx);
  A(:,1) = fx;
  for i =2:n
    for j = i:n
      A(j,i) = (p-xi(j-i+1)*A(j,i-1)-(p-xi(j))*A(j-1,i-1))/(xi(j)-xi(j-i+1));
%     endfor
%     endfor  % teacher revised
    end
  end