A major determinant of plant architecture is the arrangement of branches around the stem, known as phyllotaxis. However, the specific form of branching conditions is not known. Here we discuss this question and suggest a branching model which seems to be in agreement with biological observations.Recently, a number of models connected with the genetic network or molecular biology regulation of the processes of pattern formation appeared. Most of these models consider the plant hormone, auxin, transport and distribution in the apical meristem as the main factors for pattern formation and phyllotaxis. However, all these models do not take into consideration the whole plant morphogenesis, concentrating on the events in the shoot or root apex. On the other hand, other approaches for modeling phyllotaxis, where the whole plant is considered, usually are mostly phenomenological, and due to it, do not describe the details of plant growth and branching mechanism. In this work, we develop a mathematical model and study pattern formation of the whole, though simplified, plant organism where the main physiological factors of plant growth and development are taken into consideration. We model a growing plant as a system of intervals, which we will consider as branches. We assume that the number and location of the branches are not given a priori, but appear and grow according to certain rules, elucidated by the application of mathematical modeling. Four variables are included in our model: concentrations of the plant hormones auxin and cytokinin, proliferation and growth factor, and nutrients-we observe a wide variety of plant forms and study more specifically the involvement of each variable in the branching process. Analysis of the numerical simulations shows that the process of pattern formation in plants depends on the interaction of all these variables. While concentrations of auxin and cytokinin determine the appearance of a new bud, its growth is determined by the concentrations of nutrients and proliferation factors. Possible mechanisms of apical domination in the frame of our model are discussed.