function [ lse_out ] = simple_fit_fn_2(x) global q uy uz Vmax Ky Kz eps delta global yv nv uv cv global icount icount = icount + 1; Vmax = x(1); Ky = x(2); Kz = x(3); delta = x(4); uz = x(5); lse_out = 0; us = (1-Kz*delta/(Ky*eps))*q*Ky/(Vmax-q) +(delta/eps)*uz; wx = mean(yv); wn = mean(nv); for ii = 1:length(uv) uy = uv(ii); ss = [yv(ii);nv(ii)]; if uy