In theory, the five-axis numerical control (NC) machining of sculptured surfaces can be classified into facing milling and cylindrical milling (or side milling). In general, the first one, using flat-end cutter, is suitable for the machining of large sculptured surfaces, e.g. the blade of hydraulic turbine, whose binding relations with drive surface (DS) and check surface (CS) are simple, and the second one, using cylindrical cutter, has wide applications for the milling of small and middle dimensional surfaces whose binding relations with DS and CS are more complex, such as the milling of integral turbine wheels. In practice, the second one suffers more difficulties than the first one, which are mostly related to gouge avoidance, interference avoidance and tool strength. This paper, on the basis of the theories of differential geometry and analytical geometry, describes research on algorithms for the toolpath generation of five-axis cylindrical milling of sculptured surfaces with cylindrical cutter. The approach includes (a) single point offset (SPO) algorithm, and (b) double point offset (DPO) algorithm for the cutter location data (CLDATA) calculation of five-axis cylindrical milling.