matlab如何请输入插值节点n,功能强大的牛顿插值matlab程序代码

format long;

way_in = input('请选择输入的内容(1或2):\n1、输入为f(x)表达式，区间[a,b]及其等分数n的值\n2、输入为f(x)表达式和插值点横坐标xi的值\n');

switch way_in

case 1

f = input('请输入函数表达式:f(x) = ', 's');

a = input('请输入区间左端值a：');

b = input('请输入区间右端值b：');

n = input('请输入区间等分值n：');

np = input('请输入插值函数在区间内绘图点数(默认输入100)：');

for i=1:n+1

x(i) = a + (b-a)/n*(i-1);

y(i,1) = eval(subs(f,'x(i)','x'));

end

for j=1:n

for k=j:n

temp=y(k+1,j)-y(k,j);

y(k+1,j+1)=temp/(x(k+1)-x(k+1-j)) ;

end

c(j)=y(j,j);

end

c(j+1)=y(j+1,j+1);

for k=1:np-1

xx(k)= a + (b-a)/np*k;

yy(k) = eval(subs(f,'xx(k)','x'));

end

for k=1:np-1

xs=xx(k);

for i=1:n+1

if i==1

s(i)=c(i);

else

s(i)=c(i);

for j=1:i-1

s(i)=s(i)*(xs-x(j));

end

end

end

Nn(k)=sum(s);