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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
5 }- ?$ R) l) B! O, r( j- B%求已知数据点的第一类三次样条差值多项式及其插值点处的值
, p" S0 Z3 T- S4 A+ n%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
%左端点的一阶导数：y_1
! s$ V; M: f' G6 d- i%右端点的一阶导数：y_n
. [& h) p' f* o6 \%插值点的x坐标：x0
* d8 l, t& O6 t+ X' J- y%求得的三次样条差值多项式：f
%求得的x0处的插值：f0
syms t
f = 0.0000;
f0 = 0.0000;
if(length(x)==length(y))
. A# P- u6 J* \+ v    n = length(x);
8 F2 M5 ^) v5 [: Felse
    disp('x和y的维数不相等！');
6 G+ J* O: m  i7 _    return;
end              %维数检查
for i=1:n
    if(x(i)<=x0)&&(x(i+1)>=x0)
        ihdex = i;
        break;
$ i+ Q; U8 p6 B- E9 z3 j    end
end              %找到x0所在区间
: B# k9 W' Z/ l- XA = diag(2*ones(1,n));   %求解m的系数矩阵
" E# l3 L3 x, Ru = zeros(n-2,1);
lamda = zeros(n-1,1);
% }/ `/ ], R- \2 V; c* W6 j+ vc = zeros(n,1);
1 J* X6 L+ M2 |5 o- \% ?9 {# Zfor i = 2:n-1
1 q5 q" L5 L! V3 b    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));
" ?+ w$ d% ^! J+ {; \ A(i,i+1) = u(i-1);
A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
* h. ^2 a/ C, X" r: ^; e. I1 aend
c(1) = 2*y_1;
c(n) = 2*y_n;
$ o- U$ q2 X8 c5 Qm = 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+...
: m/ I5 p( ^  Y3 }" f    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所在区间的插值函数
% O4 M$ m9 n1 h8 K- @; ^f0 = abs(f,'t',x0)  %x0处的插值
6 W0 j" i( ?' Y6 B. B   
/ y$ }- ^4 x; Z3 w  L$ Z0 J! b