Suppose we have a ($n-1$ dimensional) unit sphere centered at the origin: $$ \sum_{i=1}^{n}{x_i}^2 = 1$$ Given some some $d \in [0,1]$, what is the probability that a randomly selected point on the sphere, $ (x_1,x_2,x_3,...,x_n)$, has coordinates such that $$|x_i| \leq d$$ for all $i$?

This is equivalent to finding the intersection of the $(n-1)$-hypersphere with the $n$-hypercube of side $2d$ centered at origin, and then taking the ratio of that $(n-1)$-volume over the $(n-1)$-volume of the $(n-1)$-hypersphere.

As there are closed-form formulas for the volume of a hypersphere, the problem reduces to finding the $(n-1)$-volume of the aforementioned intersection.

All attempts I've made to solve this intersection problem have led me to a series of nested integrals, where one or both limits of each integral depend on the coordinate outside that integral, and I know of no way to evaluate it. For example, using hyperspherical coordinates, I have obtained the following integral: $$ 2^n n! \int_{\phi_{n-1}=\tan^{-1}\frac{\sqrt{1-(n-1)d^2}}{d}}^{\tan^{-1}1} \int_{\phi_{n-2}=\tan^{-1}\frac{\sqrt{1-(n-2)d^2}}{d}}^{\tan^{-1}\frac{1}{\cos\phi_{n-1}}}\ldots\int_{\phi_1=\tan^{-1}\frac{\sqrt{1-d^2}}{d}}^{\tan^{-1}\frac{1}{\cos\phi_2}} d_{S^{n-1}}V$$ where $$d_{S^{n-1}}V = \sin^{n-2}(\phi_1)\sin^{n-3}(\phi_2)\cdots \sin(\phi_{n-2})\ d\phi_1 \ d\phi_2\ldots d\phi_{n-1}$$ is the volume element of the $(n-1)$–sphere. But this is pretty useless as I can see no way of integrating this accurately for high dimensions (in the thousands, say).

Using cartesian coordinates, the problem can be restated as evaluating: $$ \int_{\sum_{i=1}^{n-1}{x_i}^2\leq1, |x_i| \leq d} \frac{1}{\sqrt{1-\sum_{i=1}^{n-1}{x_i}^2}}dx_1 dx_2 \ldots dx_{n-1}$$ which, as far as I know, is un-integrable.

I would greatly appreciate any attempt at estimating this probability (giving an upper bound, say) and how it depends on $n$ and $d$. Or, given a particular probability and fixed $d$, to find $n$ which satisfies that probability.

Edit: This question leads to two questions that are slightly more general:

I think part of the difficulty is that neither spherical nor Cartesian coordinates work very well for this problem, because we're trying to find the intersection between a region that is best expressed in spherical coordinates (the sphere) and another that is best expressed in Cartesian coordinates (the cube). Are there other problems that are similar to this? And how are their solutions usually formulated?

Also, the problem with the integral is that the limits of each of the nested integrals is a function of the "outer" variable. Is there any general method of solving these kinds of integrals?

2more comments