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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
) P1 r& i" R% |) K5 s6 z%求已知数据点的第一类三次样条差值多项式及其插值点处的值
6 W" D' a9 I9 R3 F5 C%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
1 S1 m! e; d6 V4 K4 P2 P: f%左端点的一阶导数：y_1
%右端点的一阶导数：y_n
%插值点的x坐标：x0
%求得的三次样条差值多项式：f
: H1 ?/ K- \# a# a6 M9 `%求得的x0处的插值：f0
syms t
, S) u& W! x8 A0 f8 o) Sf = 0.0000;
f0 = 0.0000;
if(length(x)==length(y))
; h4 n% c- p0 y1 _6 M+ s; ?( ]    n = length(x);
! Z& N) s7 E/ ^0 V4 G+ xelse
3 j  S- {2 K" `" Q2 a6 i0 J    disp('x和y的维数不相等！');
    return;
7 {2 |, r% F/ @1 W1 z; fend              %维数检查
for i=1:n
2 t; j- F! t' i    if(x(i)<=x0)&&(x(i+1)>=x0)
. {  p" k! u; L5 P# h% {        ihdex = i;
6 X% h' H+ f  _. D        break;
    end
end              %找到x0所在区间
A = diag(2*ones(1,n));   %求解m的系数矩阵
, ~* k8 K4 p+ g; d. pu = zeros(n-2,1);
lamda = zeros(n-1,1);
c = zeros(n,1);
for i = 2:n-1
    u(i-1) = (x(i)-x(x-1))/(x(i+1)-x(i-1));
  F3 D/ E, |9 V1 r+ G    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(i,i+1) = u(i-1);
5 z3 f( b  C* L5 u4 L A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
end
c(1) = 2*y_1;
c(n) = 2*y_n;
m = followup(A,c);   %用追赶法求解方程组
3 E. {8 K  k' r- fh = 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+...
1 a7 @' t1 b/ J) \# Z    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处的插值
" H  q2 d: U4 t" I   
; F7 z3 _- A' t& c