三次样条插值函数matlab程序

function y = nbcspline（x0 y0 x）
n = length（x0）;   m = length（x）;
y = zeros（1 m）;
p（n） = 0;  q（1） = 0;  d（1） = 0;  d（n） = 0;
for k = 2:n-1
    h（1） = x0（k） - x0（k-1）;  h（2） = x0（k+1） - x0（k）;
    p（k） = h（1） / （h（1） + h（2））;
    q（k） = h（2） / （h（1） + h（2））;
    d（k） = 6*（（y0（k+1） - y0（k））/h（2） - （y0（k） - y0（k-1）） / h（1）） / （h（1） + h （2））;
end
b（1） = q（1） / 2; 
for k = 2:n-1
    b（k） = q（k） / （2 - p（k）*b（k-1））;
end

temp（1） = d（1） / 2;
for k = 2:n
    temp（k） = （d（k） - p（k）*temp（k-1）） /