程序有点问题？？？？？？？帮帮忙

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
7 G$ e9 J, D. T4 _  a: c& i8 J%求已知数据点的第一类三次样条差值多项式及其插值点处的值
%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
8 K  v* @# B: [%左端点的一阶导数：y_1
%右端点的一阶导数：y_n
! m; x( Q$ w5 U7 |%插值点的x坐标：x0
$ w/ Z. N1 u$ e$ U0 A9 `%求得的三次样条差值多项式：f
%求得的x0处的插值：f0
8 [1 x' Z" I) [+ e8 Psyms t
f = 0.0000;
& n: R; Z8 T8 G9 e$ O, if0 = 0.0000;
3 }* ^9 s) c6 K* y" B8 z7 Rif(length(x)==length(y))
    n = length(x);
$ x! M, Z) i/ p; ^6 [else
1 Z- s) Z! A& C* H& A    disp('x和y的维数不相等！');
    return;
. E1 u2 f1 n: N; P( e. o% V( b' @/ iend              %维数检查
for i=1:n
    if(x(i)<=x0)&&(x(i+1)>=x0)
) y) ^: \/ K; P        ihdex = i;
        break;
    end
9 g& B; P. B9 y: j3 Oend              %找到x0所在区间
A = diag(2*ones(1,n));   %求解m的系数矩阵
( t+ ^, N. I" `2 i/ \- tu = zeros(n-2,1);
lamda = zeros(n-1,1);
c = zeros(n,1);
+ a0 N+ ?: F1 f/ d$ Y  \# P6 ufor i = 2:n-1
) c, P% \7 a8 q    u(i-1) = (x(i)-x(x-1))/(x(i+1)-x(i-1));
    lamda(i) = (x(i+1)-x(i))/(x(i+1)-x(i-1));
3 L& P5 X. ~, J/ c: ]! A1 V    c(i) = 3*lamda(i)*(y(i)-y(i-1))/(x(i)-x(i-1))+3*u(i-1)*(y(i+1)-y(i))/(x(i+1)-x(i));
2 e; g8 V# t* }% i, O; m A(i,i+1) = u(i-1);
A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
end
c(1) = 2*y_1;
c(n) = 2*y_n;
m = followup(A,c);   %用追赶法求解方程组
h = x(index+1) - x(index);  %x0所在区间长度
f = y(index)*(2*(t-x(index))+h)*(t-x(index+1))^2/h/h/h+...
5 V& A2 J# e3 a/ B) w    y(index+1)*(2*(x(index+1)-t)+h)*(t-x(index))^2/h/h/h+...
    m(index)*(t-x(index))*(x(index+1)-t)^2/h/h-...
    m(index+1)*(x(index+1)-t)*(t-x(index))^2/h/h; %x0所在区间的插值函数
f0 = abs(f,'t',x0)  %x0处的插值
   
   
) k1 s' w+ i4 y' n* |
( i# t! m$ J! s/ p! i- j