gamma = 0.5 # external stimulus amplitude
kappa = 0.2 # frequency of external stimuli
# Initialize variables
t_r = np.linspace(0, T, N)
psi = np.zeros(N) # consciousness function
R = np.zeros(N) # time resonance function
# Initial value of awareness
psi[0] = 0.1
# Derivative functions of consciousness
def dpsi_dt(psi_t, t):
  return -beta * psi_t + gamma * np.sin(kappa * t)
Beri Komentar
Belum ada komentar. Jadilah yang pertama untuk memberikan komentar!