% script to run maxEig function 200 times
%% Configure parpool
myCluster = parcluster('local');
myCluster.NumWorkers = str2double(getenv('SLURM_NTASKS'));
myCluster.JobStorageLocation = getenv('TMPDIR');
myPool = parpool(myCluster, myCluster.NumWorkers);

count = 200;
dim = 5001;
sumMaxe = 0;

tic
parfor i=1:count;
  sumMaxe = sumMaxe + maxEig(i,dim);
end
toc
avgMaxEig = sumMaxe/count

delete(myPool);
exit