The coupled diffusion equations recently proposed for concentrated hard-sphere suspensions of interacting Brownian particles, the nonlinear deterministic diffusion equation with the self-diffusion coefficient DS(F(x, t)) for the average local volume fraction F(x, t), and the linear stochastic diffusion equation with DS(F(x, t)) for the density fluctuations dn(x, t) are numerically solved under a spatially inhomogeneous, nonequilibrium initial state. Thus, in a supercooled region where Fß = F g, the slow evolution of the cluster-like glassy domains with F(x, t) = Fg and the slow relaxation of the nonequilibrium density fluctuations are shown to be caused by the dynamic singularity of the self-diffusion coefficient, DS(F(x, t)) ~ (1-F(x, t)/Fg)2, where F is a particle volume fraction, Fg = (4/3)3/(7 ln 3 - 8 ln 2 + 2) is the colloidal glass transition volume fraction, and Fß is the crossover volume fraction.