fixposys.m

来自「various matlab programs to slove various」· M 代码 · 共 18 行

function [alpha, res, nit]=fixposys(dim, x0, nmax, toll, Phi, F)
x = x0; alpha=[x']; res = 0;
for k=1:dim, 
    r=abs(eval(F(k,:))); if (r > res), res = r; end 
end;
nit = 0; residual(1)=res;
while ((nit <= nmax) & (res >= toll)),
   nit = nit + 1; 
   for k = 1:dim, xnew(k) = eval(Phi(k,:)); end
   x = xnew; res = 0; alpha=[alpha;x]; x=x';
   for k = 1:dim, 
       r = abs(eval(F(k,:)));  
       if (r > res), res=r; end, 
    end
    residual(nit+1)=res;
 end
 res=residual';