Single-level Anderson quantum dot attached to BCS superconducting leads exhibits a 0 − Pi impurity quantum phase transition, which can be experimentally controlled either by the gate voltage or by the superconducting phase difference. We present two simple analytical formulae describing the position of the phase boundary in parameter space for the weakly correlated and Kondo regime, respectively. Furthermore, we show that the two-level approximation provides an excellent description of the low-temperature physics of superconducting quantum dots near the phase transition. We discuss reliability and mutual agreement of available finite temperature numerical methods (Numerical Renormalization Group and Quantum Monte Carlo) and suggest a novel approach for determination of the quantum phase boundary from measured finite temperature data.