The Integrated Force Method (IFM) has been developed in recent years for the analysis of civil, mechanical and aerospace engineering structures. In this method all independent or internal forces are treated as unknown variables which are calculated by simultaneously imposing equations of equilibrium and compatibility conditions. The solution by IFM needs the computation of element equilibrium and flexibility matrices from the assumed displacement, stress-resultant fields and material properties. This paper presents a general purpose code for the automatic generation of element equilibrium and flexibility matrices for plate bending elements using the Integrated Force Method. Kirchhoff and the Mindlin-Reissner plate theories have been employed in the code. Paper illustrates development of element equilibrium and flexibility matrices for the Mindlin-Reissner theory based four node quadrilateral plate bending element using the Integrated Force Method.