/* slow loop - counting elements > 5.99 */

x = -2*ln(rndu(1000000,3));


tim1=hsec;
count1 = 0;
i = 1;
do until i > rows(x);
	j = 1;
	do until j > cols (x);
		if (x[i,j] > 5.99); 
			count1 = count1 +1; 
		endif;
		j = j+1;
	endo;
	i=i+1;
endo;
tim1 = hsec - tim1;

tim2=hsec;
count2 = sumc(vecr(x .> 5.99));
tim2 = hsec-tim2;

"times, count" 0.01*(tim1~tim2) 100*(count1~count2)/(rows(x)*cols(x));