In this work we present a new computation of the lepton flavor violating Higgs boson decays that are generated radiatively to one-loop from heavy right-handed neutrinos. We work within the context of the inverse seesaw model with three $��_R$ and three extra singlets $X$, but the results could be generalized to other low scale seesaw models. The novelty of our computation is that it uses a completely different method by means of the mass insertion approximation which works with the electroweak interaction states instead of the usual 9 physical neutrino mass eigenstates of the inverse seesaw model. This method also allows us to write the analytical results explicitly in terms of the most relevant model parameters, that are the neutrino Yukawa coupling matrix $Y_��$ and the right-handed mass matrix $M_R$, which is very convenient for a phenomenological analysis. This $Y_��$ matrix, being generically nondiagonal in flavor space, is the only responsible for the induced charged lepton flavor violating processes of our interest. We perform the calculation of the decay amplitude up to order ${\cal O}(Y_��^2+Y_��^4)$. We also study numerically the goodness of the mass insertion approximation results. In the last part we present the computation of the relevant one-loop effective vertex $H\ell_i\ell_j$ for the lepton flavor violating Higgs decay which is derived from a large $M_R$ mass expansion of the form factors. We believe that our simple formula found for this effective vertex can be of interest for other researchers who wish to estimate the $H \to \ell_i \bar \ell_j$ rates in a fast way in terms of their own preferred input values for the relevant model parameters $Y_��$ and $M_R$.

36 pages, 15 figures. This version v3 matches the manuscript published in PRD