Unlike the previous section, it is not easy to proceed forward with analytical calculations. Hence numerical calculations are done from henceforth. BSMPT,[1]}, [2]}, is a C++ tool used for calculating the strength of the electroweak phase transition in extended Higgs sectors. This is based on the loop-corrected effective potential, also including daisy re-summation of the bosonic masses. The program calculates the vacuum expectation values of the potential (VEV) (\(v\) ) as a function of temperature and, in particular, its value at the critical temperature \(T_c\) . The models implemented within this tool include CP-conserving 2HDM and next-to-minimal 2HDM. In this work, we restrict ourselves to the real sector of CP-conserving 2HDM potential. This tool was used to calculate the VEVs, \(T_c\) s, and the effective value of potential \(V_{eff}(T)\) s for each benchmark point in Table REF . The parameters were chosen to satisfy the limiting conditions for type-I real 2HDM. We assumed that \(VEV/Tc > 0.01\) . It is to be noted that this boundary condition can be anything. The exact value of \(0.01\)  was chosen to broaden the parameter space for scanning. The differential equation in Equation (REF ) was solved numerically for all the benchmark points in Table REF  and the dilution factor as well as the influx of entropy influx for some benchmark points are shown in figure REF .
<TABLE><FIGURE>