I need to calculate the volume bounded by: $$x^2+y^2+z^2\:=\:1,\:x^2+\left(y-1\right)^2+z^2=1$$
My solution: Because the volume I want to calculate is symetric, Ill calculate only one half of it and then double it.
$$\frac{1}{2}\le y\le 1$$ $$x^2+z^2\le 1-y^2$$ Then I do the known substitution: $$x\:=\:r\cos \left(\theta \right)\cos \left(\alpha \right),\:y\:=\:r\sin \left(\theta \:\right)\cos \left(\alpha \:\right),\:z\:=\:r\sin \left(\alpha \right)$$
where the new integral is: $$\int \int \int r^2\sin \left(\alpha \right)$$
when: $$-\frac{\pi }{2}\le \theta \le \frac{\pi }{2},\:0\le \alpha \le \pi ,\:1\le r\le 2$$
A little explanation about how i get $r$:$$\frac{1}{2}\le r\sin \left(\theta \right)\cos \left(\alpha \right)\le \frac{r}{2}\le 1\:$$ because: $$\sin \left(\theta \right)\cos \left(\alpha \right)$$ has maximum when: $$\theta =\frac{\pi }{4},\:\alpha =\frac{\pi }{4}$$
Method1 So you have two spheres, both of which are positioned at $(0,0,0)$ and $(0,1,0)$ with radii $r=1$.
So we can use the formulae $x^2+y^2+z^2=1$ and $x^2+(y-1)^2+z^2=1$ to find their intersection (which you've done already:)
Solving simultaneously gives $$y^2=(y-1)^2 \Longrightarrow y=\frac{1}{2}$$
Now we can use triple integrals to find the volume, but I'm gonna just use the Formula $V=\frac{4}{3} \pi r^3$ : for each sphere this makes their volume $\frac{4\pi}{3}$.
But obviously we have to deal with the overlapping region. We employ the trick that the overlapping spheres are axisymmetric. This lets us calculate the overlapping region by integrating Solids of revolution about the y axis; then it will just be a matter of subtracting this overlap region twice to calculate the real volume contained by the overlapping spheres.
What's the curve that we need to revolve? Let's look at the sphere centred at the origin.
The blue curve is just the projection of the sphere into the xy or yz planes. Let's assume its the xy plane. This curve is described by the equation $$x^2+y^2=1$$
So clearly the solid of revolution we need to consider is of the form
$$\int \pi [x(y)]^2 dy$$.
What is our $x(y)$? Well we can rearrange the equation $$x^2+y^2=1 \Longrightarrow x(y)=\sqrt{1-y^2}$$ to get $x$ in terms of $y$.
What are the limits of integration? This is clearly when the intersection starts and the intersection ends with respect to $y$: Well the overlap for the blue region happens between $\frac{1}{2} \leq y \leq 1$. So we have to evaluate
$$\int_{\frac{1}{2}}^{1} \pi [x(y)]^2 dy=\int_{\frac{1}{2}}^{1} \pi (1-y^2)dy$$.
(if my mental arithmetic serves me well!) this integral gives us a volume of $\frac{5 \pi}{24}$.
By symmetry this is half of the duplicated region, so this means the lens shaped region has volume $\frac{5 \pi}{12}$.
Then it's an adding game: The Total region is $$\frac{4\pi}{3}+\frac{4\pi}{3}-\frac{5 \pi}{12}=\frac{27\pi}{12}=\frac{9\pi}{4}$$
For Interest, I plugged this into mathematica to check: