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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
%求已知数据点的第一类三次样条差值多项式及其插值点处的值
+ s! E  u/ v4 `%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
2 A; _( |* @- \5 k%左端点的一阶导数：y_1
%右端点的一阶导数：y_n
- y7 Z1 u5 ^$ A0 n%插值点的x坐标：x0
* f" G2 _' s  ~0 d3 B# G%求得的三次样条差值多项式：f
( O3 @  N7 J7 p7 G% |1 O%求得的x0处的插值：f0
syms t
8 j/ ]# K, z' l6 I* Sf = 0.0000;
* H3 [8 j' b- a- ^% _f0 = 0.0000;
! `, |" _6 {8 t+ Y3 P4 m3 \if(length(x)==length(y))
    n = length(x);
else
    disp('x和y的维数不相等！');
    return;
  ]* I' V7 a" ^* D7 }, Pend              %维数检查
for i=1:n
    if(x(i)<=x0)&&(x(i+1)>=x0)
        ihdex = i;
% [: d0 X1 {- V        break;
    end
/ S+ |& k+ o0 Dend              %找到x0所在区间
A = diag(2*ones(1,n));   %求解m的系数矩阵
u = zeros(n-2,1);
lamda = zeros(n-1,1);
c = zeros(n,1);
for i = 2:n-1
$ ~; w1 S$ @' ]* }    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));
    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));
: A+ O5 \/ ~% ?7 `$ F. N+ }1 b A(i,i+1) = u(i-1);
A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
end
# I  V$ X  ]' ]! [. Y% @& e5 ?0 h# f& v: Pc(1) = 2*y_1;
c(n) = 2*y_n;
7 R2 p1 R% ~& A% J9 Q, Jm = followup(A,c);   %用追赶法求解方程组
* e2 n- c$ A# t2 R- V( o, d  Yh = x(index+1) - x(index);  %x0所在区间长度
f = y(index)*(2*(t-x(index))+h)*(t-x(index+1))^2/h/h/h+...
& d% o& Z$ G1 C) x/ J    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处的插值
# |- U4 k7 u9 K  p7 x" t8 P   
   

% \! m0 D  g( ~; x/ P( U  ?