三次样条插值的Matlab

function [m,H,lambda,mu,D,A,dY,s]=ClampedSp_1 (X,Y,dyo,dyn)
m=length(X);A=zeros(m,m); n=m-1;H=zeros(1,n);lambda=zeros(1,n);
mu=zeros(1,n);lambda(1)=1;A(1,1)=2;A(1,2)=lambda(1);
D=zeros(1,n);mu(n)=1;
for k=1:n
    hk=X(k+1)-X(k);H(k)=hk;
end
for k=1:n-1
    lambdak=H(k+1)/(H(k)+H(k+1)); lambda(k+1)=lambdak;
    muk