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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
%求已知数据点的第一类三次样条差值多项式及其插值点处的值
3 J. K: \2 X$ |5 X7 ~1 e/ i. K' e%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
6 L, k" {1 Y& e%左端点的一阶导数：y_1
%右端点的一阶导数：y_n
/ @4 Z& Q) t# v# i+ Q. E. b( d0 i%插值点的x坐标：x0
%求得的三次样条差值多项式：f
%求得的x0处的插值：f0
# N* P6 e) Y& M: jsyms t
4 o& M( G9 Y5 G4 `9 x4 Nf = 0.0000;
8 c9 P' w. a  u. wf0 = 0.0000;
6 \* q8 s. L! I$ i/ v. cif(length(x)==length(y))
    n = length(x);
8 Q. m0 e) {* }! Pelse
    disp('x和y的维数不相等！');
, w; H+ U2 R5 N$ y. \: q8 y! g, I    return;
end              %维数检查
5 Z: N7 x% K% w; A8 u4 yfor i=1:n
- @! K2 x5 L  L3 w    if(x(i)<=x0)&&(x(i+1)>=x0)
6 c6 X2 z  p) h5 l, Q        ihdex = i;
+ I3 ?1 Q% c, Y8 v  T        break;
    end
end              %找到x0所在区间
* E( V8 I! h* Q$ z8 U- HA = diag(2*ones(1,n));   %求解m的系数矩阵
u = zeros(n-2,1);
lamda = zeros(n-1,1);
1 b6 R+ S$ g/ L  e+ X/ y- ?9 i: Xc = zeros(n,1);
for i = 2:n-1
+ [  ^* k) y" S: Y2 F, R" Y6 E+ X( Z    u(i-1) = (x(i)-x(x-1))/(x(i+1)-x(i-1));
8 C4 m0 W% s% M9 n  D    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));
4 t" b- V' ]) h  Y+ s3 E A(i,i+1) = u(i-1);
! y9 L- G5 @1 B/ p- q5 U A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
end
c(1) = 2*y_1;
' W5 A" ^& ?5 {5 @6 Q/ uc(n) = 2*y_n;
m = followup(A,c);   %用追赶法求解方程组
. p, f+ D& `) u" T- [+ a: Yh = x(index+1) - x(index);  %x0所在区间长度
f = y(index)*(2*(t-x(index))+h)*(t-x(index+1))^2/h/h/h+...
! X: |7 c/ J8 T! f9 S( `    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处的插值
   
0 \0 f0 U! a/ K* C' _' U9 |+ o   

) K( Y& P- w( `* m( M* G1 J