Here, we present a method of simulating Josephson junction circuits based on graph theory that eliminates all algebraic equations to create a system of only ordinary differential equations (ODEs).