added code to determine upper limit on rbf sum

When we first wrote this RBF code, the upper limit in the sum was hard-wired at n_max = 40. I've added code that determines an upper limit for the sum, once the shape parameter alpha is specified.