Network analysis PS#1, 6a, 6b

PS#1 6a-b

datdir = '/Users/clinpsywoo/Documents/2011-2016yr/2014-2015(4th_GS)/Network_modeling/Problem_sets/FB100 dataset/facebook100';
datfiles = filenames(fullfile(datdir, '*mat'), 'absolute');

for i = 1:numel(datfiles)
    
    fprintf('\nWorking on %02d/%02d...\t', i, numel(datfiles)); 
    tic; 
    
    load(datfiles{i});
    
    % number of vertex
    n = length(A);
    
    % mean degree <ku>
    mean_ku(i) = sum(sum(A,2))/n;
    
    % mean degree <kv>
    mean_kv(i) = sum((A * sum(A,2))./(sum(A,1)'))/n;
    
    expected_net_MND1(i) = sum(sum(A,2).^2)/sum(sum(A,2));
    
    expected_net_MND2(i) = sum(A * sum(A,2))/sum(sum(A));
    
    t = toc;
    fprintf('Done in %04dmin %02dsec.\n', floor(t/60),rem(t,60)); 
end

My first upload in gist, thanks to Sanghee.