  In calculating one-loop contributions to amplitudes of the lepton-flavor
violating decays of the neutral Higgses (LFVHD) to different flavor charged
leptons, the analytic expressions can be written in term of Passarino-Veltman
functions. Then, they can be computed numerically by Looptools [1]. Another
approach is using suitable analytic expressions established for just this
particular case. We compare numerical results obtained from Looptools and those
computed by different expressions that have been applied recently. Then we
derive the preferable ones that are applicable for the large ranges of free
parameters introduced in the models beyond the Standard Model (SM). For
illustration, the LFVHD in a simple model, which has been discussed recently,
will be investigated more precisely.
