# Fungsi untuk tekanan internal P
def pressure(rho, w):
  return w * rho
# Fungsi untuk menghitung H(t)
def H_dot(rho, P, E, kappa, H):
  # Fungsi untuk menghitung perubahan kecepatan ekspansi
  return (8 * np.pi * G / 3) * (P - kappa * H**2) # Simpifikasi
# Parameter awal
G = 6.67430e-11 # Konstanta gravitasi
rho_0 = 1e-9 Â # Kerapatan energi awal
w = 0.33 Â # Parameter keadaan, misalnya untuk radiasi
Beri Komentar
Belum ada komentar. Jadilah yang pertama untuk memberikan komentar!