% interpolacija funkcije f(x) = sin(2*pi*x/5)
% metode: linearna interpolacija - default ili dodati "linear"
%         kubicna interpolacija - "cubic"
%         kubicna spline interpolacija - "spline"
% koristimo f-ju interp1() sa pripadnim argumentima
% sve rezultate crtamo na istom grafu radi lakse usporedbe
% Numericka matematika

xf = [0:0.05:10];  % tocke za crtanje grafa - os x
yf = sin(2*pi*xf/5); %tocke za crtanje grafa - os y

xp = [0:10]; %tocke interpolacije - os x
yp = sin (2*pi*xp/5); %tocke interpolacije - os y

%interpolacija pomocu funkcije interp1() u tockama (xp, yp)
lin = interp1 (xp, yp, xf);
cub = interp1 (xp, yp, xf, "cubic");
spl = interp1 (xp, yp, xf, "spline");

%graficki prikaz rezultata
plot (xf,yf,"r", xf,lin,"b", xf,cub,"c", xf,spl,"m", xp,yp,"r*"); legend ("original", "linear", "cubic", "spline");