An overview of the computational methods implemented in a two-dimensional laminar flow solver is presented. The methods discussed include coordinate transformations making the code capable of solving flow in non-rectilinear domains, a discretization method implemented in the computational domain, and a pressure coupling method which is used to enforce the continuity equation. Results of the numerical solver for laminar flow are presented and discussed.