  We study one-particle spectra and a two-particle correlation function in the
  130 GeV/nucleon Au+Au collisions at RHIC by making use of a hydrodynamical
model. We calculate the one-particle hadronic spectra and present the first
analysis of Bose-Einstein correlation functions based on the numerical solution
of the hydrodynamical equations which takes both longitudinal and transverse
expansion into account appropriately. The hydrodynamical model provides
excellent agreement with the experimental data in the pseudorapidity and the
transverse momentum spectra of charged hadrons, the rapidity dependence of
anti-proton to proton ratio, and almost consistent result for the pion
Bose-Einstein correlation functions. Our numerical solution with simple
freeze-out picture suggests the formation of the quark-gluon plasma with large
volume and low net-baryon density.
