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

潇湘夜雨123456

function [f,f0] = ThrSample1(x,y,y_1,y_n,x0)
3 w1 j$ W# t: P* \2 }0 c: I%求已知数据点的第一类三次样条差值多项式及其插值点处的值
%已知数据点的x坐标向量：x
) @: S7 Z- k4 {' {%已知数据点的y坐标向量：y
( b8 L1 |% `7 `9 E& q%左端点的一阶导数：y_1
%右端点的一阶导数：y_n
' ?5 O, q0 E" k. m( o* c+ @%插值点的x坐标：x0
%求得的三次样条差值多项式：f
$ U, q: J0 z/ Q" `2 _0 U%求得的x0处的插值：f0
# D4 E& B# T* \; i6 Asyms t
f = 0.0000;
f0 = 0.0000;
5 E7 i$ V# w7 _9 Q" ?, zif(length(x)==length(y))
    n = length(x);
2 y, C4 L7 l$ l& {$ c. }else
( J# }- T3 i# M3 ]    disp('x和y的维数不相等！');
    return;
end              %维数检查
2 _5 V$ L# g  M2 J) N. ofor i=1:n
: w9 c8 m9 h5 w$ r! g    if(x(i)<=x0)&&(x(i+1)>=x0)
        ihdex = i;
  t4 S& K+ {, D8 \7 Q. u: f        break;
    end
end              %找到x0所在区间
8 g0 ?1 u- w( i& `- {# {) GA = diag(2*ones(1,n));   %求解m的系数矩阵
: {) W$ a1 n' D6 D5 @6 n1 Au = zeros(n-2,1);
6 {) e3 O* O6 j5 i3 m8 S3 Flamda = zeros(n-1,1);
3 x# ^1 A7 o6 M$ q7 Uc = zeros(n,1);
; f7 E$ ?) @3 F0 X: |. t9 N- _for 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
5 S/ A2 |8 c) jend
8 F+ t  ^0 l! t6 Cc(1) = 2*y_1;
c(n) = 2*y_n;
m = followup(A,c);   %用追赶法求解方程组
, R4 t* h0 h2 H0 p  Dh = x(index+1) - x(index);  %x0所在区间长度
1 E2 b+ N, \1 n3 o4 Z5 if = 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-...
( o) k3 T; P! M+ e& n: ~& J    m(index+1)*(x(index+1)-t)*(t-x(index))^2/h/h; %x0所在区间的插值函数
$ ?" i* {" Y! X+ ?3 q) m, ^f0 = abs(f,'t',x0)  %x0处的插值
   
" ^8 V/ k) k- k% B8 I  b# `   
  g" D1 W2 s0 g7 V- i- ~( g' _
6 ~5 [7 h$ ]* r# q1 J& j% F