The critical magnetic scattering has been investigated in EuS by means of small-angle scattering with polarized neutrons using an inclined magnetic field geometry, allowing the determination of three-spin correlation functions. Two contributions to the critical magnetic scattering I ⌺ ͑q͒ = I ↑ ͑q͒ + I ↓ ͑q͒ and ⌬I͑q͒ = I ↑ ͑q͒ − I ↓ ͑q͒ were studied for temperatures near T C = 16.52 K. The I ↑ ͑q͒ and I ↓ ͑q͒ are the scattering intensities for the incident neutron beam polarized along ͑↑͒ and opposite ͑↓͒ to the magnetic field. The symmetric contribution, namely I ⌺ ͑q͒, comes from the pair-spin correlation function. The scattering intensity is well described by the Ornstein-Zernike expression I ⌺ ͑q͒ = A͑q 2 + 2 ͒ −1 , where = −1 is the inverse correlation length of the critical fluctuations. The correlation length obeys the scaling law = a 0 − , where = ͑T − T C ͒ / T C is the reduced temperature, a 0 = 0.17 nm, and = 0.68± 0.01. The difference contribution ⌬I͑q͒ is caused by the three-spin chiral dynamical spin fluctuations that represent the asymmetric part of the polarization dependent scattering. The q dependence of ⌬I͑q͒ follows closely 1 / q 2 . ⌬I͑q͒ depends on the temperature as − with = 0.64± 0.05. The exponents as determined by means of the static measurements by and the dynamic measurements ͑using the chirality͒ are in excellent agreement with each other, demonstrating the internal consistency of the theory and the experiment. Therefore, our results confirm the principle of the critical factorization, which is known as Polyakov-Kadanoff-Wilson operator algebra.