function [Ln] = LagrangeInterpolacio(x,y)
  %{
    -- [Ln] = LagrangeInterpolacio(x,y)
    Az adott alappontokhoz és függvényértékekhez tartozó Lagrange-polinom
    kiszámítása (a bázispolinomok összegeként).
    
    Bemenõ paraméterek
    x:      alappontok
    y:      függvényértékek
    
    Kimenõ paraméter
    Ln:     a Lagrange polinom együtthatóvektora
  %}
  
  %%% Input ellenõrzése:
  x  = x(:);
  y  = y(:);
  n  = length(x);
  ny = length(y);
  if n != ny || !isfloat(x) || !isfloat(y)
    warning('LagrangeInterpolacio: hibás dimenziójú/típusú input adatok');
    Ln = Nan(n)
    return
  endif  
  
  %%% A tényleges számolás:
  LBazis = LagrangeBazisPolinomok(x);
  Ln = zeros(1,n);
  for i=1:n
   Ln = 
  endfor
  
endfunction
 
% A Lagrange polinomot a bázispolinomok segítségével írjuk föl;
% a LagrangeBazisPolinomok(x) segédfüggvény által visszaadott
% LBazis nxn-es mátrix i-dik sorában az i-dik Li(x) bázispolinom
% együtthatói állnak

function [LBazis] = LagrangeBazisPolinomok(x)
    n = length(x);
    LBazis = zeros(n,n);
    for i = 1:n
      Li = poly(x([1:i-1,i+1:n]));
      LBazis(i,1:n) = 
    endfor  
endfunction