A calculation method of the breakdown voltage for the drain region with the cylindrical structure in LDMOS is proposed. The depletion region of the drain is divided into many smaller regions and the doping concentration of each split region is assumed to be uniformly distributed. The field and potential in each split region is calculated by the integration of the Poisson equation and the ionization integral method is used to compute the breakdown voltage. The breakdown voltage resulted from the proposed method shows the maximum relative error of 2.2% compared with the result of the 2-dimensional device simulation using BANDIS.