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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
%求已知数据点的第一类三次样条差值多项式及其插值点处的值
%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
%左端点的一阶导数：y_1
5 o' ^! D% Y; t) q* A+ h, K%右端点的一阶导数：y_n
%插值点的x坐标：x0
%求得的三次样条差值多项式：f
( v3 @- O2 g0 i$ b$ o%求得的x0处的插值：f0
7 Y1 W3 o+ i! T% D, fsyms t
f = 0.0000;
1 N7 h" ]( u! T; @( @# @f0 = 0.0000;
, s, s, ?7 C2 Q1 _5 R6 Dif(length(x)==length(y))
8 Z: Z. i8 b& A3 w( q( Y    n = length(x);
else
    disp('x和y的维数不相等！');
) c. e/ r! h5 L1 H0 E9 u    return;
/ V/ P+ R9 `+ [) @' oend              %维数检查
for i=1:n
/ }6 m9 q2 O0 n' }* k% z. D    if(x(i)<=x0)&&(x(i+1)>=x0)
        ihdex = i;
4 v2 ^; B+ f3 p2 v1 K. Q2 _        break;
    end
( n+ S2 A5 o$ e/ Q; u! Z5 |end              %找到x0所在区间
A = diag(2*ones(1,n));   %求解m的系数矩阵
2 j* u. J) p7 W: v' e' P' P' Iu = zeros(n-2,1);
% X  k0 X  c( Q: t$ clamda = zeros(n-1,1);
c = zeros(n,1);
7 D7 s; s  i, ^+ 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));
, R8 C5 G" n' \0 A# `6 |" b. {$ x    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
7 ]/ ^" [5 {4 ?4 d7 ^c(1) = 2*y_1;
c(n) = 2*y_n;
m = followup(A,c);   %用追赶法求解方程组
h = x(index+1) - x(index);  %x0所在区间长度
" B( y$ h. }2 d3 ]5 T# v! k9 Uf = 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+...
    m(index)*(t-x(index))*(x(index+1)-t)^2/h/h-...
& s+ n3 z/ ~8 Q3 f4 L) q6 h    m(index+1)*(x(index+1)-t)*(t-x(index))^2/h/h; %x0所在区间的插值函数
f0 = abs(f,'t',x0)  %x0处的插值
# }5 D* E+ I/ S' l: `   
7 I  Q2 q$ r  r$ e' z   

, [6 }* H8 @5 y* Q! `7 t