The conventional shape function for the finite-element method (FEM) is a linear one, which is inadequate for analyzing numerically complex flows for high Reynolds numbers. In the present study, we propose GSMAC-CIP scheme using 3rd order shape function which requires continuity of the value of the function and its first space derivative in whole space and is formulated a finite element method for the cubic interpolated pseudo particle (CIP) method. We verified the effectiveness of this new scheme by the forced convection in a square cavity at Re=1 000, 5 000, 10 000. The numerical results obtafined by the present scheme is compared with those of GSMAC-FEM using rough meshes and it is showed that present scheme is superior to GSMAC-FEM on the space accuracy. Moreover, it is shown that the numerical results obtained by present scheme using fine meshes are agreement with those obtained by Ghia et al.