An efficient, robust, ODE-based formulation of the continuation problem is presented for the automatic construction of all branches of the pinch point curves in any homogeneous distillation system (even highly nonideal ones) that can be represented using an activity coefficient model, with any number of components and for any adiabatic, nonreactive column configuration. Constant molar overflow is assumed. The curve is parametrized by reflux ratio, making it particularly applicable to distillation design and synthesis. Naturally, it can be used to find pinch point locations at a given finite reflux, as this is simply a subset of the pinch point curve. The method can also automatically locate all N-component azeotropes very efficiently (orders of magnitude faster than commercial process simulators). It has the additional advantage of being solvable using tools already available in typical mathematical software packages. A computer implementation of the proposed method is available.