matlab中的牛顿下山法

function [gentimer]=xiashan（faTOL1TOL2rminN）
%TOL2为残量精确度，TOL1为根的误差限
df=diff（sym（f））;
x1=a;
time=0;
r=1; %r为下山因子
if abs（subs（fx1））    gen=x1; return 
end
while time    time=time+1;
    r=1;
    df1=subs（dfx1）;
    fx=subs（fx1）;
    if df1==0
        error
    end
    gen=x1-r*fx