x=[0 1 -1 3];
f=[1 2 2 0];
newtonuv_polynom(1,x,f);