Constrained random points on a circle: Difference between revisions

Content added Content deleted
Line 3,491: Line 3,491:
val lim = ( ~350.0/900.0, ~225.0/900.0 ) ; (* limits to usum *)
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 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 ) (* draw 2500 couples *)
val uv = ListPair.zip ( urandomlist (#1 realseeds) 2500 , urandomlist (#2 realseeds) 2500 ) (* take 2500 couples *)
in
in
List.take ( List.filter select uv , 1000 )
List.take ( List.filter select uv , 1000 )