function fmax=quest1_SteadyState_util(x, pname, eq_ind, exs_)
global M_

if nargin==0,
  x=[];
end
for j=1:length(x),
set_param_value(pname{j},x(j));
end

ys = quest1_steadystate;

exs = M_.exo_nbr;
ex_ = zeros(1,exs);

yc = quest1_static(ys, ex_, M_.params);
% fmax = yc(10);
if nargin<3,
  [fmax, iff] = max(abs(yc));
  if length(x)>1,
    fmax=log10(fmax);
  else
    fmax=yc(iff);    
  end
else
  fmax=yc(eq_ind);
end