/* ev.tc - expected value - 1/30/10 - tct
/* ev=20
/* linux interpreter required

int last
rand [ /* returns random number in 1..2^31-1
int m;m=2147483647
int a;a=16807	
last=last*a%m
if (last<0) last=-last
return last
]

ev [
int i,c,j,k
/* 107374182.4 = (2^31-1)/20
/* see ev.html, evc.html
pl "";ps "ev.tc - 1/30/10 - tct";pl ""
pl "";ps "seed? ";last=gn;pl ""
for (k=1;k<=10;++k) [
 c=0
 for (j=1;j<=1000;++j) [
  for (i=1;i<=1000;++i) [
   if (rand()<107374182) [c=c+i;break]
   ]
  ]
 ps "expected value * 1000 = ";pn c;pl ""
 ]
]