$A = r \int_{0}^{\pi}\int_{0}^{2\pi} e^{i (\alpha+\theta)} d\alpha d\theta = r \int_{0}^{\pi} [\frac{-i}{\alpha+\theta} e^{i(\alpha + \theta)}]_{0}^{2\pi} d\theta = ... = r \int_{0}^{\pi} e^{2i\alpha} ( \frac{1}{\alpha} - \frac{1}{\alpha+2\pi} ) d\theta$
and according to wolfram it is wrong, here. So how do you deduce the area of a sphere in polar coordinates?
[Edit] Trevor gave us formula $r^2 \int_{2π}^{0}\int^{π}_{0}sin(\phi) d\phi d\theta$ without deduction, by deduction, I mean things on which it is based on.
- Why is it odd function $sin(\phi)$ and not instead an even function $cos(\phi)$? (I know that it is and I can calculate it but some proof here missing)
- Why no term with $\theta$ but with $\phi$?
- Why $r^{2}$? -Okay it is area and by symmetry it sounds logical but why $r^{2}$ and not for example $rd$ where $d$ is some lenght or even a path?
- What about an arbitrary case?
- What determines the order in integration?
Many things open to really deduce it.
[Edit] Zhen Lin has an excellent answer where s/he noted that 2-sphere, i.e. the sphere in three dimensional euclidean space, is two-dimensional. Next deductions should be on 0-sphere and then incrementally on others. More about n-sphere here.
- How do you deduce area in
0-sphere,1-sphere,2-sphere, andn-sphere? - Why does the Jacobian matrix in
2-sphere-in-3-dim-euclidean-casehave only two two variables $(\phi, \theta)$? - Are $\theta$ and $\phi$ non-free-variables because they are bounded somehow by the manifold (maybe misusing terminology)?
- What is $r$ in the sphere? For example in the
2-sphere? - I can sense that the var $r$ is somehow different to $\phi$ and $\theta$ because what is it if you derivate it, zero so perhaps it is a trivial case (again probably misusing terminology, sorry)?
- What about if $r$ depends on $t$, time, what is it then? $\frac{\partial r}{\partial t}$ is now something. But would the extra dimension be now $t$ instead of $r$?
Please, use the term dimension in your replies, I feel it is very important to even understand what a sphere is, let alone a dimension.
[Conclusion] The general solution is still open. The $r$ is a homothetic transformation. The terms $(\theta, \phi)$ are apparently tied to the manifold, (ALERT! perhaps misusing terminology). The acceptance of a question does not mean that all things appeared here are solved. Zhen`s solution is really the solution to look at, big thanks for it!

The easiest proof is to use the formula $$A = \int_S \sqrt{\det g} \, du \, dv $$ from Riemannian geometry, where $g = \begin{pmatrix} E & F \\ F & G \end{pmatrix}$ is the metric tensor / first fundamental form and $(u, v)$ are coordinates for the manifold $S$. Unfortunately this doesn't really explain anything, because now the question is how to justify this formula!
So I'll attempt to justify it now. First, we need to restrict to the case where we already know how to calculate areas, namely subsets of the Euclidean plane. We know that the area of a subset $S \subset \mathbb{R}^2$ is just $$A = \int_S dx \, dy$$ where $(x, y)$ are the standard Cartesian coordinates. So, for example, the area of a square with corners $(0, 0), (a, 0), (a, b), (0, b)$ is just $$A = \int_0^b \int_0^a dx \, dy = ab$$ But what about shapes which we can't easily parametrise using Cartesian coordinates? The change of variables formula from multivariable calculus tells us we need to introduce a factor — the Jacobian determinant — when we do this. For example, if we use plane polar coordinates $(r, \theta)$, then $$\begin{aligned} x & = r \cos \theta \\ y & = r \sin \theta \end{aligned}$$ so the Jacobian matrix is $$J = \begin{pmatrix} \cos \theta & -r \sin \theta \\ \sin \theta & r \cos \theta \end{pmatrix}$$ and the Jacobian determinant is $\det J = r$. Hence, the polar-coordinate form of the general formula is $$A = \int_S r \, dr \, d\theta$$ and we can use this to calculate, say, the area of a circle of radius $a$: $$A = \int_0^{2\pi} \int_{0}^{a} r \, dr \, d\theta = \pi a^2$$
But what does any of this have to do with $g$? Well, it turns out that $g = J^T J$, so for example $$g = \begin{pmatrix} 1 & 0 \\ 0 & r^2 \end{pmatrix}$$ for polar coordinates. Now, by the various properties of the determinant, it follows that in general $\det g = (\det J)^2$, i.e. $\det J = \sqrt{\det g}$. So this justifies the formula at the very beginning of my post.
Now let's apply this formula to the sphere. We have spherical polar coordinates $(\phi, \theta)$ such that $$\begin{aligned} x & = r \cos \theta \sin \phi \\ y & = r \sin \theta \sin \phi \\ z & = r \cos \phi \end{aligned}$$ and this gives the Jacobian $$J = \begin{pmatrix} r \cos \theta \cos \phi & -r \sin \theta \sin \phi \\ r \sin \theta \cos \phi & r \cos \theta \sin \phi \\ -r \sin \phi & 0 \end{pmatrix}$$ so the metric is $$g = J^T J = \begin{pmatrix} r^2 & 0 \\ 0 & (r \sin \phi)^2 \end{pmatrix}$$ hence the area of the sphere is just $$A = \int_0^{2\pi} \int_0^\pi r^2 \sin \phi \, d\phi \, d\theta = 4 \pi r^2$$ as we all know.