st.m

来自「无约束优化中的共轭梯度算法程序」· M 代码 · 共 15 行

function s=st(x,eps)
if nargin==1
    eps=1.0e-6;
elseif nargin<1
    error
    return
end
xx=myfun(x);
while norm(xx-x)>=eps
    x=xx;
    xx=myfun(x);
end
s=xx;
return