# Euler Simulation
for i in range(1, N):
  psi[i] = psi[i-1] + dpsi_dt(psi[i-1], t_r[i-1]) * dt
  phi = 1 + alpha * psi[i]
  R[i] = np.cos(omega * t_r[i] * phi)
# Plot the simulation results
plt.figure(figsize=(12, 6))
plt.subplot(2,1,1)
plt.plot(t_r, psi, label='Consciousness $\psi(t_r)$')
plt.xlabel('Relative Time $t_r$')
plt.ylabel('Consciousness')
Beri Komentar
Belum ada komentar. Jadilah yang pertama untuk memberikan komentar!