A new method for determining various quantities describing the radiation field in an inhomogeneous, plane-parallel atmosphere is proposed in this two-part paper. The essence of this method is the reduction of the boundary value problems which arise during the customary statement of various astrophysical problems associated with solving the radiative transfer equations to initial value problems. Compared to previous attempts in this area, the proposed method is universal and simple. The first part of this paper deals with one-dimensional media. Scalar, as well as vector-matrix problems relating to the diffusion of radiation in spectral lines with frequency redistribution are examined.