Configuracion Spline

%%SPLINE%%
X=[0.0 4.5 6.8 9.5 11.1 12.4 13.2 13.8 15.0 15.7 16.2 17.5 18.8 19.5 20.8 22.5 25.2 28.5 32.6 35.0];
Y=[0.00.6 1.4 3.1 5.1 8.1 10.2 12.1 15.5 17.3 18.8 20.2 20.8 21.0 21.1 21.2 21.3 21.5 21.9 22.1];Figure,plot(X,Y,’g’),title(‘tarea spline’),xlabel(‘x’),ylabel(‘F(x)’);
L=length(X); a=Y’;
%%calculo de h y diferencias de a%%
For i=1:L-1;h(i,1)=X(1,i+1)-X(1,i);
da(i,1)=Y(1,i+1)-Y(1,i);
end
%%matriz de diagonal 1%%
m=zeros(L,L);
n=seros(L,1);
for i=1:L;
m(i,i)=1;end
for i=1:L-2;
sh(i,1)=2*(h(i+1,1)+h(i,1));
sa(i,1)=(-3/h(i,1))*da(i,1)+(3/h(i+1,1))*da(i+1,1);
end
for i=2:L-1;n(i,1)=sa(i-1,1);
m((i,i)=sh(i.1,1);
end
j=0;
for i=1:2:L;
m(2,i)=h(i-j,1);
j=1;
end
c=inv(m)*n;
d=zeros(L,1);
fori=1:L-1;
d(i,1)=(c(i+1,1)-c(i,1))/(3*h(i,1));
end
b=zeros(L,1);
for i=1:L-1;b(i,1)=(1/h(i,1))*da(i,1)-(h(i,1)/3)*(2*c(i,1)+c(i+1,1));
end
disp(‘m yn:’);[a,b,c,d]
disp(‘el valor que desea aproximar:’);
X0=input(‘X0’);
IfX0>=X(1,1)&X0<=X(1,L);
for i=1:L-1;
if X0>=X(1,i)&X0<X(1,i+1);
k=I;
end
end
elseif X0>X(1,L);
disp(‘elvalor esta fuera del intervalo’)
end
vy=a(k,1)+b(k,1)*(X0-X(1,k))+c(k,1)*(X0-X(1,k))˄2+d(k,1)*(X0-X(1,k))˄3