kn = 5000;
  nn = 60;
  randn('state',sum(100*clock))
    rnv = unifrnd(0,10,kn,nn);
for kk =1:kn
 rn = rnv(kk,:);
sumexp(1) = rn(1);
for ii=2:nn
 sumexp(ii) = sumexp(ii-1)+rn(ii);
end;
 j1 = max(find(sumexp<60));
 j2 = min(find(sumexp>60));
 if (isempty(j1)==1) 
 dif1(kk) = 0;
else
 dif1(kk) = sumexp(j1);
end;
 dif2(kk) = sumexp(j2);
end;
avg =sum((dif2-dif1))./kn;