notas-fnys/deca2.py

22 lines
451 B
Python
Raw Normal View History

2023-10-05 23:04:53 +02:00
from random import random
import matplotlib.pyplot as plt
import numpy as np
#Constantes
NTl = 1000 #numero inicial de thalios
tau = 3.053*60 #vida media del thalio en seg.
#los puntos a graficar
temps = []
#Determinando tiempo de decaimiento
for i in range(NTl):
t = -tau*np.log(1-random()/np.log(2))
temps.append(t)
#La grafica
(counts, bins, patches) = plt.hist(temps)
plt.xlabel("Tiempos")
plt.ylabel("Eventos")
plt.show()