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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
%求已知数据点的第一类三次样条差值多项式及其插值点处的值
" O( \4 D$ L; S( g%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
%左端点的一阶导数：y_1
%右端点的一阶导数：y_n
5 O- ]5 O9 c* [: M% s% X%插值点的x坐标：x0
4 ?0 s; V9 }& e: g5 f' _%求得的三次样条差值多项式：f
%求得的x0处的插值：f0
syms t
f = 0.0000;
' e+ y: Y! `' yf0 = 0.0000;
if(length(x)==length(y))
    n = length(x);
else
    disp('x和y的维数不相等！');
( ~! q/ \4 G8 [+ K/ w' X( u    return;
end              %维数检查
  T" h  z4 l+ i: B, _for i=1:n
    if(x(i)<=x0)&&(x(i+1)>=x0)
        ihdex = i;
7 K8 L; W# v, ]) i3 k. C        break;
    end
4 P+ H2 e: u$ o. ]) C- N" ~end              %找到x0所在区间
9 L& _% z; m% X- @: k0 `/ r! c$ O4 dA = diag(2*ones(1,n));   %求解m的系数矩阵
u = zeros(n-2,1);
/ I- d! L, G) g% j& m' w: d  v* llamda = zeros(n-1,1);
c = zeros(n,1);
for i = 2:n-1
9 b' E2 ~1 F- h$ v$ e& E    u(i-1) = (x(i)-x(x-1))/(x(i+1)-x(i-1));
. y5 R' I& N" ~    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);
A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
end
c(1) = 2*y_1;
# G% l" G& a2 v+ A$ a5 x, E6 Q  sc(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+...
    y(index+1)*(2*(x(index+1)-t)+h)*(t-x(index))^2/h/h/h+...
  t# j# _; ~8 A. Y6 u# e4 R! Q    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所在区间的插值函数
* Z6 f" S* N! i' @/ M, |f0 = abs(f,'t',x0)  %x0处的插值
   
' \3 j' O: K0 ]2 ^, p2 w   
# i4 v8 o& @7 B2 u; X  T" d, E0 a