% Cvorovi interpolacije
x = [1 2 3 4]; y = [0 1 2 0];

% Konstruisemo splajnove
spl_per = csape(x, y, 'periodic');          %periodicni
spl_nat = csape(x, y, 'second');            %prirodni
spl_es  = csape(x, [2 y 2], 'complete');    %zadati izvodi

% Vrednost periodicnog splajna u 2.5
ppval(spl_per, 2.5)

% Crtanje splajnova na [1, 4]
xx = linspace(1, 4);
plot(xx, ppval(spl_per, xx), ...
     xx, ppval(spl_nat, xx), ...
     xx, ppval(spl_es, xx), ...
     x, y, 'o')
legend('periodicni', 'prirodni', 'zadati izvodi', 'cvorovi')