We introduce an accurate non-Hermitian Schrodinger-type approximation of Bloch optical equations for two-level systems. This approximation provides a complete description of the excitation, relaxation and decoherence dynamics in both weak and strong laser fields. In this approach, it is sufficient to propagate the wave function of the quantum system instead of the density matrix, providing that relaxation and dephasing are taken into account via automatically-adjusted time-dependent gain and decay rates. The developed formalism is applied to the problem of scattering and absorption of electromagnetic radiation by a thin layer comprised of interacting two-level emitters.