In this paper, a method of forming stiffness matrices for a quarter of slab and subgrade has been developed, which may be used for analysing arbitrarily loaded slabs on elastic solid foundation with finite element method.