In order to calculate laminar skin friction more accurately, the laminar flat-plate skin friction was calculated by solving the N-S equations. For solving the N-S equations, 4 advanced upwind schemes, 14 different grid distributions and 5 limiters were employed. All the results were compared with Blasius solution. The results show that the van Lee′s FVS(flux vector splitting)scheme has considerable numerical dissipation, it is also sensitive to grid quality and limiter, so it is unsuitable for friction calculation. The solution is more accurate with a higher limiter and scheme viscous resolution. The calculation precision strongly depends on the height of the first level mesh. Considering both of the calculation ability and the grid aspect ratio, the height of the first level grid can be reduced to improve precision. For skin friction calculation, there must be no less than 4 grid points in the boundary layer.