We consider multiple resonance scattering with complete frequency redi
stribution (CFR) in a semi-infinite conservative atmosphere (photon de
struction probability epsilon(1) = 0) with the sources at infinite dep
th. The polarization arising in resonance scattering is completely acc
ounted for. The problem we consider is the resonance-scattering counte
rpart of the Chandrasekhar-Sobolev problem of Rayleigh scattering in t
he: conservative atmosphere. The numerical data on the matrix source f
unction S(tau) in the atmosphere with conservative dipole resonance sc
attering (the depolarization parameter W = 1 are presented: we assume
Doppler profile. The source matrix is found by a non-iterative numeric
al solution of the matrix Wiener-Hopf integral equation with the matri
x A-operator. Depth dependence of the elements of the source matrix S(
tau) is discussed. Some unexpected peculiarities are revealed in the b
ehavior of its polarization terms. The matrix I(z) which is the genera
lization of the Chandrasekhar H-function to the case of polarized reso
nance scattering is found by the iterative solution of the Chandrasekh
ar-type nonlinear matrix integral equation. We present high-accuracy (
5 s.f.) numerical data on I(z) for dipole conservative scattering with
the Doppler profile. The center-to-limb variation of the degree of po
larization in the core of a Doppler broadened resonance line is found.
In conservative case. the limiting limb polarization delta(0) in the
core of such a line is 9.4430% (for W = 1). The dependence of bo on th
e depolarization parameter W is found. Simple interpolation formula, d
elta(0) = (9.443 - 38.05 root epsilon(1))%, is suggested for the limb
polarization of the radiation emerging from an isothermal nearly conse
rvative atmosphere (epsilon(1) much less than 1. W = 1). The data on I
(z) are used to find the polarization line profiles and to trace their
center-to-limb variation. The asymptotic expansions of S(tau) for tau
--> infinity (deep layers) and of I(z) for z --> infinity (line wings
) are found for the case of the Doppler profile. The coefficients of t
he expansions are determined by recursion relations. The numerical dat
a on the accuracy and the domain of applicability of the asymptotic th
eory lire presented.