func13.asv

来自「数值分析最常用的四十种算法」· ASV 代码 · 共 11 行

x=linspace(1,8,8)
y=[15.3,20.5,27.4,36.6,49.1,65.6,87.8,117.6];
y = log(y);
[p,S] = polyfit(x,y,1)
f = polyval(p,x);
y=exp(y)
f=exp(f)
plot(x,y,'o',x,f,'-')
disp(abs(norm(y,inf)-norm(f,inf)))

%the error