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


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$.


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.

  • $\begingroup$ Do you mean $\int_0^r (r^m-x^m)^\frac{1}{m} dx$ ? $\endgroup$ – Delta-u Apr 21 '18 at 13:16
  • $\begingroup$ Yes. Thank you. $\endgroup$ – Mason Apr 21 '18 at 13:17
  • $\begingroup$ desmos.com/calculator/u8uvtw7cnc If you want to play with my empirical evidence. $\endgroup$ – Mason Apr 21 '18 at 13:45
  • $\begingroup$ I checked up to $m=50$ for 50 decimal places. It works ! Beautiful. This is also the result from a CAS. $\endgroup$ – Claude Leibovici Apr 21 '18 at 13:49
  • 1
    $\begingroup$ 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. $\endgroup$ – Arthur Apr 21 '18 at 13:57

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.

  • 1
    $\begingroup$ @Mason, should be $p+q$. $\endgroup$ – xpaul Apr 21 '18 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.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.