function [mean_dt_array, diff_dt] = mean_dt(dt1,dt2)



% [mean_dt_array, diff_dt] = (dt1,dt2)

%

% Note that the two input arrays dt1 and dt2 must be of the same size

%

% outputs:

%

% mean_dt_array - mean packet departure time difference [sec]

% diff_dt - difference in input dt arrays [sec]

% std_out - display successful queries per hop




[M,N]=size(dt1);

mean_dt_array = zeros(M,1);

dt_array = (dt1 - dt2);

disp(dt_array)



% Find lost packets given by NaN

lost = isnan(dt_array(:,:));

dt_array(lost) = -1;

disp(dt_array)



for I=1:M,

    i=0;

    dt_sum=0;

    for J=1:N,

        % If not lost packet

        if (dt_array(I,J) ~= -1)

           
i=i+1;

           
dt_sum = dt_sum + dt_array(I,J);

        end

    end

    % Calculate mean for this row of queries

    mean_dt_array(I,1)=(dt_sum/i);

    disp(i)

end