For a general class of nonlinear (possibly higher index) differential-algebraic equations we show existence and uniqueness of solutions. These solutions are regular in the sense that Newton's method will converge locally and quadratically. On the basis of the presented theoretical results, numerical methods for the determination of consistent initial values and for the computation of regular solutions are developed. Several numerical examples are included. ma@