df2(j+2)=(y(j)-4*y(j+1)+3*y(j+2))/(2*10); end

%*****************增长率计算属于拟合函数的一阶导数 syms z1;

p1=(c(1)+c(2)*z1)./(1+c(3)*z1); pf=diff(p1,'z1',1); p2=inline(pf,'z1'); subplot(122)

plot(z,p2(z),'mp',x,df1,'g--',x,df2,'r*',x,df2,'r-'); legend('拟合函数图像','二点公式','三点公式'); title('人口年增长率'); function err=fitfun2(d,x,y) a=d(1); b=d(2); c=d(3);

err=y.*(1+c*x)-a-b*x;

err=err'*err;