function gaussPhases3(Db,eta,lambda,filename,clr,solver)
y0min = 0;
y0max = 1.2;
Kc0min = 2;
Kc0max = 2.5;
Ksteps = 4;
Kinc = (Kc0max-Kc0min)/Ksteps;
yinc = 0.01;
lf = 150;
N = 3;

for j = 1:Ksteps
    Kc0(j) = Kc0min + Kinc*j;
    ini = vertcat([Kc0(j); 1; y0max; eta; lambda],Db*ones(2*N,1));
    y0c(j) = gaussSepDuch3(ini,y0max-y0min,yinc,lf,solver);
    j
end

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

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