d=2
h=0.3
epsilon=1e-6
Klimit=round(40*sqrt(d)/h)

[K,p_max,r]=ImprovedFastGaussTransformChooseParameters(d,h,epsilon,Klimit)
d =

     2


h =

    0.3000


epsilon =

  1.0000e-006


Klimit =

   189


K =

          13


p_max =

          17


r =

    1.1151