%uloha 122

k=0:10; x=k-5; y=1./(1+x.^2);   %tecky nutne
p=polyfit(x,y,length(x)-1)  %interpolacni polynom
x_jemne = [0:.25:10]-5;
y1=polyval(p,x_jemne);   %hodnoty interpolacniho polynomu
y2=spline(x,y,x_jemne);  %hodnoty splajnu
y_fce=1./(1+x_jemne.^2); %presne funkcni hodnoty

plot(x_jemne,y_fce,'r',x,y,'r+',x_jemne,y1,'m',x_jemne,y2)
legend('funkce','uzlove body','interp.polynom','splajn')

x2=[.5:4.5]'; porovani=[x2 1./(1+x2.^2) polyval(p,x2), spline(x,y,x2)]