  We demonstrate that our model [Phys.Rev. E91, 032312 (2015)] serves as a
useful tool to trace the evolution of equilibrium configurations of
one-component charged particles confined in a disk. Our approach reduces
significantly the computational effort in minimizing the energy of equilibrium
configurations and demonstrates a remarkable agreement with the values provided
by molecular dynamics calculations. We show that the Comment misrepresents our
paper, and fails to provide plausible arguments against the formation hexagonal
structure for n>200 in molecular dynamics calculations.
