Constrained random points on a circle: Difference between revisions

m
Line 3,491:
val lim = ( ~350.0/900.0, ~225.0/900.0 ) ; (* limits to usum *)
val select = fn i => usum i <= #2 lim andalso usum i >= #1 lim ; (* select according to inequalities *)
val uv = ListPair.zip ( urandomlist (#1 realseeds) 2500 , urandomlist (#2 realseeds) 2500 ) (* drawtake 2500 couples *)
in
List.take ( List.filter select uv , 1000 )
Anonymous user