function [pmax] = bootfun(submax);
pmax= nansum(submax)./length(submax);
end