function rigdPlot(l,A,clr)

Kc = A(:,1); vc = A(:,2); Ks = A(:,3); vs = A(:,4); 
y = A(:,5); Db = A(:,6); w = A(:,7); lam = A(:,8); L = A(:,9);
        
F = 1./sqrt((1-exp(-w)));
g = 1 + (2./w).*log(1+1./F);

r1 = (vs./vc).^(2-Kc);
r2 = (vs./vc).^(-Kc);
r3 = 1+(1./Kc.^2);
r4 = 1-(1./Kc.^2);
r3s = 1+(1./Ks.^2);
r4s = 1-(1./Ks.^2);
f  = -0.25*(r1.*r4.*Kc+r2.*r4s.*Ks);
dw = abs(L.*Db.*(1-cosh(w)).*(1-f).*F.^3);

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,w,clr,'LineWidth',1.5);
plot(l,dw,strcat('-.*',clr),'LineWidth',1.5)
xlabel('l');
ylabel('\omega');
grid on;

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

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