# Generate flux integers
n_i = poisson.rvs(mu=5, size=100)
n_i = np.clip(n_i, -10, 10)
# Solve inflaton dynamics
phi = [phi0]
for _ in range(10000):
dW = np.random.normal(0, np.sqrt(dt))
dphi = (-V_prime(phi[-1])/(3*H) * dt) + (H**(3/2)/(2*np.pi)) * dW
phi.append(phi[-1] + dphi)
```
Simulation Conclusion
Beri Komentar
Belum ada komentar. Jadilah yang pertama untuk memberikan komentar!