reset

p = 0.25 (slider from 0.1 to 0.9, step 0.05)
n = 1600 (slider from 100 to 2000, step 10)
wurm = Sequence(If[random() < p, 1, 0], i, 1, n)
wurmrh = Sequence(Sum[wurm, i] / i, i, 1, n)
ywerte: wurmrh
xwerte = Sequence(1, n)
punkte = Sequence((Element[xwerte, i], Element[ywerte, i]), i, 1, n)
Q = Element(punkte, 100)
R = Element(punkte, 400)
P = Element(punkte, 25)
Liste1: {}
o(x) = p + (2 * sqrt((p * (1 - p) / x)))
u(x) = p - (2 * sqrt((p * (1 - p) / x)))
wurmflop = Sequence(If[u(i) < Element[wurmrh, i] < o(i), 0, 1], i, 1, n)
wurmimtrichter = If(Sum[wurmflop] ≟ 0, 1, 0)
listewurmtrichter: {1, 1, 1, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0, 0}
wurmimtrichterah = Sum(listewurmtrichter)
wurmimtrichterrh: wurmimtrichterah / Length[listewurmtrichter]
wurmimtrichterrhdezimal: round(wurmimtrichterrh, 3)
S = Element(punkte, 1600)
Simulationen = Length(listewurmtrichter)
Text1: "
Simulationen: " + Simulationen + "
Anteil im Wurzeltrichter geblieben:  " + wurmimtrichterrh + ""
rh100: (Sum[wurmflop, 100] * 0.01)
rh400: (Sum[wurmflop, 400] * 0.0025)