The neutron diffusion parameters of water and ice were measured by the pulsed source method at two temperatures; 1.0°C. and -19°C. Neutron pulses were obtained at one millisecond intervals by modulating the beam in a Cockcroft-Walton type accelerator. The ₁H³(d,n)₂He⁴ reaction was used to generate neutrons. The samples were contained in cylindrical aluminum cans covered with cadmium. The experiment was conducted inside a large paraffin block which served as a neutron shield and thermal insulator. The temperature of the samples was maintained constant to within ±1°C. Neutrons leaving one surface of the sample were counted in a BF₃ proportional counter. The time distribution of these neutrons was recorded by an eighteen channel time analyzer. The width of each channel was 20 microseconds. The opening of the first channel was delayed 100 microseconds with respect to the start of the neutron burst to minimize harmonics in the neutron decay. A geometric buckling was calculated for each sized sample from B²=[2.405/(R+∈)]²+[π/(H+∈)]² where B² = geometric buckling 2.405 = first zero of J_o Bessel Functions R = radios of cylinder H = height of cylinder ∈ = extrapolation distance The extrapolation distance ∈ was calculated from ∈ =0.71 λ_tr where λ_tr = mean free path of neutrons in water The extrapolation distance was assumed to vary as T^½ where T is the temperature in degrees Kelvin. The measured decay constants, α, were fitted by the method of least squares to a polynomial in B² of the form α = (∑_av) ÷ D_oB² - CB⁴ where ∑_a = the macroscopic absorption cross-section v = the neutron velocity D_o = diffusion coefficient C = diffusion cooling coefficient The resultant values of (∑_av) and D_o for each temperature are below. The data did not permit a determination of C. 1.0°C. (∑_av) = 4595 ± 365 sec⁻¹ D_o = 29600 ± 840 cm²/sec -19°C. (∑_av) = 4355 ± 263 sec⁻¹ D_o = 27050 ± 630 cm²/sec / Ph. D.