Tunable difference equations for the time-domain simulation of power system operating states