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 ) (* |
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 ) |