Piecewise quasilinearization methods for singular boundary-value problems in second-order ordinary differential equations are presented. These methods result in linear constant-coefficients ordinary differential equations which can be integrated analytically, thus yielding piecewise analytical solutions. The accuracy of the globally smooth piecewise quasilinear method is assessed by comparisons with exact solutions of several Lane-Emden equations, a singular problem of non-Newtonian fluid dynamics and the Thomas-Fermi equation. It is shown that the smooth piecewise quasilinearization method provides accurate solutions even near the singularity and is more precise than (iterative) second-order accurate finite difference discretizations. It is also shown that the accuracy of the smooth piecewise quasilinear method depends on the kind of singularity, nonlinearity and inhomogeneities of singular ordinary differential equations. For the Thomas-Fermi equation, it is shown that the piecewise quasilinearization method that provides globally smooth solutions is more accurate than that which only insures global continuity, and more accurate than global quasilinearization techniques which do not employ local linearization.