Model Kosmologi Berbasis Fisika Ketupat Lebaran

# Loop untuk melakukan integrasi numerik

for i in range(1, len(time)):

    H = H_values[i-1]

    rho = rho_0  # Asumsi kerapatan konstan untuk simplifikasi

    P = pressure(rho, w)

    H_values[i] = H_values[i-1] + H_dot(rho, P, E_0, kappa_0, H) * dt

# Plot hasil simulasi

plt.plot(time, H_values)

plt.xlabel('Waktu (t)')

plt.ylabel('Kecepatan Ekspansi H(t)')

plt.title('Simulasi Ekspansi Semesta')