Numerical evaluation of the response of a Josephson tunnel junction in an arbitrary circuit

A numerical technique is presented for calculating the response of a Josephson tunnel junction when connected to an arbitrary circuit. The approach includes all of the details of the microscopic theory including sine and cosine terms as well as the quasiparticle current together with its reactive art. An example is given, and techniques for reducing computation time are discussed.