  Based on lattice QCD results of equation of state (EOS) and baryon number
susceptibility at zero baryon chemical potential, and supplemented by machine
learning techniques, we construct the analytic form of the holographic black
hole metric in the Einstein-Maxwell-Dilaton (EMD) framework for pure gluon,
2-flavor, and 2+1-flavor systems, respectively. The dilaton potentials solved
from Einstein equations are in good agreement with the extended non-conformal
DeWolfe-Gubser-Rosen (DGR) type dilaton potentials fixed by lattice QCD EOS,
which indicates the robustness of the EMD framework. The predicted critical
endpoint (CEP) in the 2+1-flavor system is located at $(T^c$=0.094GeV,
$\mu^c_B$=0.74GeV), which is close to the results from the realistic
Polyakov-Nambu-Jona-Lasinio(PNJL) model, the functional renormalization group,
and the holographic model with extended DeWolfe-Gubser-Rosen dilaton potential.
