% Archivo Script para Splines lineales

% Datos
x=-1:0.5:1; % Modificar en cada caso
y=1./(1+25*x.^2); % Modificar en cada caso

% Calculamos el spline
n=length(x); 
a=y(1:n-1);
b=diff(y)./diff(x);
s=[b',a']; % Coeficientes de los q_k
ps=mkpp(x,s); % Codificacion del Spline

% Dibujamos el Spline
xx=x(1):.1:x(n); % Abscisas para dibujar el spline(Modificar en cada caso)
yy=ppval(ps,xx); % Evaluamos el Spline
plot(x,y,'*',xx,yy)


hold on
yy=1./(1+25*xx.^2);
plot(xx,yy)