I. INTRODUCTION
For to calculate the transient electromagnetic responses of
horizontal layered media which is excited by a large circular
loop antenna by the transformation method , the double-order
double-argument Hankel transformation which is defined in (1)
are encountered[1], the efficiency of its numerical algorithm is
very important because numerous of (1) must to be computed
in frequency domain[1-3],