For the positive solution of the higher order differential-algebra equations, predictor-corrector numerical integration methods with error feedback control were proposed.