The displacements in a laminated composite are represented as products of two sets of unknown functions, one of which is only a function of the thickness coordinate and the other is a function of the in-plane coordinates (i.e., separation of variables approach), and the minimization of the total potential energy is reduced to a sequence of iterative linear problems. Analytical solutions are developed for cross-ply and angle-ply laminated composite rectangular plates. The solution for simply-supported cross-ply plates under sinusoidal transverse load reduces to that of Pagano. Numerical results for stresses and displacements for antisymmetric angle-ply laminates are presented. The three-dimensional elasticity solutions developed are important because they can be used to study the behavior of composite laminates, in addition to serving as reference for approximate solutions by numerical methods and twodimensional theories.