Example of a problem with two superconducting leads and some phase
difference
Rok Zitko, Aug 2011

The model is defined as a Mathematica script siamphase.m

Caveat: in F.dat, only imaginary part of the delta peak weights is
taken into account. One should also compute the real parts
(reim=true), do a Kramers-Kronig transformation of the real part, and
add that to the anomalous spectral function.
