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

function [x_star,k]=my_iterate(fun,a,ep)


%用非线性求解迭代法通过 x=ψ(x),求解ψ(x)-x=0的根
%N最大迭代次数
% a是不动点的近似
%ep容许误差


N=500
my_fun=fun;
k=1;
x_0=a
x_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)