Numerical Methods for Stiff Index-3 DAEs