Numerical simulations of an unsteady laminar lid-driven cavity Newtonian flow are executed to demonstrate formation of eddies under different aspect ratios and steady-state criteria. A couple of nonlinear unsteady partial differential equations (PDEs) of Navier Stokes satisfying a set of boundary conditions is examined with both inclusions of no-slip and free-slip effects. Finite difference method (FDM) is employed to solve the vorticity-stream function using a self-developed Matlab® code embedded in a GUI to ease the cavity problem analysis by the end users. Four models are studied where Ansys Fluent® employing finite element method is used to verify the present FDM steady-state results produced for the first model. It is observed that the value of stream function begins to drop and streamline distribution changes shape when