function a = example1_serial()

% adjust m and N according to your resources 
N = 2000; 
m = 50; 
a = zeros(N,1); 
seeds = (1:N);

tic; % serial (regular) for-loop 
for I = 1:N 
    a(I) = largestEigenvalue(m,seeds(I)); 
end
t = toc;

disp(['Time for serial calculation: ', num2str(t), ' seconds.'])
hist(a)