function statPhases3(Db,filename,clr)
y0min = -0.7;
y0max = 1.4;
Kc0min = 2;
Kc0max = 3;
Ksteps = 30;
Kinc = (Kc0max-Kc0min)/Ksteps;
yinc = 0.005
lf = 100;

for j = 1:Ksteps
    Kc0(j) = Kc0min + Kinc*j;
    ini = [Kc0(j) 1 y0max Db 0];
    y0c(j) = statSepDuch3(ini,y0max-y0min,yinc,lf);
end

if strcmp(filename,'nf') == 0
	dlmwrite(filename,[Kc0' y0c'],' ');
end

if strcmp(clr,'np') == 0
    plot(Kc0,y0c,clr);
    grid on;
end