An accurate analysis procedure to solve the flow about a flat plate at various incidences has been developed. The Navier-Stokes equations of stream function and vorticity form are solved in a sufficiently large computational domain, in which the grid lines are mutually orthogonal. The details of the flow near the singularity at the tip of the plate is well captured by the analytic solution which is asymptotically matched to the numerically generated outer solution. The solution for each region is obtained iteratively : the solution of one (inner or outer) region uses that of the other as the boundary condition after each cycle. The resulting procedure is accurate everywhere and also computationally efficient as the singularity has been removed. It is applied to the flat plate for a wide range of Re : the results agree very well with the existing computation and experiment.