function flopPlot(l,A,clr)

Kc = A(:,1); vc = A(:,2); Ks = A(:,3); vs = A(:,4); 
y = A(:,5); Db = A(:,6); s = A(:,7); lam = A(:,8);
        
figure(1)
subplot(2,2,1)
hold on;
plot(l,Kc,clr,'LineWidth',1.5);
%title('K_{c}-D_{b} phase space');
xlabel('l');
ylabel('Kc');
grid on;

subplot(2,2,2)
hold on;
plot(l,Db,clr,'LineWidth',1.5);
xlabel('l');
ylabel('Db');
grid on;

subplot(2,2,3)
hold on;
plot(l,Ks,clr,'LineWidth',1.5);
xlabel('l');
ylabel('Ks');
grid on;

subplot(2,2,4)
hold on;
plot(l,y,clr,'LineWidth',1.5);
%plot(l,Ks.*vs-Kc.*vc,strcat('-.*',clr),'LineWidth',1.5)
xlabel('l');
ylabel('y');
grid on;

figure(2);
% subplot(2,1,1)
hold on;
plot(l,lam,clr,'LineWidth',1.5);
xlabel('l');
ylabel('\lambda');
grid on;

% subplot(2,1,2)
% hold on;
% plot(l,s,clr,'LineWidth',1.5);
% xlabel('l');
% ylabel('s');
% grid on;

figure(3)
% subplot(2,1,1)
hold on;
plot(Kc,Db,clr,'LineWidth',1.5);
xlabel('Kc');
ylabel('Db');
grid on;

% subplot(2,1,2)
% hold on;
% plot(Ks,y,clr,'LineWidth',1.5);
% xlabel('Ks');
% ylabel('y');
% grid on;