rosen.m

来自「LMF方法解非线性方程组」· M 代码 · 共 8 行

function r = rosen(x)
% Rosenbrock's valey with a constraint
R = sqrt(x(2)^2+x(1)^2)-.5;
% Residuals:
r = [10*(x(2)-x(1)^2) % first part
1-x(1) % second part
(R>0)*R*1000 % penalty
];