The M1 strengths (or level density of 1+ states) are of importance for estimation of interaction strengths between neutrinos and nuclei for the study of the supernova neutrino-process. In 1957, Agodi predicted theoretically angular distribution of neutrons emitted from states excited via dipole transitions with linearly polarized gamma-ray beam at the polar angle of θ=90° should be followed by a simple function, a + b cos(2φ), where φ, is azimuthal angel. However, this theoretical prediction has not been verified over the wide mass region except for light nuclei as deuteron. We have measured neutron angular distributions with (polarized gamma, n) reactions on Au, Nal, and Cu. We have verified the Agodi's prediction for the first time over the wide mass region. This suggests that (polarized gamma, n) reactions may be useful tools to study M1 strengths in giant resonance regions.