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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
%求已知数据点的第一类三次样条差值多项式及其插值点处的值
* a* A# y( j- I2 h6 B3 D%已知数据点的x坐标向量：x
5 q( @7 O- A% z%已知数据点的y坐标向量：y
%左端点的一阶导数：y_1
* w& @# s1 Z' R; C, c%右端点的一阶导数：y_n
%插值点的x坐标：x0
%求得的三次样条差值多项式：f
0 v& I9 I, _: L: G%求得的x0处的插值：f0
syms t
. H9 C! ]' b3 Z# d6 o9 xf = 0.0000;
f0 = 0.0000;
( F" H0 U' w: xif(length(x)==length(y))
    n = length(x);
else
    disp('x和y的维数不相等！');
    return;
end              %维数检查
2 b: g# Y+ l6 e; i$ T" ^% F! ufor i=1:n
8 y. ^. \6 h1 ^0 @" b    if(x(i)<=x0)&&(x(i+1)>=x0)
6 q* N& ]0 B' q, c' T4 ?6 @        ihdex = i;
, a: d$ z" y: }/ ]+ F& N        break;
    end
8 i; _; Y1 G) g3 lend              %找到x0所在区间
A = diag(2*ones(1,n));   %求解m的系数矩阵
u = zeros(n-2,1);
lamda = zeros(n-1,1);
c = zeros(n,1);
2 B; ~2 I; ~% c+ p' e% Hfor i = 2:n-1
    u(i-1) = (x(i)-x(x-1))/(x(i+1)-x(i-1));
: g/ r0 N! A% w4 [    lamda(i) = (x(i+1)-x(i))/(x(i+1)-x(i-1));
$ K. q7 k- b2 w3 [7 \, I    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(i,i+1) = u(i-1);
A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
end
  Q- h+ L4 q' i& W$ x+ P! rc(1) = 2*y_1;
c(n) = 2*y_n;
& p6 D7 D" P( }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+...
    y(index+1)*(2*(x(index+1)-t)+h)*(t-x(index))^2/h/h/h+...
2 ?: B& T3 x( r$ U7 v    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处的插值
' A: C* y3 |0 V4 a; t. P5 `, Q, k