d=2
h=0.3
epsilon=1e-6
rx=0.1

[p_max]=ImprovedFastGaussTransformChooseTruncationNumber(d,h,epsilon,rx)
d =

     2


h =

    0.3000


epsilon =

  1.0000e-006


rx =

    0.1000


p_max =

           9