方程求根（牛顿迭代法）

function x = agui_newton(fname,dfname,x0,e)
N=10;
x=x0;
x0=x+2*e;
k=0;
while abs(x0-x)>e&&k<N
    k=k+1;
    x0=x;
    x=x0-feval(fname,x0)/feval(dfname,x0);
    disp(x)
end
if k==N 
    warning('已达最大迭代次数');end

end

命令行

>> fun=inline('exp(x)+10*x-2')

fun =

     内联函数:
     fun(x) = exp(x)+10*x-2

>> x=agui_bisect(fun,0,1,0.5e-3)

x =

   0.090332031250000