function y0c = statSepDuch3(ini,range,step,lf)
prec = 0.001;

iniL = length(ini);
yI = iniL-2;
DbI = iniL-1;
if length(ini) > 5
    eqmatrix = 'statFlowMatrix';
else
    eqmatrix = 'statRedFlowMatrix';
end

lspan = 0:0.1:lf;
yspan = fliplr(ini(yI)-range:step:ini(yI));

j = 0;
delocal = 1;
while delocal
    j = j+1;
    ini = [ini(1:yI-1) yspan(j) ini(yI+1:iniL)];
    [l y] = ode45(eqmatrix, lspan, ini);
    [last vars] = size(y);
    [y(last,1) y(last,yI) y(last,DbI) min(y(:,DbI))]
    delocal = abs((y(last,DbI) - min(y(:,DbI)))) < prec
end

y0c = yspan(j);