1667877a-9fc4-4a16-b61c-1aa34d0b287a 28 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

'N' random in critical region N 1000 :

t0 1 time :

data spec pts i 1 data rows range p i el a b c spec : for v p 1 el : Unest 'p' j 2 p rows range v v p j el stack : for data v augment 51 line :

size/define critical region.
Plot first before filling xmin 1 - : xmax 5 : ymin 5 - : ymax 0 : 22 mat space_x xmax xmin - (: space_y ymax ymin - (: 21 line 21 sys

A function that tests an arbitrary point 
for its membership in the critical region. x y Test x 2 y * + 1 ≤ (5 x * 2 y * - 5 ≥ (* x 2 y * - 9 ≤ (* :

data spec pts ~ 12 times faster ≡

A routine that finds an array of N random 
pairs in the critical region. Pts i 1 N range x 1 Random space_x * xmin + 1 el : y 1 Random space_y * ymin + 1 el : 1 x y Test - x 1 Random space_x * xmin + 1 el : y 1 Random space_y * ymin + 1 el : 21 line while P i 1 el x : P i 2 el y : 51 line for P 21 line :

a b c spec . 6 #808080 13 mat : fill Pts a b c spec pts : 21 line

Filling a range of a PDF

1 x - 2 / 5 5 x * - ( 2 / - 9 x - ( 2 / - 3 1 sys fill 2 1 sys

1 time t0 - 21 s *

1 x - 2 / 55 x * - (2 / - 9 x - (2 / - 31 sys fill 21 sys

Solving for the X's, Y's Apex

I hate FindRoot ! f init GivenFind f init FindRoot :

f1 1 x - 2 / 55 x * - (2 / - : (0 ≡ f2 55 x * - (2 / 9 x - (2 / - : (0 ≡ f3 9 x - (2 / 1 x - 2 / + : (0 ≡ 31 sys

common 'x' init x 1 ≡ 11 mat :

x sol f1 init GivenFind : X1 x sol : 21 line x sol f2 init GivenFind : X2 x sol : 21 line x sol f3 init GivenFind : X3 x sol : 21 line 31 sys

Apex1 X1 1 X1 - 2 / 12 mat : Apex2 X2 55 X2 * - (2 / - 12 mat : Apex3 X3 9 X3 - (2 / - 12 mat : 31 line

Coordinates of each Apex of the critical region Apex Apex1 Apex2 Apex3 stack : 111 mat 011 mat 1 - 11 mat 5 - 11 mat 511 mat 2 - 11 mat 32 mat

Validate the search in the critical domain x y region x 2 y * + 1 ≤ (5 x * 2 y * - 5 ≥ (* x 2 y * - 9 ≤ (* :

1 ≡ valid search
0 ≡ invalid search 22 - region 2 - 2 - region 21 mat 1021 mat