Pre-main-sequence evolutionary tracks of massive stars have been calculated under the canonical theory and the mass accretion paradigm. The canonical tracks for 15, 20, 25, 30, 45, and 60 M. stars have been calculated with both Cox-Stewart (1970) and the Rogers and Iglesias (1992a, b) opacities. Stellar models less massive than about 30 M. develop convective cores only during the onset of central hydrogen burning. Because of the increasing importance of radiation pressure, stars more massive than 30 M. develop convective cores before the onset of hydrogen burning. A generalized Naur-Osterbrock ( 1953) criterion has been used to understand the development of the convective core. It is estimated that stars more massive than approximate to 95 M. evolving canonically to the main sequence will have a convective core throughout the whole pre-main-sequence phase. The canonical 30 M. model has been evolved through the main-sequence phase using the mass-loss algorithm of de Jager, Nieuwenhuijzen, and van der Hucht( 1988). The pre-main-sequence evolution has also been calculated under the accretion paradigm: beginning with a collapsed core, mass was accreted until a final stellar mass was reached. The initial hydrodynamic evolution of the protostar was not followed. The calculation was begun with a 1 M. core, and a constant accretion rate of 10(-5) M. yr(-1) was used. The initial model is fully convective and burning deuterium in the core. The accretion track intersected the canonical main sequence at M = 8.5 M.. For masses from 8.5 to approximate to 15 M., the accretion and canonical zero-age main-sequence models were essentially the same. Beyond approximate to 17 M., the accretion models became increasingly more luminous and cooler with smaller convective cores. The main-sequence evolution of the 30 M. star, begun with an accretion model, occurs at a slightly lower luminosity and a lower effective temperature compared with canonical evolution, The accretion sequence was calculated up to the final mass of 45 M.. The locus where massive stars first become optically visible, the upper stellar birth line, was determined from these models. It is pointed out that the upper mass limit of stable stars must be determined by the formation environment of the star rather than the onset of nuclearly energized pulsation instability (Stothers 1992), since the increasingly large chemical inhomogeneity that develops as a massive star accretes matter evolving along the upper stellar birth line will stabilize nuclearly driven pulsations.