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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
4 d' y: e* c1 {% g  u: b%求已知数据点的第一类三次样条差值多项式及其插值点处的值
%已知数据点的x坐标向量：x
2 D. ~( a3 N% [" `; @: S1 M- X, |%已知数据点的y坐标向量：y
%左端点的一阶导数：y_1
%右端点的一阶导数：y_n
; M* ], D. M1 E* l9 }; |7 p%插值点的x坐标：x0
! h( d. ^* Q' y4 ?, }) s4 Q+ C%求得的三次样条差值多项式：f
%求得的x0处的插值：f0
syms t
f = 0.0000;
f0 = 0.0000;
if(length(x)==length(y))
    n = length(x);
/ N/ V0 k' K' e" {else
    disp('x和y的维数不相等！');
    return;
6 w( f+ Q2 J( E+ b+ |  Pend              %维数检查
$ i! @  D) ]3 {# l, ^for i=1:n
2 s+ R" x4 v  x% T; J    if(x(i)<=x0)&&(x(i+1)>=x0)
3 s- y+ w1 b. v! H        ihdex = i;
        break;
    end
. E/ L2 B) a' e# S8 W. _' Dend              %找到x0所在区间
: c9 P$ A5 B) t; oA = diag(2*ones(1,n));   %求解m的系数矩阵
u = zeros(n-2,1);
/ Z* K) u/ @. a2 m* @! Plamda = zeros(n-1,1);
c = zeros(n,1);
for i = 2:n-1
' M$ E5 Q0 W% _7 n0 P    u(i-1) = (x(i)-x(x-1))/(x(i+1)-x(i-1));
, Y' j$ B% t( c3 R  D( J5 t/ ?4 r5 {    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# S4 {& }/ a A(i,i-1) = lamda(i);   %构造系数矩阵及向量c
end
c(1) = 2*y_1;
c(n) = 2*y_n;
m = followup(A,c);   %用追赶法求解方程组
* K) e" O: E; D  F/ Ih = x(index+1) - x(index);  %x0所在区间长度
1 N1 C+ Q& d* V# Ff = y(index)*(2*(t-x(index))+h)*(t-x(index+1))^2/h/h/h+...
' p% R+ L" ], L( |    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-...
% U. V' d3 q( D0 d/ T* ~    m(index+1)*(x(index+1)-t)*(t-x(index))^2/h/h; %x0所在区间的插值函数
5 {6 f9 M0 u7 ]5 U  J/ Hf0 = abs(f,'t',x0)  %x0处的插值
   
   
! Q6 I" c8 O5 i; Q9 i% ?