This code supports the study of Andreev bound states (ABS) in a hybrid junction composed of a quantum spin Hall (QSH) system with helical edge states, contacted by two superconducting electrodes.The normal region includes an arbitrary ferromagnetic domain (F). The code is structured around two main components: A numerical solver for the ABS energy spectrum as a function of the superconducting phase difference. To use it, the user must specify all the other parameters entering the ABS equation. A post-processing routine that organizes the numerical output to facilitate visualization and plotting, tailored to the specific physical configuration analyzed.