In finite element method, the order of complete polynomial of the interpolation function is related to the number of nodes in the element. This paper presents a four-node quadrilateral element with quadratic function on it. The presented displacement functions maintain C0-continuity. Meanwhile, the element stiffness matrix is derived from the displacement functions. Test problems show that high accuracy can be achieved by the use of the new displacement function on the element.