#highest of 50n observations from geometric

$q=17/18;
$p=1/18;
$n=6000;
$s=$n*50;

for $i (0..6000) {
$j=$i+1;
print $i." ".(1-$q**$j)**$s."\n";
}