Carlettos

23.c

Oct 23rd, 2020
744
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. import random as rnd
  2. import matplotlib.pyplot as plt
  3. # import ganas.de.vivir as suicidio
  4.  
  5.  
  6. def browniano(n, μ, σ):
  7.     x = [0]
  8.     y = [0]
  9.  
  10.     for i in range(n):
  11.         dx = rnd.normalvariate(μ, σ)
  12.         dy = rnd.normalvariate(μ, σ)
  13.         x.append(x[i] + dx)
  14.         y.append(y[i] + dy)
  15.         pass
  16.     return x, y
  17.  
  18.  
  19. σ = 1
  20. μ = 0
  21. tiempos = [100, 500, 1000, 5000, 10_000]
  22. k = 100
  23. for n, t in enumerate(tiempos):
  24.     promedio = 0
  25.     for j in range(k):
  26.         x, y = browniano(t, μ, σ)
  27.         # plt.subplot(k, len(tiempos), n*k + 1 + j)
  28.         # plt.plot(x, y)
  29.         # plt.title(f"t = {t}, intento = {j + 1}")
  30.         promedio += sum([x1**2 + y1**2 for x1, y1 in zip(x, y)])/t
  31.     print(f"promedio en t = {t}: {promedio/k}")
  32.  
  33. # plt.show()
  34.  
RAW Paste Data

Adblocker detected! Please consider disabling it...

We've detected AdBlock Plus or some other adblocking software preventing Pastebin.com from fully loading.

We don't have any obnoxious sound, or popup ads, we actively block these annoying types of ads!

Please add Pastebin.com to your ad blocker whitelist or disable your adblocking software.

×