text: This article provides an error analysis of time discretization applied to spatially discrete approximation of the stationary and nonstationary Navier-Stokes equations. The nonlinearity of the convection term is the main problem in solving a stationary or nonstationary Navier-Stokes equation or Euler equation problems.  Stoke incorporated ‘the method of artificial compressibility’ to solve these problems.