The state-of-art software for simulating classical Josephson junction circuits represents their dynamics using a system of differential-algebraic equations (DAEs).