Abstract. In this study, a finite volume method (FVM) is suitably used for solving the problem of a fully coupled fluid flow in a rectangular domain with slip boundary conditions. Numerical solutions for the flow variables, viz. velocity, and pressure have been computed. The FVM, with an upwind scheme, has been implemented to discretize the governing equations of the present problem. The well known SIMPLE algorithm is employed for pressure-velocity coupling. This was executed with the aid of a computer program developed and run in a C-compiler. Computations have been performed for unknown variables with Reynolds numbers (Re) = 50, 100, 250, 500, 750 and 1000. The behavior of steady-state solutions of velocity and pressure of the fluid along horizontal and vertical through geometric center of the rectangular domain have been illustrated. We observed that, with the increase of the Reynolds number, the absolute value of velocity components decreases whereas the pressure value increases.MSC 2010: 35Q30, 65M08, 65N08, 74S10, 76D05, 76M12