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

iniL = length(ini);
yI = iniL-4;
DbI = iniL-3;
if length(ini) > 7
    eqmatrix = 'flopFlowMatrix';
else
    eqmatrix = 'flopRedFlowMatrix';
end

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

j = 0;
TS = 1;
while TS & (j < (length(yspan)))
    j = j+1;
    ini = [ini(1:yI-1) yspan(j) ini(yI+1:iniL)];
    switch solver
	case 0
        [l,y]=ode45(eqmatrix, lspan, ini);
	case 1
        [l,y]=ode113(eqmatrix, lspan, ini);
	case 2
        [l,y]=ode23s(eqmatrix, lspan, ini);
	case 3
        [l,y]=ode23t(eqmatrix, lspan, ini);
	case 4
        [l,y]=ode23(eqmatrix, lspan, ini);
	end
    [last vars] = size(y);
    %[y(last,1) y(last,yI) y(last,DbI) min(y(:,DbI))]
    [ini(1) yspan(j)]
    TS = ((abs((y(last,DbI) - min(y(:,DbI)))) < prec) & (y(last,yI) >= 0));
end

y0c = yspan(j);
