function F = steadystate1(x)
global U1;global U2;
X = [U2-U1*(x(1)+x(1)^2)];
F = sqrt(X'*X);
F = F*1.e4;

