a = -5;
b = 5;
k = 0.5;
x1 = (a:k:b); %ekvidistantni cvorovi
n = numel(x1);
x2=[]; % ovde stavljamo cebisevljeve cvorove
y1=[]; % vektori za y vrednosti u ekvidistantnim tackama
y2=[]; % -||- u cebisevljevim tackama
for i=1:n
    x2(i) = (a+b)/2+(b-a)/2*cos((n-i+0.5)/n*pi);
    y1(i)=1/(1+25*x1(i)^2);
    y2(i)=1/(1+25*x2(i)^2);
end

p1 = polyfit(x1,y1,2); %polinomi drugog stepena
p2 = polyfit(x2,y2,2);
p1s = spline(x1,y1); %splajnovi
p2s = spline(x2,y2);

%iscrtavamo sve zajedno
plot(x1,polyval(p1,x1),x2,polyval(p2,x2),x1,ppval(p1s,x1),x2,ppval(p2s,x2));
hold on

%prava funkcija (na jako finom (sitnom) intervalu)
x = (a:0.01:b);
y = [];
for i=1:numel(x)
   y(i)=1/(1+25*x(i)^2); 
end
plot(x,y);

legend('Ekvidistantni pol2','Cebisevljevi pol2','Ekvidistantni splajn','Cebiseljevi splajn','Prava funkcija');