% Test 1, calculo numerico (521230), S2-2011
% Tema 5: Interpolacion
% Pregunta 3

% construyendo spline cubico natural
n=10;
i=0:n;
x = -1+2*i/n;
fx = 1./(25*x.^2 + 1);
cs = csape(x,fx,[2 2]);

% puntos para evaluar spline y funcion
t = linspace(-1,1);
% evaluar f
ft = 1./(25*t.^2 + 1);
% evaluar spline
cst = fnval(cs,t);
% error de interpolacion
plot(t,abs(ft-cst));