A random sphere containing the center of the unit cube Inspired by a Putnam problem, I came up with the following question:
A point in randomly chosen in the unit cube, a sphere is then created using the random point as the center such that the sphere must be contained inside the cube (In other words, the largest sphere that fits). What is the probability that the center of the cube is contained inside the sphere created?
No real idea how to approach this one but thought some of you might find this interesting.
 A: I'll try an approach using order statistics. Suppose the coordinates of the chosen random point are $(X_1, X_2, X_3)$. We can assume that the coordinates of the chosen point are all positive. (If they are not, we can reflect the point into the first octant.) Hence, $X_1, X_2, X_3$ are independent uniform random variables on $(0, 1)$.
Now, we define the order statistics $Y_1, Y_2, Y_3$ so that $Y_1$ is the smallest of the $X_i$ values, $Y_2$ is the middle value, and $Y_3$ is the largest value. Note that the $Y_i$ variables are neither uniform on $(0, 1)$ nor independent of one another.
There are two variables of interest: $R$, the radius of the sphere, and $D$, the distance from the chosen point to the origin. Note that $R = \min\{1 - X_1, 1 - X_2, 1 - X_3\} = 1 - Y_3 $ and that $D = \sqrt{X_1^2  + X_2^2 + X_3^2} = \sqrt{Y_1^2 + Y_2^2 + Y_3^2}$. The operative question is: what is $\mathbb P(D < R)$? That is, what is $$\mathbb P \left(\sqrt{Y_1^2 + Y_2^2 + Y_3^2} < 1 - Y_3 \right)?$$
First, a bit of algebra to clean this up:
\begin{align*}
  \mathbb P \left(\sqrt{Y_1^2 + Y_2^2 + Y_3^2} < 1 - Y_3 \right) &= \mathbb P \left( Y_1^2 + Y_2^2 + Y_3^2 < (1 - Y_3)^2 \right) \\
&= \mathbb P \left(Y_1^2 + Y_2^2<1-2 Y_3 \right)
\end{align*}
We know the joint pdf of these order statistics to be
$$f(y_1, y_2, y_3) = \begin{cases}
  3!, & 0 < y_1 < y_2 < y_3 < 1 \\
0, & \text{otherwise} \end{cases}$$
so we just need to integrate that density over the set $\{y_1^2 + y_2^2 < 1 - 2 y_3\}$ in the cube $[0, 1]^3$. Note that this requires in particular that $y_3 \leq 1/2$. I claim that this triple integral can be expressed as
\begin{align*}
  \int_0^{1/2} \int_0^{\min\{y_3, \sqrt{1 - 2 y_3}\}} \int_0^{\min\{y_2, \sqrt{1 - 2y_3 - y_2^2}\}} 6 \, \textrm d y_1 \, \textrm d y_2 \, \textrm d y_3.
\end{align*}
I don't know the analytic value of that integral (I haven't tried very hard yet), but Wolfram Alpha estimates it to be $\fbox{0.0868}$.

When I do a really long probability computation like this, I always assume I screwed up somewhere and verify my work with a Monte Carlo simulation. Here's that work in R:
spherecube <- function(){
  center <- runif(3, min=-1, max=1)
  radius <- min(abs(1 - center), abs(1 + center))
  sum(center * center) < radius^2
}

mean(replicate(100000, spherecube()))
# 0.08674

Aside from being an interesting problem, it's a great advertisement for the power of Monte Carlo simulations!
A: Here's the solution, confirming the previous answers, using symmetry as I suggested. I'll reorder my variables to conform with previous discussions. So we're going to consider the cube $[-1,1]^3$ and restrict to centers of the sphere in the pyramid $0\le x\le y\le z\le 1$. This means that the face $z=1$ will be the closest. The origin will be (on or) inside such a sphere if and only if $x^2+y^2+z^2\le (1-z)^2$, i.e., $2z\le 1-(x^2+y^2)$. 
Over what region in the $xy$-plane does our region project? Since $x\le y\le z$, we must have $2y\le 2z\le 1-(x^2+y^2)$, which means $x^2+y^2+2y\le 1$, or $x^2+(y+1)^2\le 2$. This results in the portion of $0\le x\le y$ lying inside the circle $x^2+(y+1)^2\le 2$. Note that $0\le x\le \dfrac{\sqrt3-1}2$.
Setting up the triple integral, the volume we desire is
$$\int_0^{\frac{\sqrt3-1}2}\int_x^{\sqrt{2-x^2}-1}\int_y^{\frac12(1-x^2-y^2)} dz\,dy\,dx,$$
and since we're comparing to the volume of the full pyramid, which is $1/6$, we take $6$ times this answer. 
(We can also set this up nicely in polar coordinates:
$$\int_0^{\pi/4}\int_0^{\sqrt{\cos^2\theta+1}-\cos\theta}\int_{r\sin\theta}^{\frac12(1-r^2)}\,r\,dz\,dr\,d\theta.)$$
The integral, multiplied by $6$, turns into
\begin{align*}
\int_0&^{\frac{\sqrt3-1}2} \big({-}5+4\sqrt{2-x^2}+3x^2-2x^2\sqrt{2-x^2}-(3x-3x^2-4x^3)\big)dx\\ &= \int_0^{\frac{\sqrt3-1}2} \big({-}5-3x+4x^3+4\sqrt{2-x^2}+6x^2-2x^2\sqrt{2-x^2}\big)dx \\ &= \frac54-\frac98\sqrt3+\frac{\pi}4.
\end{align*}
