We have developed a mathematical model for the plant vascular system that enables the prediction of a compound's phloem systemicity as a function of its partition coefficients and acid dissociation constants. The mathematical model can account for the sensitivity of systemicity to plant parameters such as plant size and pH of the phloem sap. This paper reviews this model and demonstrates how it accounts for the phloem systemic properties of most herbicides as well as that of many endogenous substances such as plant hormones. The model also can be used to design phloem systemic pesticides as illustrated for a pronematicide that successfully controls nematodes when applied foliarly to transgenic tobacco plants capable of regenerating the parent nematicide in a root specific fashion.