Problem with decision logic for 5 layers with the forward modeling method

Good afternoon everyone. Can anyone answer this question for me? I have a problem with decision logic for 5 layers with the forward modeling method.

The syntax below is for 3 layers. Would you mind for helping me with the algorithm of a 5 layer forward modeling?

Thank you.

x_hole_1 = 0 ; x_hole_2 = 1000; y_hole = 500;
y_batas = [100 220 320 420 500]; 
v_lap = [1000 1500 2000 2500 3000];
sou_y = [50 80]; 
rec_y = [50 300 480]; 

%%%FUNGSI  HITUNG  TIME	MODEL (DIRECT )
%%%menghitcng panjang direct ray dari source ke receiver 
n_ray = length(sou_y) * length(rec_y) ;
pan_ray_total = zeros(length(sou_y),length(rec_y));
for i = 1 : length( sou_y)
    for j = 1 :length(rec_y)
        pan_ray_total(i,j)= sqrt((x_hole_1 - x_hole_2).^2 + (sou_y(i) - rec_y(j)).^2) ;
    end
end
%%menghitung sudut tembak ray dari titik shoot
sudut_tembak = zeros(length(rec_y),length(sou_y));
for i = 1 : length(sou_y)
    for j = 1 :length(rec_y)
        sudut_tembak(i,j) = acosd((x_hole_2/pan_ray_total(i,j))) ;
    end
end
%%%menghitung panjang  tiap segmen ray yang melewati tiap blok
k = 1;
for i = 1 :length(sou_y) 
    for  j = 1 :length(rec_y)
        %%%dasar pembutaan syntax memakai simulasi sthooting dari source 1
        tikpot_bts1 = ((y_batas(1)-sou_y(i))/tand(sudut_tembak(i,j)))+ x_hole_1; 
        pan_ray_bag(k,1) = sqrt((x_hole_1 - tikpot_bts1).^2 + (sou_y(i) - y_batas(1)).^2);
        tikpot_bts2 = ((y_batas(2) - sou_y(i))/tand(sudut_tembak(i,j))) + x_hole_1;
        pan_ray_bag(k,2) = sqrt((tikpot_bts1 - tikpot_bts2).^2 + (y_batas(1) - y_batas(2)).^2); 
        pan_ray_bag(k,3) = pan_ray_total(i,j) - (pan_ray_bag(k,1) + pan_ray_bag(k,2));
        %%%loop hubungan kondisi posisi source dan receiver 
        if (sou_y(i) < y_batas(1)) && (rec_y(j) < y_batas(1))
            %%%direct ray di blok 1
            pan_ray_bag(k,1) = pan_ray_total(i,j);
            pan_ray_bag(k,2) = 0;
            pan_ray_bag(k,3) = 0;
        elseif (sou_y(i)< y_batas(1)) && (rec_y(j) > y_batas(1)) && (rec_y(j) < y_batas (2)) 
            pan_ray_bag(k,2)= pan_ray_total(i,j)- pan_ray_bag(k,1);
            pan_ray_bag(k,3) = pan_ray_total(i,j) -(pan_ray_bag(k,1) + pan_ray_bag(k,2));
        elseif (sou_y(i) > y_batas(1)) && (sou_y(i) < y_batas(2)) && (rec_y(j) < y_batas (1))
            pan_ray_bag(k,2) = pan_ray_bag(k,1);
            pan_ray_bag(k,1) = pan_ray_total(i,j) - pan_ray_bag(k,1);
            pan_ray_bag(k,3) = pan_ray_total(i,j) - (pan_ray_bag(k,1) + pan_ray_bag(k,2));
        elseif (sou_y(i) > y_batas(1)) && (sou_y(i) < y_batas(2)) && (rec_y(j) > y_batas(1)) && (rec_y(j) < y_batas(2))
        %%%direct ray di blok 2 
            pan_ray_bag(k,1) = 0;
            pan_ray_bag(k,2) = pan_ray_total(i,j);
            pan_ray_bag(k,3) = 0;
        elseif (sou_y(i) > y_batas(1)) && (sou_y(i) < y_batas(2)) && (rec_y(j) > y_batas(2)) 
            pan_ray_bag(k,1) = 0;
            pan_ray_bag(k,2) = sqrt((x_hole_1 - tikpot_bts2).^2 + (sou_y(i) - y_batas(2)).^2);
            pan_ray_bag(k,3) = pan_ray_total(i,j) - (pan_ray_bag(k,1) + pan_ray_bag(k,2));
        elseif (sou_y(i) > y_batas(2)) && (rec_y(j) < y_batas(1))
            pan_ray_bag(k,3) = sqrt((x_hole_1 - tikpot_bts2).^2 + (sou_y(i) - y_batas(2)).^2);
            pan_ray_bag(k,2) = pan_ray_bag(k,1) - pan_ray_bag(k,3);
            pan_ray_bag(k,1) = pan_ray_total(i,j) - abs(pan_ray_bag(k,2) + pan_ray_bag(k,3));
        elseif (sou_y(i) > y_batas(2)) && (rec_y(j) > y_batas(1)) && (rec_y(j) < y_batas (2))
            pan_ray_bag(k,1) = 0;
            pan_ray_bag(k,2) = sqrt((x_hole_1 - tikpot_bts2).^2 + (sou_y(i) - y_batas(2)).^2);
            pan_ray_bag(k,3) = pan_ray_total(i,j) -(pan_ray_bag(k,1) + pan_ray_bag(k,2));
        elseif (sou_y(i) > y_batas(2)) && (rec_y(j) > y_batas(2))
            %%%direct ray di blok 3
            pan_ray_bag(k,1) = 0;
            pan_ray_bag(k,2) = 0;
            pan_ray_bag(k,3) = pan_ray_total(i,j) - (pan_ray_bag(k,1)+pan_ray_bag(k,2));
        end
    k = 1 + k;
    end
end
%%%menghitung matriks forward 
dt_cal = pan_ray_bag *(1./v_lap');
for l = 1 : length (sou_y)
    d(l,:) = dt_cal((length(rec_y)*(l-1))+1 : length (rec_y)*l,1);
end
G  = pan_ray_bag ;