function alfa = InterpolaceMonomialy(f,x)

n = length(x)-1;
A = zeros(n+1,n+1);
b = zeros(n+1,1);
for i=0:n
    for j=0:n
        A(i+1,j+1) = x(i+1)^j;
    end
    b(i+1) = f(x(i+1));
end
cond(A)
alfa = A\b;