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
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