Compute bounds as a function of one parameter for more stable computations

[buntonj]

Currently, bound relies on dividing powers of $\delta$ with $(d+1)!$ which is numerically unstable. Could try and connect some of the parameters and make the bound computations more numerically stable in practice.