The non-dimensional Navier stokes equations in rotational form for the boundary layer flow are solved using direct numerical simulation. The length scale and velocity scale of the base flow the boundary layer thickness and the inviscid velocity outside the layer are used as the length and velocity scales at the inlet boundary of the computational domain are used as two characteristics to define the flow Reynolds number. The governing equations are discritised in the streamwise direction using a sixth order compact finite difference scheme, and in the cross-direction using a mapped compact finite difference scheme. An algebraic mapping is used to map the physical domain into the computational domain .The compact third order of Runge-kutta method is used for the time-advancement purpose. The convective outflow boundary condition is employed to create a non-reflective type boundary condition at the outlet. The simulation results of this flow were compared by Blasius solution that show accuracy program. In this study, also, characteristics of laminar boundary layer flow verification for accuracy program with divided the lengths and velocity by length of plane and uniform velocity of environment respectively. Profiles and contours of velocity and vorticity have planed in flow arrow and self- similar have seen.