计算方法单点截弦matlab实现

 function Dan(fun,a,b,e)
%  使用单点截弦法求方程 fun=0 在区间[x1,x2]处的根 e为误差限，默认为 eps
%  对于任意的 x 属于 [a,b] ， diff(fun,x) 与 diff(fun,x,2) 连续且不变号
%  选择x0为不动点，xn为动点， 要求 fun(x_0)*diff(fun,x_0,2) > 0
%  其中迭代公式为 y=a-(x-a)/(fun(x)-fun(a))*fun(a)

% 作者：野渡无人

% 最后修改日期：2008.3.31

% Dan('x^3-2*x-5',2,3,10^(-5))
% n=  0       x0= 3.00000000000000000      
% n=  1       x1= 2.00000000000000000      
% n=  2       x2= 2.05882352941176450    
% n=  3       x3= 2.08126365984502290    
% n=  4       x4= 2.08963921009084740    
% n=  5       x5= 2.09273957431800550    
% n=  6       x6= 2.09388370846184820    
% n=  7       x7= 2.09430545112525970    
% n=  8       x8= 2.09446084576648770    
% n=  9       x9= 2.09451809338575720    
% n= 10       x10= 2.094539182