求Gauss插值多项式

function f = GaussInter(x,y,x0)
%求已知数据点的Gauss插值多项式
if(length(x) == length(y))
    n = length(x);
else
    disp('x和y的维数不相等！');
    return;
end

xx =linspace(x(1),x(n),(x(2)-x(1)));
if(xx ~= x)
    disp('节点之间不是等距的！');
    return;
end

if( mod(n,2) ==1)                  %n为奇数时用斯特林插值
    if(nargin == 2)
        f = GStirling(x,y,n);
    else if(nargin == 3)
        f = GStirling(x,y,n,x0);
        end
    end
else                              %n为偶数时用贝赛尔插值
    if(nargin == 2)
        f = GBessel(x,y,n);
    else if(nargin == 3)
        f = GBessel(x,y,n,x0);
        end
    end
end