r = 10000
kreistreffer = 0
          quadrattreffer = (2 * r + 1) * (2 * r + 1)
for y = -r to r
    for x = -r to r
        if x * x + y * y <= r * r then
            kreistreffer = kreistreffer + 1
            ausgabe 4 * kreistreffer / quadrattreffer