Exponential Integral and Euler-Mascheroni Constant

limx0+(Ei(x)ln(x))=γ\large \lim_{x\to 0^{+}} \left( \text{Ei}(-x) - \ln(x)\right) = \gamma

EDIT: I figured out a solution; see my comment below.

To clarify:

  • Ei(y)=yettdt=yettdt\text{Ei}(y) = \displaystyle\int_{-y}^{\infty} \dfrac{-e^{-t}}{t} dt = \displaystyle\int_{-\infty}^{y} \dfrac{e^t}{t}dt denotes the exponential integral;

  • γ \gamma denotes the Euler-Mascheroni constant, γ0.5772\gamma \approx 0.5772 .

  • Notation: Hn H_n denotes the nthn^\text{th} harmonic number, Hn=1+12+13++1n H_n = 1 + \dfrac12 + \dfrac13 + \cdots + \dfrac1n.

I came across this limit when I was trying to derive the derivative of the factorial function: d(n!)dn=n!(Hnγ)\dfrac{d(n!)}{dn} = n!(H_n - \gamma) for all whole numbers nn. Let's denote the derivative as D(n)D(n) We know that (n+1)!=(n+1)n!(n+1)! = (n+1)n!. Therefore, differentiating gives us:

D(n+1)=(n+1)D(n)+n!D(n+1) = (n+1)D(n) + n!

D(n+1)(n+1)!=(n+1)D(n)(n+1)!+n!(n+1)!=D(n)n!+1n+1\dfrac{D(n+1)}{(n+1)!} = \dfrac{(n+1)D(n)}{(n+1)!} + \dfrac{n!}{(n+1)!} = \dfrac{D(n)}{n!} + \dfrac{1}{n+1}

It follows that D(n)n!=Hn+a\dfrac{D(n)}{n!} = H_n + a for some constant aa. What I can't figure out is how to prove that a=γa = -\gamma. I tried to find the value of D(0)D(0):

D(n)=ddn(0xnexdx)=0d(xnex)dndxD(n) = \dfrac{d}{dn}\left(\int_{0}^{\infty} x^n e^{-x} dx \right) = \int_{0}^{\infty} \dfrac{d(x^n e^{-x})}{dn} dx

D(n)=0xnln(x)exdxD(n) = \int_{0}^{\infty} x^n \ln(x) e^{-x} dx

Therefore, D(0)=0ln(x)exdxD(0) = \int_{0}^{\infty} \ln(x) e^{-x} dx Using integration by parts, this becomes:

D(0)=exln(x)00exxdx D(0) = -e^{-x} \ln(x) \rvert_{0}^{\infty} - \int_{0}^{\infty} \dfrac{-e^{-x}}{x}dx

=limy0+[exln(x)yyexxdx]=\lim_{y\to 0^{+}} \left[ -e^{-x} \ln(x) \rvert_{y}^{\infty} - \int_{y}^{\infty} \dfrac{-e^{-x}}{x}dx \right]

D(0)=limy0+(eyln(y)Ei(y))D(0) = \lim_{y\to 0^{+}} (e^{-y} \ln(y) - \text{Ei}(-y))

We can simplify this a little by observing the following: limy0+(ey1)ln(y)=limy0+ey1ylimy0+(yln(y))=10=0\lim_{y\to 0^{+}} (e^{-y} - 1) \ln(y) = \lim_{y\to 0^{+}} \dfrac{e^{-y} - 1}{-y} \lim_{y\to 0^{+}} (-y\ln(y)) = 1 * 0 = 0

Hence, the equation for D(0)D(0) can simplify to:

D(0)=limy0+(ln(y)Ei(y))D(0) = \lim_{y\to 0^{+}} (\ln(y) - \text{Ei}(-y))

Now this is supposed to equal γ-\gamma, but I cannot figure out how to prove it, and can't seem to find the proof anywhere either.

Note by Ariel Gershon
5 years, 2 months ago

No vote yet
1 vote

  Easy Math Editor

This discussion board is a place to discuss our Daily Challenges and the math and science related to those challenges. Explanations are more than just a solution — they should explain the steps and thinking strategies that you used to obtain the solution. Comments should further the discussion of math and science.

When posting on Brilliant:

  • Use the emojis to react to an explanation, whether you're congratulating a job well done , or just really confused .
  • Ask specific questions about the challenge or the steps in somebody's explanation. Well-posed questions can add a lot to the discussion, but posting "I don't understand!" doesn't help anyone.
  • Try to contribute something new to the discussion, whether it is an extension, generalization or other idea related to the challenge.
  • Stay on topic — we're all here to learn more about math and science, not to hear about your favorite get-rich-quick scheme or current world events.

