R = 40 cm

R = 20 cm

R = 10 cm

R= 30 cm

Estimación de la frecuencia en la zona verde.

Ver

¿Qué observamos?

R_2 = 1 (slider from -5 to 5, step 0.1)
R_4 = 4 (slider from 1 to 8, step 1)
n = 630 (slider from 1 to 1000, step 1)
L1 = Sequence((R_4 * sqrt(random() + 0.0 * k); (random() + 0.0 * k) * 360.0°), k, 1.0, n)
R_1 = 3 (slider from 1 to 5, step 1)
score = Sequence(If[x(Element[L1, k])^(2.0) + y(Element[L1, k])^(2.0) ≤ R_2^(2.0), 1.0, 0.0], k, 1.0, n)
A = Intersect(xAxis, yAxis)
c = Circle(A, R_1)
d = Circle(A, R_4)
T1: "Fréquence théorique : " + ((R_1 / R_4)^(2.0))
Freq: Sum[score] / n
a = Sum(score)
R_3 = 2 (slider from 0 to 3, step 0.1)
e = Circle(A, R_3)
f = Circle(A, R_2)
g = Circle(A, R_1)
h = Circle(A, R_4)
k = Circle(A, R_3)
p = Circle(A, R_2)
score2 = Sequence(If[x(Element[L1, k])^(2.0) + y(Element[L1, k])^(2.0) ≤ R_3^(2.0), 1.0, 0.0], k, 1.0, n)
Freq2: Sum[score2] / n
NB2 = Sum(score2)
q = Circle(A, R_4 + 0.3)
B, C = Intersect(q, yAxis)
D = Intersect(d, h, 1.0)
E = Intersect(c, g, 1.0)
H = Intersect(f, p, 1.0)
I = Intersect(d, h, 1.0)
L = Intersect(d, h, 1.0)
M, N = Intersect(h, xAxis)
O, P = Intersect(g, xAxis)
Freq3: R_1^(2.0) / R_4^(2.0)
text11: "f ≅ " + Freq3
b = false
text3: "Número total de lanzamientos: " + n
text2: "Lanzamientos en la zona roja: " + a
text4: "Lanzamientos en la zona amarilla: " + NB2
text5: "Frecuencia experimental: " + Freq2
text1: "Frecuencia experimental : " + Freq