We formulate linear and nonlinear exact and approximate mathematical computational models for spherical and hemispherical densities and fluxes of solar radiation in scattering and absorbing inhomogeneous plane-parallel layers. In particular, we consider Rayleigh and conservative layers. We place emphasis on the formulation of boundary conditions for these models.