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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
%求已知数据点的第一类三次样条差值多项式及其插值点处的值
%已知数据点的x坐标向量：x
3 g" ~2 g* m# D7 `8 X$ f%已知数据点的y坐标向量：y
, z' r8 h: p4 g7 Y%左端点的一阶导数：y_1
%右端点的一阶导数：y_n
3 y& Q; {& i; v/ h4 x! Y: v' J  I9 y%插值点的x坐标：x0
* A' B: u6 L- ~9 x! x$ |%求得的三次样条差值多项式：f
' Y3 P# L2 S% A" c- B%求得的x0处的插值：f0
syms t
2 n* t* c6 `0 D! [$ @4 Mf = 0.0000;
f0 = 0.0000;
if(length(x)==length(y))
# F, \- Z/ b7 R  M    n = length(x);
% D2 d6 q5 A1 h+ E8 l. q  T0 ielse
7 s3 D: a4 L# u) f5 ]3 \1 }8 i    disp('x和y的维数不相等！');
, }& Z8 ]2 _4 a1 _  _  Q" G( Z) b    return;
end              %维数检查
for i=1:n
    if(x(i)<=x0)&&(x(i+1)>=x0)
        ihdex = i;
; T8 t' r. J3 b6 O6 R, y2 m        break;
9 k+ U! q( t6 J    end
end              %找到x0所在区间
A = diag(2*ones(1,n));   %求解m的系数矩阵
u = zeros(n-2,1);
$ F6 u3 P: u) s. o9 E( m$ r# [) E* u3 ?lamda = zeros(n-1,1);
% y: n  U/ D5 u) e4 m0 q2 O: b* ~c = zeros(n,1);
3 N: ?) [2 u: `  Ffor i = 2:n-1
    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(i,i+1) = u(i-1);
A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
4 u9 l5 o2 `4 t% \$ Q9 B  L" v9 Vend
c(1) = 2*y_1;
# Y' ~8 D7 J4 C- N' U5 Ec(n) = 2*y_n;
) |. z# J0 n% J4 um = followup(A,c);   %用追赶法求解方程组
" H) Y0 Q& L8 \5 Z1 n  \8 [* yh = x(index+1) - x(index);  %x0所在区间长度
f = y(index)*(2*(t-x(index))+h)*(t-x(index+1))^2/h/h/h+...
: \- h3 J6 `; |. \    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处的插值