Langkah 6.1: Definisikan Fungsi-fungsi yang Diperlukan
Fungsi untuk kerapatan energi: Menyatakan perubahan kerapatan energi dalam semesta berdasarkan hukum termodinamika dan elastisitas.
Fungsi untuk tekanan internal: Menyatakan hubungan antara PPP dan \rho berdasarkan persamaan keadaan P=wP = w \rhoP=w.
-
Fungsi untuk elastisitas: Fungsi untuk menghitung respons elastisitas materi pengisi dan elastisitas wadah ruang-waktu.
Fungsi untuk kecepatan ekspansi H(t)H(t)H(t): Berdasarkan persamaan Friedmann dan elastisitas.
Langkah 6.2: Implementasi Metode Euler atau Runge-Kutta
Misalnya, dengan metode Euler, kita akan melakukan iterasi untuk menghitung nilai H(t)H(t)H(t) untuk setiap langkah waktu t\Delta tt.
python
SalinEdit
import numpy as np
import matplotlib.pyplot as plt