计算方法习题四4.6

若f(x)=x^7 + x^3+1，求f[1,2,4…128]和f[1,2,4…128,256]

clear
clc
x=[1 2 4 8 16 32 64 128 256]
y=[3 137 16449 2097665 268439553 3.435977113700000e+10 4.398046773249000e+12 5.629499555184650e+14 7.205759405470515e+16]
deltx=diff(x);
delty=diff(y);
firstorder=delty./deltx %一阶
for i=1:length(x)-2
    delt2x(i)=x(i+2)-x(i);
end
delt2y=diff(firstorder);
secondorder=delt2y./delt2x %二阶
for i=1:length(x)-3
    delt3x(i)=x(i+3)-x(i);
end
delt3y=diff(secondorder);
thirdorder=delt3y./delt3x %三阶
for i=1:length(x)-4
    delt4x(i)=x(i+4)-x(i);
end
delt4y=diff(thirdorder);
fourorder=delt4y./delt4x %四阶
for i=1:length(x)-5
    delt5x(i)=x(i+5)-x(i);
end
delt5y=diff(fourorder);
fif=delt5y./delt5x
for i=1:length(x)-6
    delt6x(i)=x(i+6)-x(i);
end
delt6y=diff(fif);
six=delt6y./delt6x
for i=1:length(x)-7
    delt7x(i)=x(i+7)-x(i);
end
delt7y=diff(six);
sev=delt7y./delt7x
for i=1:length(x)-8
    delt8x(i)=x(i+8)-x(i);
end
delt8y=diff(sev);
eig=delt8y./delt8x