smain.m

Lagrange插值、Newton插值函数等插值函数

%--首先，要求建立方程所需的各个系数；先求这几个区间；再依次求得各个系数。最后求解方程，带入方程求解。

function   smain(n)

hold off 
x=[-5:0.01:5];
y=f(x);

sx=[-5:10/n:5];
ys1=f(sx);

h=cdd(sx);    % 等分个个区间，求出系数h；
[u,v]=xsm(h); % 根据已经求得的系数求解u,v；
d=cs(ys1,h);  % 根据已知的节点函数值和h求得d；
m=solve(u,v,d); % 带入边界条件进行求解，得系数m；
ys=lq(x,sx,ys1,h,m);  %根据m,节点，以及节点函数值插值；


plot(x,y,'r',x,ys,'b')
hold on
plot(sx,ys1,'go')


xlabel('x')
ylabel('y=1/1+x2')