Several scales' density fluctuations which exist in the early universe will first gravitationally collapse along one axis and make pancake-like structures. If the collapsed baryonic pancake heats up over 104K by shock formation, radiative cooling begins to work and mass accretion toward the central region will advance. Because of this effect, mass fraction of the high density layer becomes large. Densities and widths of the layers will reflect masses of structures (e.g. galaxy) which will be formed after caustics. In this respect, we assumed an Einstein-de Sitter universe dominated by cold dark matter (ΩDM = 0.9) and investigated the evolutions of fluctuations numerically using one-dimensional hydrodynamic plus N-body codes. We applied a new method for larger fluctuation scales; it is a hybrid method of Eulerian PPM and Zeldovich approximation and it can simulate around the central pancake region with high accuracy.