晕!
慢慢!待我将相关MATLAB程序贴到我BLOG

function [x,flag]=SecDiv(fun,a,b,ep)
%二分法求方程之根
%求根函数：fun
%求根区间：[a,b]
%精度设置：ep
flag=1;%求解成功
fa=fun(a);
fb=fun(b);
if a>=b |fa*fb>0
    flag=0;
    x=inf;
   return ;
end

while abs(b-a)/2>=ep
    x=(a+b)/2;
    fx=fun(x);
    if fa*fx<0
        b=x;
        fb=fx;
    else
        a=x;
        fa=fx;
    end
end
x=(a+b)/2;

function x=SecDivDev(fun,a,b,ep,k)
%二分法求方程之根
%求根函数：fun
%求根区间：[a,b]
%精度设置：ep
%可能存在的根数:k
x=zeros(1,k);
step=1/k;
ab=a:step:b;
n=length(ab);
m=0;
for i=1:n-1
[xx,flag]=SecDiv(fun,ab(i),ab(i+1),ep);
if flag==1
    m=m+1;
    x(m)=xx;
    if m==k
        break;
    end
end
end
x=x(1:m);

function x=Steffensen(fun,x0,ep)
%Steffensen法求方程之根
%求根函数：fun
%与根相近点：x0
%精度设置：ep
x=x0;
x1=x0;
while 1
x1=x;
y=fun(x)+x;
z=fun(y)+y;
x=x-(y-x)*(y-x)/(z-2*y+x);
if abs(x-x1)<ep
    break;
end
end