In this paper we develop a general technique for establishing analyticity of solutions of partial differential equations which depend on a parameter e. The technique is worked out primarily for a free boundary problem describing a model of a stationary tumor. We prove the existence of infinitely many branches of symmetry-breaking solutions which bifurcate from any given radially symmetric steady state; these asymmetric solutions are analytic jointly in the spatial variables and in e. 1. THE MODEL AND MAIN RESULT In this paper we present a general technique for establishing analyticity of solutions of systems of partial differential equations which depend analytically on a parameter e. The method works not only for boundary value problems but also for free boundary problems. In this latter context it can be used to establish long time existence of transient solutions, and also to study the existence of spatially asymmetric steady solutions. Since free boundary problems are typically more challenging than their boundary-value counterparts, we shall concentrate here on a free boundary problem from developmental biology, namely, a model of tumor growth. To further exemplify the generality of our approach an instance of a boundary value problem (in a fixed domain) is presented in the last section of the paper. A variety of other problems are amenable to the same analysis, including, in particular, the Hele-Shaw model of fluid flow [11]. Within the last several decades a number of mathematical models have been developed that aimed at describing the evolution of carcinomas (see. e.g., [1, 5, 6, 8, 12, 13] and the references cited there). The main objective of these models has been to qualitatively describe, under various simplifying assumptions, the growth and stability of tumor tissue. Analysis and simulations of such models are helping to assess the relative importance of various mechanisms affecting tumor growth as well as the efficacy of certain cancer treatments. On the other hand, the description of the stationary (dormant) configurations that arise from the models has only been addressed in the case of spherical tumors, but otherwise it remains largely unexplored. In this paper we develop a method for establishing analyticity of Received by the editors August 17, 1999. 1991 Mathematics Subject Classification. Primary 35B32, 35R35; Secondary 35B30, 35B60, 35C10, 35J85, 35Q80, 92C15, 95C15.