Gap equation and GL expansion coefficient calculation code and the data of figures

The code was compiled and debugged using the Wolfram Mathematica software. The code generation process proceeded as follows: first, the complete forms of the gap equation and the Ginzburg–Landau (GL) expansion coefficients appearing in the section “Explicit evaluations and regularizations” were implemented in the code.

After substituting the relevant model parameters, the code was compiled and executed, yielding the numerical results for the gap equation and the GL expansion coefficients, which serve as the data required for plotting. Subsequently, these numerical data were used to generate the figures shown in Section III “NUMERICAL RESULTS”, and the plotting routines were optimized so that the figures could be directly inserted into the paper.