迭代法求方程的根 （ψ(x)-x=0的根）

function [x_star,k]=my_iterate(fun,a,ep)%用非线性求解迭代法通过 x=ψ(x),求解ψ(x)-x=0的根%N最大迭代次数% a是不动点的近似%ep容许误差N=500my_fun=fun;k=1;x_0=ax_1=feval(my_fun,a)    while abs(x_1-x_0) && N >= 0        k=k+1;        N=N-1        x_0=x_1        x_1=feval(my_fun,x_0)           end    x_star=x_1;end

如求解x=exp(-x)

fun=inline('exp(-x)')

[x_star,k]=my_iterate(fun,0.5)