# Proving $\int_0^r{(r^m-x^m)^{1/m}dx}=\frac{\Gamma\left(\frac{1}{m}+1\right)\Gamma\left(\frac{1}{m}+1\right)}{\Gamma\left(\frac{2}{m}+1\right)}r^2$

First let's put the question succinctly. How can I go about showing the following?

$$\int_0^r{(r^m-x^m)^{1/m}dx}=\frac{\Gamma\left(\frac{1}{m}+1\right)\Gamma\left(\frac{1}{m}+1\right)}{\Gamma\left(\frac{2}{m}+1\right)}r^2$$

Now for some exposition: I am a math enthusiast and this result kind of fell into my lap after playing around a bit with "circles"... This is my first encounter with the $\Gamma$ function. I am not quite sure how one goes about establishing such a claim. At this point I am at the "I better look into this $\Gamma$ function" part of my research but I figured I would document the question and take any input offered.

Consider the equation $|x|^m+|y|^m=1$, for $m \in {1,2,3}$

The $m=1$ case then corresponds to the square in the picture which has side lengths $\sqrt{2}$. The whole square has area $2$ and therefore the area of the square limited to the first quadrant is $1/2$.

$$\int_0^1{(1-x)dx}=\frac{\Gamma(2)\Gamma(2)}{\Gamma(3)}=\frac{(2-1)!(2-1)!}{(3-1)!}=\frac{1}{2}$$

I only invoke the idea that over the whole numbers $\Gamma(n+1)=n!$ here because I found the formula by examining this in the case when my inputs for $\Gamma$ were whole numbers. Then I replaced my factorial symbols with $\Gamma$s to get the claim above which I have only verified empirically.

For the $m=2$ case. We have the unit circle. The area in the first quadrant should be $\pi/4$. And indeed:

$$\int_0^1{(1-x^2)^{1/2}dx} =\frac{\Gamma\left(\frac{3}{2}\right)\Gamma\left(\frac{3}{2}\right)}{\Gamma(2)} =\frac{ \sqrt{\pi}}{2}\frac{ \sqrt{ \pi} }{2}=\dfrac\pi4$$

Cool! So now I was excited to see that this worked not only in the cases with whole number inputs to $\Gamma$.

$m=3$ Well then what's the area under the curve $|x|^3+|y|^3=1$? This corresponds to the outermost curve in the diagram. Well... I assume this value must be some transcendental number. It's construction is similar to the way we think about $\pi$. But what is it?

\begin{align*}\int_0^1{(1-x^3)^{1/3}dx}&=\frac{\Gamma(\frac{1}{3}+1)\Gamma(\frac{1}{3}+1)}{\Gamma(\frac{2}{3}+1)}\\ &\approx 0.883319375142724978656844749824219351285934269101278765063\end{align*}

Which matches up with numerical integration. Wolfram alpha can present this number in a few other ways. For example, $$\frac{\Gamma(1/3)^3}{4\sqrt{3}{\pi}}$$ These other representations all seem to invoke the Gamma function.

• Do you mean $\int_0^r (r^m-x^m)^\frac{1}{m} dx$ ? Apr 21, 2018 at 13:16
• Yes. Thank you. Apr 21, 2018 at 13:17
• desmos.com/calculator/u8uvtw7cnc If you want to play with my empirical evidence. Apr 21, 2018 at 13:45
• I checked up to $m=50$ for 50 decimal places. It works ! Beautiful. This is also the result from a CAS. Apr 21, 2018 at 13:49
• Don't know whether it helps, but I find $\frac{\Gamma(\frac{1}{m})\Gamma(\frac{1}{m})}{2m\Gamma(\frac{2}{m})}$ to look a bit nicer. Apr 21, 2018 at 13:57

## 3 Answers

Under $x\to rx$ and $x^m\to x$, one has \begin{eqnarray} &&\int_0^r{(r^m-x^m)^{1/m}dx}\\ &=&\int_0^1{(r^m-r^mx^m)^{1/m}rdx}\\ &=&r^2\int_0^1(1-x^m)^{1/m}dx\\ &=&r^2\frac1m\int_0^1(1-x)^{1/m}x^{\frac1m-1}dx\\ &=&r^2\frac1m\frac{\Gamma(\frac1m+1)\Gamma(\frac1m)}{\Gamma(\frac2m+1)}\\ &=&\frac{\Gamma^2(\frac{1}{m}+1)}{\Gamma(\frac{2}{m}+1)}r^2. \end{eqnarray} Here $$\int_0^1x^{p-1}(1-x)^{q-1}dx=\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+q)},\Gamma(x+1)=x\Gamma(x)$$ are used.

• @Mason, should be $p+q$. Apr 21, 2018 at 17:00

Use Beta function and with substitution $x=r\cos^{\frac2m}t$ write \begin{align} \int_0^r(r^m-x^m)^{1/m}dx &= \dfrac{r^2}{m} 2\int_{0}^{\frac{\pi}{2}}\sin^{\frac2m+1}t\cos^{\frac2m-1}t dt\\ &= \dfrac{r^2}{m} \beta\left(\frac2m+1,\frac2m\right)\\ &= \frac{\Gamma\left(\frac{1}{m}+1\right)\Gamma\left(\frac{1}{m}+1\right)}{\Gamma\left(\frac{2}{m}+1\right)}r^2 \end{align}

After playing a little with xpaul's proof I think I can take on the generalized squircle $|x|^a+|y|^b=1$. Claim: In general the area under this curve is $\frac{4}{a+b}\beta(a^{-1},b^{-1})$.

I will take on the general radius soon. But for now while I am relearning calculus:

$$\int_{0}^{1} ({1-x^a})^{1/b}dx$$

Let $u=x^a$ and note that this means that $(u^{1/a})^{a-1}=x^{a-1}$

$du=ax^{a-1}=au^{1-\frac{1}{a}}dx$ and this means that $dx=\frac{u^{\frac{1}{a}-1}du}{a}$. In this notation

$$\int_{0}^{1} ({1-x^a})^{\frac{1}{b}}dx=\frac{1}{a}\int_{0}^{1} ({1-u})^{\frac{1}{b}}u^{\frac{1}{a}-1}du$$

Using the identity from xpaul's post. We find this is equal to $$\frac{a^{-1}\Gamma(a^{-1})\Gamma(b^{-1}+1)}{\Gamma(a^{-1}+b^{-1}+1)}$$ using the identity $\Gamma(x+1)=x\Gamma(x)$ this simplifies to $$\frac{a^{-1}b^{-1}\Gamma(a^{-1})\Gamma(b^{-1})}{(a^{-1}+b^{-1})\Gamma(a^{-1}+b^{-1})}$$ Multiplying by $\frac{ab}{ab}$ we arrive at $$\frac{\Gamma(a^{-1})\Gamma(b^{-1})}{(a+b)\Gamma(a^{-1}+b^{-1})}$$

Then multiplying by $4$ get the area of the whole squircle. Note that we take $a=b=m$ this simplifies to the solution desired in this post's question.

• This is also given here Feb 2, 2020 at 5:44