clear all; close all

m = 16;   h=1/m;   x=0:h:1; x=x';   y=cos(pi*x);

plot(x,y,'o')
hold
n = 2; 

p= my_ls(x,y,m,n),

yl = zeros(m+1,1);

for i=1:m+1
  yl (i) = 0;
  for j=1:n+1
     yl(i) = yl(i) + p(j) *x(i)^(j-1);
%  yl(i), j, yl(i),p(1)
  end
end

plot(x,yl)