We consider the Stokes equationin a domain Ω in R n (n ≥ 2) with smooth boundary. Here u = (u 1 , . . . , u n ) are unknown velocity field and p is unknown pressure field. Initial data u 0 is assumed to satisfy a compatibility condition : div u 0 = 0 in Ω and the normal component of u 0 equals zero on ∂Ω. This system is a typical parabolic equation and it has several properties resembling the heat equation.If Ω = R n , u is reduced to a solution of the heat equation with initial data u 0 because there is no boundary condition. For example, a regularity-decay estimateholds for all 1 ≤ p ≤ ∞ with C independent of t and u 0 , where f (t) p := Ω |f (t, x)| p dx 1/p and ∇ denotes the gradient in the space variables. If p = 2, the estimate (0.2) is still valid for any domain.