clear
close all
clc

x=[0:20];
y=sin(x)+sin(2*x);
plot(x,y,'ro');
hold on

%Interpolacja funkcj¹ liniow¹
xi=[0:0.2:20];
yiL=interp1(x,y,xi,'linear');

% Interpolacja funkcj¹ sklejan¹
yiS=interp1(x,y,xi,'spline');

plot(xi,yiL);
pause;
plot(xi,yiS,'black--');