The dual reciprocity boundary element method (DRBEM) formulation is presented for solving incompressible magnetohydrodynamic (MHD) flow equations. The combination of Navier-Stokes equations of fluid dynamics and Maxwell's equations of electromagnetics through Ohm's law is considered in terms of stream function, vorticity and magnetic potential in 2D. The velocity field and the induced magnetic field can be determined through the relations with stream function and magnetic potential, respectively. The numerical results are visualized for several values of Reynolds (Re), Hartmann (Ha) and magnetic Reynolds number (Rem) in a lid-driven cavity, and in a channel with a square cylinder. The well-known characteristics of the fluid flow and MHD flow are exhibited. These are the shift of the core region of the flow and the development of the main vortex in the vorticity through the center of the cavity as Re increases. An increase in Ha causes Hartmann layers for the flow at the bottom and top walls. Higher values of Rem result in circulation of the magnetic potential at the center of the cavity. An increase in Re causes symmetric vortices behind the cylinder to elongate through the channel, and an increase in Hartmann number suppresses this elongation.