  We use the BBGKY hierarchy equations to calculate, perturbatively, the lowest
order nonlinear correction to the two point correlation and the pair velocity
for Gaussian initial conditions in a critical density matter dominated
cosmological model. We compare our results with the results obtained using the
hydrodynamic equations which neglect pressure and we find that the two match,
indicating that thare are no effects of multistreaming at this order of
perturbation. We analytically study the effect of small scales on the large
scales by calculating the nonlinear correction for a Dirac delta function
initial two point correlation. We find that the induced two point correlation
has a $x^{-6}$ behaviour at large separations. We have considered a class of
initial conditions where the initial power spectrum at small $k$ has the form
$k^n$ with $0 < n \le 3$ and have numerically calculated the nonlinear
correction to the two point correlation, its average over a sphere and the pair
velocity over a large dynamical range. We find that at small separations the
effect of the nonlinear term is to enhance the clustering whereas at
intermediate scales it can act to either increase or decrease the clustering.
At large scales we find a simple formula which gives a very good fit for the
nonlinear correction in terms of the initial function. This formula explicitly
exhibits the influence of small scales on large scales and because of this
coupling the perturbative treatment breaks down at large scales much before one
would expect it to if the nonlinearity were local in real space. We physically
interpret this formula in terms of a simple diffusion process. We have also
investigated the case $n=0$ and we find that it differs from the other cases in
certain respects. We investigate a recently proposed scaling property of
