A method is presented for calculating neutron cross sections averaged over the partial width fluctuations in the unresolved resonance range. The N dimensional integral (N ≤ 4) which appears in the averaging procedure is mathematically reduced to integration in a single dimension. The resultant expression is numerically evaluated by a Gaussian quadrature technique.