In this paper boundary elements method (BEM) is equipped with finite difference approximation (FDA) to solve two-dimensional Navier–stokes (N–S) equation. The N–S equation is converted to a system of ordinary differential equations (ODEs) respect to time in streamfunction-vorticity formulation.