In this paper we investigate the stiffness matrix of a general quadrilateral element in closed form using n x n Gauss-Legendre quadrature rule. For this, we propose four types of nodal coordinate transformation. The terms of the matrix are divided into two groups, namely - diagonal and non-diagonal. Only one term (called leading) from each group is computed, and then the remaining fourteen terms are computed from these two leading terms exploited one of the proposed types of coordinate transformation. This leads us a great savings in computational time and memory space. In order to compute the matrix we use these transformations in two ways, and thus two algorithms are given to generate the matrix. Finally, numerical example is given to verify the effectiveness of the present formulation. Keywords: Gauss-Legendre quadrature; Numerical integration; Quadrilateral finite element; Stiffness matrix; Closed form. DOI: http://dx.doi.org/10.3329/bjsir.v46i3.9050 BJSIR 2011; 46(3): 399-405