function fFlow(ini,filename,clr);
Y = flopEvalFlow(ini,100,0);
S = size(Y);

l = Y(:,1); Kc = Y(:,2);
Db = Y(:,S(2)-3); Df = Y(:,S(2)-2);
s = Y(:,S(2)-1); lambda = Y(:,S(2));

if strcmp(filename,'nf') == 0
	dlmwrite(filename,Y,' ');
end

hold on;
plot(Kc,Db,clr);
title('K_{c}-D_{b} phase space');
xlabel('K_{c}');
ylabel('D_{b}');
grid on; hold on;

figure;
subplot(2,1,1)
hold on;
plot(l,s,clr);
title('s(l)');
xlabel('s');
ylabel('l');
grid on;

subplot(2,1,2)
hold on;
plot(l,lambda,clr);
title('\lambda (l)');
xlabel('l');
ylabel('\lambda');
grid on;
