function P = LagrangeRossz1 (X, Y)
#  X = [-1 0 1 2 3]
#  Y = [8 6 7 -4 0]
  P = zeros(1,length(X));
    
  for i=1:length(X)   
      szamlalo = poly(X([1:end]));
      nevezo = polyval(szamlalo, X);
      Li = szamlalo ./ nevezo;
      P = P+Li;
  end
  
#  ans = [1.75000  -6.00000  -0.25000   5.50000   6.00000]
endfunction