67a08bca-ed3a-4912-9289-81332ceab1c1 15 4 5 decimal

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

N 10000 : X 4 N 1 Random (* 2 - : Y 4 N 1 Random (* 2 - : 31 line

r 1 : R 1.5 : 21 line

t0 1 time :

Disk i 1 N range U i el X i el X i el (2^Y i el (2^+ r 2^≥ (X i el (2^Y i el (2^+ R 2^≤ (& 0 cases : V i el Y i el X i el (2^Y i el (2^+ r 2^≥ (X i el (2^Y i el (2^+ R 2^≤ (& 0 cases : 21 line for U V augment 21 line :

Disk Disk eval :

R 0 + 100 blue 1 5 mat Disk 2 1 sys

1 time t0 - 36 s *

a# b# start end N S u v N var u v u v u - N / + (range : uv start end N var : j 1 uv rows range y j el i 2 * π * uv j el * exp eval : for a# y Re * b# y Im * augment 41 line :

start 0 : end 1 : N 200 : 13 mat a R : b R : α r : β r : 14 mat ρ 0.5 r R + (* : 31 sys

x y axes x y / - :

x a b start end N S 1 col : y a b start end N S 2 col : circle_1 x y augment : u α β start end N S 1 col : v α β start end N S 2 col : circle_2 u v augment : ω1 ρ ρ start end 12 S 1 col : ω2 ρ ρ start end 12 S 2 col : rod ω1 ω2 augment : rod rod 112 range 12 range el : 101 line

Disk circle_1 circle_2 axes rod 5 1 sys

1 time t0 - 46 s *

Washer style ...

R 1.5 :

r 1 :

rad start end N S u v N var u v u v u - N / + (range : uv start end N var : j 1 uv rows range y j el i 2 * π * uv j el * exp eval : for rad y Re * rad y Im * augment 41 line :

start 0 : end 1 : N 200 : 13 mat a R : α r : ρ 0.5 r R + (* : 13 mat 21 sys

x a start end N S 1 col : y a start end N S 2 col : circle_1 x y augment : u α start end N S 1 col : v α start end N S 2 col : circle_2 u v augment : ω1 ρ start end 12 S 1 col : ω2 ρ start end 12 S 2 col : rod ω1 ω2 augment : rod rod 112 range 12 range el : 101 line

j 1 rod rows range Ω j el rod j 1 el rod j 2 el o 10 black 15 mat : for

plot circle_1 circle_2 Ω 1 el Ω 2 el Ω 3 el Ω 4 el Ω 5 el Ω 6 el Ω 7 el Ω 8 el Ω 9 el Ω 10 el Ω 11 el Ω 12 el 141 sys :

plot