-
import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint def sfr_reheat(y, t, a11, a12, a21, a22, b1, b2, Pd): omega, z = y dydt = [a11*omega + a12*z + b1*Pd, a21*omega + a22*z + b2*Pd] return dydt H = 5 D = 1 Tr = 8 Fh = 0.3 R = 0.05…
-
Hello world!
Welcome to User’s blog Sites. This is your first post. Edit or delete it, then start writing!