function rigidPhases3(Db,filename,clr,solver)
y0min = -0.5;
y0max = 1.5;
Kc0min = 1.75;
Kc0max = 2.15;
Ksteps = 25;
Kinc = (Kc0max-Kc0min)/Ksteps;
yinc = 0.01;
lf = 150;

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

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

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