Need help plotting across more than one layer

Here is the code I currently have. However, I need to plot my first layer from x = 0:0.1:0.4 and my second layer from x = 0.4:0.1:1.

My first layer equation needs to be y(x) = B(1)C(x)a(1) where a(1) = inv(B(1))*V(:,1).

For layer two, my equation needs to be y(x) = B(2)C(x)a(2) where a(2) = inv(C(2))*inv(B(2))*T(1)*B(1)*a(1)

clear all;
format long;
im = sqrt(-1);
CellLength = 1;
ibeta = 1;

%Define materal properties 
CellLength = 1;
layers = 2;
d = [0.4;0.6];
dTotal = d(1,1)+d(2,1);

xc = [0;0.4];

Ef = 12;
pf = 3;
cf = sqrt(Ef/pf);

Em = 1;
pm = 1;
cm = sqrt(Em/pm);

w = 5;

T1 = [cos(.2*w) (1/(6*w))*sin(.2*w); -6*w*sin(.2*w) cos(.2*w)];
T2 = [cos(.6*w) (1/w)*sin(.6*w); -w*sin(.6*w) cos(.6*w)];
T = T2*T1;

Z1 = 6*w;
Z2 = w;
Z = [Z1;Z2];

%Solve eigenvalue problem for k
[V,D] = eig(T);                 %D = eigenvalues, %V = eigenvectors
k1 = log(D(1,1))/(im*dTotal);
k2 = log(D(2,2))/(im*dTotal);
k = [k1;k2];

for j = 1:layers
  B = @(j)([1 1; im*Z(j,1) -im*Z(j,1)]);
  B = B(j);
  C_a = @(j)([exp(im*k(j,1)*xc(j,1)) 0; 0 exp(-im*k(j,1)*xc(j,1))]);
  C_a = C_a(j);

      if j == 1
          a = (inv(B)*V(:,1));
          alpha = a;
          beta = B;
      else
          a = inv(C_a)*inv(B)*T*beta*alpha;
      end

      L=[];
      M=[];

      for x = 0.5:0.1:5.5
          C = @(x)([exp(im*k(j,1)*x) 0; 0 exp(-im*k(j,1)*x)]);
          C = C(x);
          y = @(x)(B*C*a);
          y = y(x);

          L(end+1)=y(1);
          M(end+1)=y(2);
      end
  end

plot(0.5:0.1:5.5,real(L))
hold on
plot(0.5:0.1:5.5,real(M))
hold off
xlabel('Position, x')
ylabel('Displacement, u')