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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
%求已知数据点的第一类三次样条差值多项式及其插值点处的值
%已知数据点的x坐标向量：x
8 e% L0 M6 g! l; G' U%已知数据点的y坐标向量：y
%左端点的一阶导数：y_1
2 R: Y5 g2 F* O: T%右端点的一阶导数：y_n
( p0 f. I4 G* R; P%插值点的x坐标：x0
%求得的三次样条差值多项式：f
, B( f( Q2 o3 m' b/ b. f; _! M%求得的x0处的插值：f0
syms t
0 B4 `) c4 T$ E/ ^" `3 z+ Gf = 0.0000;
, C( f" m6 O, |7 q2 x: F$ hf0 = 0.0000;
1 K3 @' a( b0 v% m' j: V$ S" Dif(length(x)==length(y))
/ S. \1 w2 O2 W! r2 W2 @    n = length(x);
else
    disp('x和y的维数不相等！');
' Q) Z# A! e- `% {" F, U    return;
end              %维数检查
for i=1:n
    if(x(i)<=x0)&&(x(i+1)>=x0)
        ihdex = i;
        break;
3 c# h1 [2 [* o# n- Q7 w& E6 g    end
" G  `' T2 @) i9 s  {end              %找到x0所在区间
A = diag(2*ones(1,n));   %求解m的系数矩阵
u = zeros(n-2,1);
* @) M7 g+ Q; V; ^& O7 y- Flamda = zeros(n-1,1);
3 {  F3 k- w0 D% c2 bc = zeros(n,1);
for i = 2:n-1
    u(i-1) = (x(i)-x(x-1))/(x(i+1)-x(i-1));
  j, f: K" d4 X8 O8 _# @0 D/ C    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));
& k: c/ G8 [1 o$ V A(i,i+1) = u(i-1);
0 w& F; t/ f' ^0 Y0 t A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
end
5 B0 l( g: R& a$ Dc(1) = 2*y_1;
c(n) = 2*y_n;
m = followup(A,c);   %用追赶法求解方程组
) U/ J+ \- z( bh = x(index+1) - x(index);  %x0所在区间长度
f = y(index)*(2*(t-x(index))+h)*(t-x(index+1))^2/h/h/h+...
! L: t6 v8 d" m( b3 t* A6 o    y(index+1)*(2*(x(index+1)-t)+h)*(t-x(index))^2/h/h/h+...
5 J9 ]8 k; a' N3 `    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处的插值
   
5 ]2 Z' s. ^1 F9 r( |   
- E$ L" a- _% |: D/ \$ u- {" u