One-velocity angular flux of a neutron field generated by a polarly anisotropic point source in an infinite homogeneous region is calculated in P3approximation, using the spherical harmonics operator formalism adapted to two-dimensional spherical geometry. The problem is reduced to simple algebraic equations. For weak absorption and linearly anisotropic scattering, the solution is given in closed form. (orig.).