dN/dt = R - N/τ - σcNφ; dφ/dt = cσNφ - φ/τ_c + spontaneous; threshold population inversion N_th = 1/(σcτ_c).