fdifference.m

这是我从一本基于MATLAB实现的数值分析算法的书中自带的光盘程序

%2-2 向前差分表
function f=fdifference(xi,fi)
xi=[1,3,5,7,9,11,13];
fi=[1.0,0.5,0.3333333,0.25,0.2,0.1666666,0.14285714];
n=length(fi);
deltaf=zeros(size(xi));
fprintf('xi=%12.8f\n',xi);
fprintf('fi=%12.8f\n',fi);
for k=n-1:-1:1
    fprintf('k=%12.8f\n',k);
    for j=1:k
        deltaf(j)=fi(j+1)-fi(j);
        fi(j)=deltaf(j);
        fprintf('deltaf(j)=%12.8f\n',deltaf(j));
   end 
end
return