A new algorithm and software for numerical sensitivity analysis of differential-algebraic equations are presented. The staggered corrector method proposed has lower computational complexity than both the staggered direct and the simultaneous corrector methods. The results here are especially significant for, though not limited to, large-scale sparse differential-algebraic systems.