One-loop contributions to the amplitudes are evaluated in terms of the class of Passarino-Veltman functions (PV functions) generally dened in [1]. Numerical calculations related to these functions can be done directly with either LoopTools package or more popular ones such as Mathematica implementing analytic expressions found by hand (in possible cases). Through numerical investigations to compare analytical formulas built and coded by us with the equivalent ones dened in LoopTools, we have found the very high consistence between them. Especially if masses in the loop exceed the value about 108 GeV, some PV functions of Loop- Tools appear unstable while ours still work well.