Here follow Orlando et al., PRB 60, 15398 (1999) and numerically calculate the wavefunctions and energy-levels of the Flux qubit device using a two-dimensional potential of the superconducting Flux qubit circuit. The two generalized coordinates are the gauge invariant phases across the Josephson junctions in the device.