用变步长中点方法和三点求导公式求exe ^ xex在x=1初期的导数值,比较步长变化对解的影响
functiont=var_size(fname,h )格式长度=10; forI=0:km=h/(2^I; t=feval(fname,m ); DISP([I,T] )结束名称=在线) ) (exp )-exp(1-x ) )/(2* x ) ) ); h0=1; h1=0.1; h2=0.01; var_size(fname,h0 ) var_size(fname,h1 ) var_size(fname,h2 )取小数点后五位,如果h=1,则结果在第九次迭代时收敛; 另一方面,当h=0.1时,结果在第五次迭代中收敛; h=0.01时,第二次开始接近真值。 由此可见,h越小,收敛速度越快。