function [koef,m]=oscilirajociPolinom(F,a,b,n)
format long
vektor=a:((b-a)/(n+1)):b;

y=F(vektor);
y=y(:);
A=ones(n+2,n+2);
for i = 1:(n+2)
    for j=2:(n+1)
    A(i,j)=vektor(i)^(j-1);
    end
end
for i=1:(n+2)
    A(i,n+2)=(-1)^(i-1);
end
koef=A\y;
m=koef(n+2);
koef=koef(1:end-1);


end