360Hz

#include #include #include // Define complex number type using namespace std; using complexd = complex; // Quantum coherence (wavefunction at time t) complexd psi_total(double t, vector& a, vector& E, vector& phi, double hbar) { complexd psi = 0.0; for(size_t n = 0; n > rho_t(double t, vector>& rho0, vector>& Gamma) { size_t N = rho0.size(); vector> rho(N, vector(N)); for(size_t i = 0; i 1e-6; // small threshold for numerical check } // Probability distribution from wavefunction double probability_density(complexd psi) { return norm(psi); // |psi|^2 }