function [ lse_out ] = steady_state_fit_fn(x)
global q un uc Kn Kc Vn Vc r
global yv nv uv cv
global icount

icount = icount + 1;
Vn = x(1);
Vc = x(2);
Kn = x(3);
Kc = x(4);
r  = x(5);

lse_out = 0;

for ii = 3:length(uv)
    
    un      = uv(ii);
        
    ss      = [cv(ii);nv(ii);yv(ii)];
    
    sso     = chemostat_dynamics(0,ss);
    
    lse_out = lse_out + sum(sso.^2);
end

[x' , lse_out]

end