Matlab数值分析编程：牛顿下山法解方程

%通用代码
fun=inline('x^0.5-x^3+2','x');%手动输入函数表达式和自变量
dfun=inline('0.5*x^(-0.5)-3*x^2','x');%手动输入函数一阶导数表达式和自变量
x0=1.5;  %手动输入初值
eps=0.00001;%手动输入题目允许最大误差

%以下无需修改
Lambda=1;
xk=x0;
xk_plus=xk-Lambda*(feval(fun,xk))/(feval(dfun,xk));%因为循环语句中需要判断循环条件，所以先计算一次xk+1
flag=0;%引入一个标志位，下山因子过大时置1
while (abs(xk_plus-xk)>=eps)
    if flag==0
      xk=xk_plus;
    end;
      xk_plus=xk-Lambda*(feval(fun,xk))/(feval(dfun,xk));%牛顿下山法迭代
      Dk=xk;
    if abs(feval(fun,xk_plus))>=abs(feval(fun,xk))%判断下山因子是否满足条件
      Lambda=Lambda*0.5;
      flag=1;
    end;
    if flag==1
      xk=x0;
    end;
    flag=0;
end
Answer=xk_plus;