Postup konstrukce s využitím Maplu

Jsou dána čísla x₁ ... x_k  a  f₁ ... f_k .

> k:=5:
    body:=array(1..k,[[x[1],f[1]],[x[2],f[2]],[x[3],f[3]],[x[4],f[4]],[x[5],f[5]]]);

Hledáme polynom L(x) takový, pro který platí

 L(x_i) = f_i

Lagrangerův polynom L(x) bude ve tvaru

> a:=array(1..k):
    for i from 1 by 1 to k do
        a[i]:= product( (x-body[j][1])/(body[i][1]-body[j][1]), j=1..i-1 )*
        product( (x-body[j][1])/(body[i][1]-body[j][1]), j=i+1..k )
    od;

Pro dané k bude mít výsledný polynom tvar

> L(x):=sum(body[n][2]*a[n], n=1..k);

procedura Lagrange()