A Dynamic Relaxation (DR) program based on finite differences has been developed for small and large deflection analysis of rectangular laminated plates using first order shear deformation theory (FSDT). The displacements are assumed linear through the thickness of the plate. Dynamic Relaxation (DR) method is presented for the geometrically linear and nonlinear laterally loaded, rectangular laminated plates. The analysis uses the Mindlin plate theory with first order shear deformation theory (FSDT) which accounts for transverse shear deformation. A computer program has been compiled using a FORTRAN program. The convergence and accuracy of the DR solutions for elastic small and large deflection response are established by comparison with various exact and approximate solutions. New numerical results are generated for uniformly loaded square laminated plates which serve to quantify the effects of shear deformation, material anisotropy, fiber orientation, and coupling between bending and stretching. It was found that linear analysis seriously over-predicts deflections of plates. The shear deflection depends greatly on a number of factors such as length/ thickness ratio.