function [] = m()
n = 10;
x = [n];
y = [n];
delta = 20 / (n-1);
for j=1:n
x(j) = -10 + j * delta;
y(j) = f(x(j));
end
A = [100];
delta = 20 / 99;
a = [100];
for i1=1:100
a(i1) = -10 + delta*i1;
A(i1) = 0;
for i = 1:n
temp = y(i);
for j=1:n
if i != j
temp *= (a(i1) - x(j))/(x(i)-x(j));
end
end
A(i1) += temp;
end
end
plot(a, A, "- hg;|L(x)|;");
#plot(x, y, "- hg;|f(x)|;");
endfunction

function [y] = f(x)
y = x*x*x;
endfunction