One loop contributions to neutral Higgs decay h-->mu tau

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.