The boundary element method was used for incompressible viscous fluid flow past a circular cylinder near a plane. The various phases of the flow passing the cylinder are studied numerically with different g/d. The numerical simulation of velocity fields is in good agreement with Davis' theoretical results and Taneda's experiments. In particular, the distributive curves of the surface force and the pressure on the circular cylinder were calculated numerically.