Efficient algorithm for computing large scale systems of differential algebraic equations