A model of thermoelasticity in laminated thin films is derived from three-dimensional considerations. Its modification to accommodate wrinkling and slackening is described and a procedure for numerical solution is developed. This is used to solve a number of equilibrium problems that exhibit the effects of heating and mechanical loading on the stress field and consequent wrinkling pattern.