% Archivo Script para Splines cubicos con condiciones

% Datos
x=[];
y=[];
dy1=;
dyn=;
% Construimos la matriz de coeficientes
n=length(x)
h=diff(x);
M=2*diag([h 0]+[0 h])+diag(h,1)+diag(h,-1);
% Construimos los terminos independientes
q=diff(y)./h;
p=3*diff([dy1 q dyn]);
% Resolvemos el sistema para obtener los c's
c=p/M;
% Obtenemos b's y d's
cl=c(1:n-1);
cr=c(2:n);
b=q-h.*(2*cl+cr)/3;
d=(cr-cl)./(3*h);
% Coeficientes del Spline
a=y(1:n-1);
s=[d' cl' b' a'];
% Calculo del Spline
ps=mkpp(x,s);
% Evaluacion y representacion
xx=x(1)::x(n);
yy=ppval(ps,xx);
plot(x,y,'*',xx,yy)