% ZAD 2
clear;
close all
clc

a=2;
b=12;
n=3;
x = linspace(a,b,n+1);
xbis = linspace(a,b,500);
y=sin(x)+cos(x);
ybis=sin(xbis)+cos(xbis);
y11= interp1(x,y,xbis,'spline');

% hold on
% title('DRUGI lin')
plot(x,y,'o',xbis,y11);
blad_s11 = max(abs(ybis-y11))


y12= interp1(x,y,xbis,'linear');

hold on
plot(x,y,'o',xbis,y12, 'b');
blad_s12 = max(abs(ybis-y12))
d1=ybis-y12;
d2=ybis-y11;
subplot(211), plot(xbis, d1)
subplot(212), plot(xbis, d2)