MarkdownAppears as
*italics* or _italics_ italics
**bold** or __bold__ bold

- bulleted
- list

  • bulleted
  • list

1. numbered
2. list

  1. numbered
  2. list
Note: you must add a full line of space before and after lists for them to show up correctly
paragraph 1

paragraph 2

paragraph 1

paragraph 2

[example link](https://brilliant.org)example link
> This is a quote
This is a quote
    # I indented these lines
    # 4 spaces, and now they show
    # up as a code block.

    print "hello world"
# I indented these lines
# 4 spaces, and now they show
# up as a code block.

print "hello world"
MathAppears as
Remember to wrap math in \( ... \) or \[ ... \] to ensure proper formatting.
2 \times 3 2×3 2 \times 3
2^{34} 234 2^{34}
a_{i-1} ai1 a_{i-1}
\frac{2}{3} 23 \frac{2}{3}
\sqrt{2} 2 \sqrt{2}
\sum_{i=1}^3 i=13 \sum_{i=1}^3
\sin \theta sinθ \sin \theta
\boxed{123} 123 \boxed{123}


Sort by:

Top Newest

By definition, Ei(x)=xettdt\displaystyle \text{Ei}(x) = -\int_{-x}^\infty \dfrac{e^{-t}}t \, dt and et=k=0(1)k+1tkk!\displaystyle e^{-t} = \sum_{k=0}^\infty \dfrac{(-1)^{k+1} t^k }{k!} . Integrating the Taylor series of ett \dfrac {e^{-t}}t , we get

Ei(x)=γ+lnx+k=1xkkk!x0. \text{Ei}(x) = \gamma + \ln| x| + \sum_{k=1}^\infty \dfrac{x^k}{k \cdot k!} \qquad x\ne 0 .

Take the limit and the result follows.

Pi Han Goh - 5 years, 2 months ago

Log in to reply

Here's the thing though - how do you know the first term is γ\gamma ??

Ariel Gershon - 5 years, 2 months ago

Log in to reply

Hey Ariel, I see that you like to solve difficult calculus questions. We got a group on another social platform (online forum) that discusses all kinds of difficult calculus questions as well. Plus, most of the active members on Brilliant are on this online forum as well (even the staffs!). We discuss about everything related to math and science. Plus, it's free! Would you like to join us?

If yes, register here

See you!

Aditya Kumar - 5 years, 2 months ago

Log in to reply

Just joined! Thanks! :)

Ariel Gershon - 5 years, 2 months ago

Log in to reply

I've added you!

Aditya Kumar - 5 years, 2 months ago

Log in to reply

I've figured out the solution after a while. Here I will use the generalized harmonic number: Hx=k=1(1k1x+k)H_x = \sum_{k=1}^{\infty} \left(\dfrac{1}{k} - \dfrac{1}{x+k}\right) Now let's start from the point where D(n)n!=Hn+a\dfrac{D(n)}{n!} = H_n + a. Let's integrate this on both sides from 00 to 11: 01D(n)n!dn=01(Hn+a)dn\int_{0}^{1} \dfrac{D(n)}{n!} dn = \int_{0}^{1} (H_n + a) dn ln(n!)01=a+01k=1(1k1n+k)dn\ln(n!) \large{\rvert_0^1} \normalsize{ = a + \int_{0}^{1} \sum_{k=1}^{\infty} \left(\dfrac{1}{k} - \dfrac{1}{n+k}\right) dn} 0=a+k=101(1k1n+k)dn0 = a + \sum_{k=1}^{\infty} \int_{0}^{1} \left(\dfrac{1}{k} - \dfrac{1}{n+k}\right) dn a=k=1(nkln(n+k))01-a = \sum_{k=1}^{\infty} \left(\dfrac{n}{k} - \ln(n+k)\right) \large{\rvert_0^1} a=k=1(1kln(1+k)+ln(k))-a = \sum_{k=1}^{\infty} \left(\dfrac{1}{k} - \ln(1+k) + \ln(k)\right) a=limm(k=1m1kln(1+m))a = -\lim_{m \to \infty} \left(\sum_{k = 1}^{m} \dfrac{1}{k} - \ln(1+m)\right) a=γa = -\gamma

Ariel Gershon - 5 years, 2 months ago

Log in to reply


Problem Loading...

Note Loading...

Set Loading...