1 .高斯插值
功能f=gauss (x，y，x0 ) )。
if (长度(x )==长度(y ) y ) )
n=长度(x；
else
DISP(x和y的维数不相等！ ' )；
返回；
结束
xx=linspace(x(1)、x (n )、(x ) )2- x (1) )；
if(xx )=x ) ) )
disp '节点之间不是等间隔的！' )；
返回；
结束
if(mod(n，2 )==1) ) ) ) ) ) )。
if(nargin==2) ) ) )。
f=gstirling(x，y，n )；
(elseif(nargin==3) ) ) ) ) ) ) )。
f=gstirling(x，y，n，x0 )；
结束
结束
else
if(nargin==2) ) ) )。
f=gbessel(x，y，n )；
(elseif(nargin==3) ) ) ) ) ) ) )。
f=gbessel(x，y，n，x0 )；
结束
结束
结束
功能f=gstirling (x，y，n，x0 ) )。
syms t；
nn=(n1 )/2；
f=y(nn；
for(I=1:n-1 )。
for(j=I1:n )。
y1(j )=y ) j )-y ) j-1；
结束
if(mod(I，2 )==1) ) ) ) )。
c(I )=) y1(In )/2 (y1 ) ) In2 )/2 )/2；