this is part of the program that I am working on and I have a problem when I run the program ... can anyone tell me a good solution for my program?

dt_cal = pan_ray_bag *( 1 ./ v_lap') ;
for l = 1 : length ( sou_y )
    d(l,:) = dt_cal((length(rec_y)*(l-l))+l : length(rec_y) * l,l ) ;
end
a = max (d);
b = min (d);
dt_noise = rand(1)*(a-b)*0.21;
t_obs = d+dt_noise;
%%%  membuat kecepatan awal model linier
vsem = zeros (length(sou_y), length(rec_y) ); 
for i =1 :length(sou_y)
    for j = 1 :length( rec_y)
        vsem(i,j) = pan_ray_total(i,j) ./ t_obs (i,j); 
    end
end

when running....i have a problem like that

Attempted to access t_obs(2,1); index out of bounds because size(t_obs)=[1,3].

Error in inversi_linier (line 14)
        vsem(i,j)= pan_ray_total(i,j)./t_obs(i,j);

Error in SA (line 18)
[ v0_linier,dv_linier,v_linier ] = inversi_linier(d,n_ray,sou_y,rec_y,pan_ray_total,G);

[KALYAN ACHARJYA]

Can you post the complete code with all data?

dt_cal 1D array or 2D?

[No Freed]

%%%FUNGSIHITUNG TINE MODEL	(DIRECT)
%%%menahitung 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=acosd(x_hole_2./pan_ray_total);
    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 penbutaan syntax memakai sImulasi shooting dari source 1
        tikpot_btsl=((y_batas(1)-sou_y(i))/tand(sudut_tembak(i,j))) + x_hole_1;
        pan_ray_bag(k,1)=sqrt((x_hole_1 - tikpot_btsl)^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_btsl - 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))
            %%%direct  ray di  blak	2
            pan_ray_bag(k,1)=0;
            pan_ray_bag(k,2)=pan_ray_tota1(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_tota1(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_tota1(i,j) - (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_tota1(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  blak	3
            pan_ray_bag(k,1)=0;
            pan_ray_bag(k,2)=0;
            pan_ray_bag(k,3)=pan_ray_tota1(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-l))+l : length(rec_y) * l,l ) ;
end
G =pan_ray_bag ;
end

[No Freed]

sorry sir... i"m new in matlab

[No Freed]

and this the input argument

x_hole_1 = 0; x_hole_2 = 100; y_hole = 50;
y_batas = [20 30 50]; 
v_lap = [1000 1200 1500]; 
sou_y = [5 10]; 
rec_y = [5 25 40]; 
n_ray = length(sou_y) * length(rec_y);