The generalised Dirichlet integral: integral of (sinx)^n/x^n from zero to infinity

Поділитися
Вставка
  • Опубліковано 27 лис 2024

КОМЕНТАРІ • 124

  • @rajendramisir3530
    @rajendramisir3530 Рік тому +27

    I am impressed by your generalized solution of this integral. This is the first generalized solution of this integral I have seen so far. That was a good exercise for you. Thanks for sharing.

  • @holyshit922
    @holyshit922 Рік тому +115

    You could cancel Gamma function with factorial

    • @GreenMeansGOF
      @GreenMeansGOF Рік тому +12

      Yeah. n!/Γ(n)=n

    • @gumkiller9734
      @gumkiller9734 Рік тому +10

      Ong

    • @ShanBojack
      @ShanBojack Рік тому +2

      Can you please elaborate

    • @pinsonraphael4873
      @pinsonraphael4873 Рік тому +6

      @@ShanBojack For any integer n>1, gamma(n) = (n-1)! It's pretty easy to prove by induction

    • @tsa_gamer007
      @tsa_gamer007 Рік тому +4

      Don't spoil the video
      You spoiled the video by saying that the video features gamma function😂

  • @ignaciorodriguez639
    @ignaciorodriguez639 Рік тому +12

    Instead of integrating by parts, you can expand ( sin x ) ^ n as a sum of complex exponentials. This way, you can find a closed formula
    integral ( ( sin ( x ) / x ) ^ n , x = 0 , x = infinity ) = n * pi * sum ( ( -1 ) ^ k * ( n - 2 * k ) ^ ( n - 1 ) / ( ( n - k ) ! k ! ) , k = 0 , k = floor ( ( n - 1 ) / 2 ) ) / ( 2 ^ n )

    • @Anonymous-Indian..2003
      @Anonymous-Indian..2003 Рік тому +1

      Bruh....
      I did in a same way
      I generated the exact same to same formula in November 2022...........

  • @homerthompson416
    @homerthompson416 Рік тому +10

    Wow never in my life would I have guessed you could reduce the integrals of (sin(x)/x)^n to a sum of simple integrals of the form ∫ dt/(a^2+t^2)

  • @vladimirlucic1276
    @vladimirlucic1276 Рік тому +12

    This was feature as Problem 1064 in Mathematics Magazine (in 1979). Two solutions were given.

  • @Sugarman96
    @Sugarman96 Рік тому +15

    I see sinx/x, my mind immediately goes to the Fourier transform. The Fourier transform of a window of 1/2 from -1 to 1 is sinw/w, meaning you can use the pattern of the simple convolution to find the nth convolution of said window with itself, at which point the integral just becomes the inverse Fourier transform of a simple, band limited polynomial function.

    • @emericgent5106
      @emericgent5106 Рік тому +1

      I thought about the exact same thing, do you know a pattern for the n th convolution of the fonction ?

    • @JobBouwman
      @JobBouwman Рік тому +2

      Yes, it's the central frequency of the fourier transform.
      For n = 1 the frequency spectrum is the rect function (one constant function).
      For n = 2 the spectrum is a triangle function (two linear functions).
      For n = 3 the spectrum is a composite of three quadratic functions.
      For n = 4 it's a composite of four cubic functions
      For n it's a piece wise composite of n parts which are (n-1)tic functions.
      I have to think about how to derive the central frequencies.

    • @emericgent5106
      @emericgent5106 Рік тому +1

      @@JobBouwman I also saw that, and calculated the first terms, and understood we just want the value in 0, though we can only do an induction where the whole n-th function is known in the induction hypothesis

  • @illumexhisoka6181
    @illumexhisoka6181 Рік тому +8

    I always wanted someone to do this
    Thank you very much

    • @maths_505
      @maths_505  Рік тому +5

      Me too
      So I decided why not just take it up myself 😂

    • @illumexhisoka6181
      @illumexhisoka6181 Рік тому +1

      @@maths_505 can I ask
      How old are you ?

    • @maths_505
      @maths_505  Рік тому +2

      25

    • @illumexhisoka6181
      @illumexhisoka6181 Рік тому +1

      @@maths_505 I hope that I will be at that level when I am 25
      I mean I understood everything be there is no way I would have thought of doing that
      Great work

    • @erictrefeu5041
      @erictrefeu5041 Рік тому

      ua-cam.com/users/shortsOadiTfmwjTI

  • @ichwillfrieden1635
    @ichwillfrieden1635 Рік тому +3

    This video is so awesome,i learned so much from you thank you

  • @illumexhisoka6181
    @illumexhisoka6181 Рік тому +1

    I turned the integrals into sums where you just need to substitute
    I used a limit so I don't need to change the value of n everywhere in the sum
    The sum doesn't work when n is 1 or 2
    Maybe if you can simplify the sum it would be the key to a more general formula
    This is the sum when n is odd
    DLimit((-1)^((t-1)/2)*t*Pi/2*(1/DProduct((2*y+1)^2-1,y,1,(t-1)/2)+t^(t-2)/DProduct((2*y+1)^2-t^2,y,0,(t-3)/2)+DSum((2*x+1)^(t-2)/(DProduct((2*y+1)^2-(2*x+1)^2,y,0,x-1)*DProduct((2*y+1)^2-(2*x+1)^2,y,x+1,(t-1)/2)),x,1,(t-3)/2)),t,n)
    And this is the sum when n is even
    DLimit((-1)^((t-2)/2)*t*Pi/2*(2^(t-3)/DProduct((2*y)^2-2^2,y,2,t/2)+t^(t-3)/DProduct((2*y)^2-t^2,y,1,(t-2)/2)+DSum((2*x)^(t-3)/(DProduct((2*y)^2-(2*x)^2,y,1,x-1)*DProduct((2*y)^2-(2*x)^2,y,x+1,t/2)),x,2,(t-2)/2)),t,n)
    If you need the sums in a different way of writing or photo of them written tell me
    I wouldn't have done this without your video and without blackpenredpen's way to do partial fraction (the cover up method)

  • @zunaidparker
    @zunaidparker Рік тому +40

    I would be very interested to see a plot of I(n) vs n. How does it behave? Is it monotonic? Does it exhibit interesting patterns? How do the even vs odd cases compare?
    Great video!

    • @erictrefeu5041
      @erictrefeu5041 Рік тому

      ua-cam.com/users/shortsOadiTfmwjTI

    • @svencollister2355
      @svencollister2355 8 місяців тому

      First of all, there are papers on this. Nevertheless, I did exactly this as I found it very interesting as well. I evaluated the integrals up until 70 numerically (because my code runs into some problems I dont understand) and plotted them. Maybe you could do this analytically, but I didnt try at all. If you are interested in this, you can read about that in papers, which reduce this all down to one sum.
      Anyways, when you plot the integral solutions, its not too interesting at all. Its monoton decreasing, maybe it is converging to some value or to 0 for n->inf, but I dont know. Also, it is decreasing slower than 1/x, so the convergence is rather slow. But numerically there are no real patterns for odd and even

    • @svencollister2355
      @svencollister2355 8 місяців тому

      Nvm on the convergence part. For n->inf the integrals converge to 0 as sin(x)/x is bounded by 1, the integral is obviously convergent and therefore, you can interchange the limits. Then the integrand is zero besides at x=0, which does not contribute to the integral

  • @VittorinoPata
    @VittorinoPata Рік тому +1

    Take a look at: A direct computation of a certain
    family of integrals, by
    L. Fornari, E. Laeng and V. Pata
    You will find a more general formula, more explicit, and with a simpler proof.

  • @alankuo5579
    @alankuo5579 Рік тому +1

    You have made my day.

  • @davidblauyoutube
    @davidblauyoutube Рік тому +3

    I have been able, through some small amount of algebra, to reduce the integral to a sum for even n (odd n is no more difficult). This is:
    integral{x=0 to infinity} (sin x/x)^(2n) dx = n pi sum{k=1 to n} (-1)^(n-k) k^(2n-1) / ( (n-k)! (n+k)! ).
    This is always a rational multiple of pi. Multiplying and dividing by (2n)! turns the factorials into a binomial coefficient, which shows that the denominator of the rational multiple (when reduced) always divides (2n)!.
    In particular, when n=1 (so 2n=2), the sum is the single term (1^1/0!2!) = 1/2 and the integral equals pi/2. When n=2 (so 2n=4), the value works out to 2 pi (-1^3/1!3! + 2^3/0!4!) = pi/3. When n=3 (so 2n=6), the value is 3 pi (1^5/2!4! - 2^5/1!5! + 3^5/0!6!) = 11 pi/40. And so on.
    I'm happy to post a PDF with the details.

    • @maths_505
      @maths_505  Рік тому

      I'd love to read that PDF
      My email is in the about section

    • @erictrefeu5041
      @erictrefeu5041 Рік тому

      ua-cam.com/users/shortsOadiTfmwjTI

    • @joaomatos6598
      @joaomatos6598 Рік тому

      And I wold be happy to read it

    • @ShanBojack
      @ShanBojack Рік тому +1

      Share it bro share it

  • @ahsgdf1
    @ahsgdf1 Рік тому +3

    Excellently laid out, as usual, thanks a lot.
    Just a question: I like your path integral solution for n=1, and succeeded applying the same ideas to the case n=2, but I could not solve the case n=3 in this manner. I'd greatly appreciate your comment.

    • @maths_505
      @maths_505  Рік тому +1

      There's a qncubed3 video solving this using contour integration. You should check that out....its quite a nice video

    • @wolfganghintze732
      @wolfganghintze732 Рік тому

      @@maths_505 Thank you for the hint. I am amazed about the high speed with which you use to reply to questions. Great!

  • @Kapomafioso
    @Kapomafioso Рік тому +3

    Homework for n = 69:
    I_69 = 998,343,250,657,696,659,388,623,720,040,379,470,133,597,913,727,156,038,875,228,541,953,757,055,024,051,950,731,143,915,115,744,965,383,517,459,741*pi / 12,027,626,526,020,745,490,674,841,999,023,506,972,927,778,751,265,048,277,098,155,425,840,767,196,171,570,579,020,944,634,806,272,000,000,000,000,000
    I totally didn't use Mathematica to get this. Pinky promise!

  • @jayaadithya
    @jayaadithya Рік тому +4

    Integral 0 to infinity sin^69(x)/x^69 = 0.260765

  • @skyethebi
    @skyethebi 8 місяців тому +1

    A year or so ago when I first watched bprp’s video about the 3rd Dirichlet integral I decided to try and solve a generalized form of the problem and I managed to come up with a formula involving some finite summations. My solution development was pretty bad and not remotely rigorous (I’d only just started learning multivar within the pasts few months and most of what I knew of it still came from UA-cam) and the end results was still pretty ugly but I was pretty happy with myself considering I’d only just learned Feynman’s technique. I should try it again and see if I can get a bit more rigor involved and then maybe try solving generalized fresnel integrals.

    • @skyethebi
      @skyethebi 8 місяців тому

      I might try looking back through my math notes from a year ago to find my solution but I’m not sure I’ll be able to because that’s a lot of notes to check through and I didn’t organize them at all

    • @skyethebi
      @skyethebi 7 місяців тому

      I managed to find it. The answer was:
      For odd n:
      \frac{\pi n}{2^n} \sum_{k=0}^{\frac{n-1}{2}} \frac{(-1)^k(n-2 k)^{n-1}}{k !(n-k) !}
      For even n:
      \frac{\pi n}{2^n} \sum_{k=0}^{\frac{n}{2}} \frac{(-1)^k(n-2 k)^{n-1}}{k !(n-k) !}

    • @skyethebi
      @skyethebi 7 місяців тому

      This wasn’t done at all rigorously so I can leave it to you to actually prove it but it works for all n that I’ve tested. It basically results from doing a power reduction followed by Feynman’s technique going to the n-1 derivative and assuming that the constant after antidifferentiating is always 0 (which is true but I haven’t proved it for the cases where you end up with a cosine in the numerator). Maybe one day I’ll go ahead and do it rigorously but honestly I’m pretty happy to have a solution for the nth Dirichlet integral that only has finite summations.

    • @erictrefeu5041
      @erictrefeu5041 6 місяців тому

      @@skyethebi ua-cam.com/users/shortsDvT_6yRrSOk

  • @erivaldolopes632
    @erivaldolopes632 3 місяці тому

    You could go even further by realising that the integral on the right is a sum of a product of varying coefficients with arctg or pi's

  • @VerSalieri
    @VerSalieri Рік тому

    I really enjoyed this.. thank you.

  • @digxx
    @digxx Рік тому +6

    I think you should be able to calculate the t-integral for general n by the residue theorem.

    • @davidblauyoutube
      @davidblauyoutube Рік тому +4

      I agree, although computing the residues is nasty algebra. I was thinking partial fractions might be easier, since all of the 1/(t^2+n^2) factors integrate nicely to arctans.

    • @attica7980
      @attica7980 Рік тому +5

      @@davidblauyoutube Actually. the residue calculations are fairly easy, since all the zeros in the denominator are simple. In any case, the starting Dirichlet integral can be directly evaluated by the residue theorem. The technique to be used is explained in Ahlfors, Complex Analysis, Section 5.3 (pp. 154-159 in the third edition) Evaluation of Definite Integrals.

  • @nicholaselias9312
    @nicholaselias9312 16 днів тому

    When you isolated the x integral from the t integral, the x integral looked like a Laplace transform.

  • @امینظاهرزاده
    @امینظاهرزاده Рік тому +1

    Please solve
    Integral(-1)^[x]
    Correct component

  • @PyarMatKaro
    @PyarMatKaro Рік тому +2

    I expect that the value (fraction) gets complicated with increasing n but for n=6 it's simply 11 pi / 40

    • @maths_505
      @maths_505  Рік тому

      The partial fractions are mostly arctans so yeah one can expect nice results

    • @erictrefeu5041
      @erictrefeu5041 Рік тому

      ua-cam.com/users/shortsOadiTfmwjTI

  • @shubhamkumar-vx4ld
    @shubhamkumar-vx4ld Рік тому +2

    Hello sir can you do some jee advanced calculus problems these are some of the toughest undergraduate problems

  • @shanmugasundaram9688
    @shanmugasundaram9688 Рік тому

    Excellent calculation.

  • @manstuckinabox3679
    @manstuckinabox3679 Рік тому +7

    YOOOOOOOO WE MADE IT BROOOOO, but we takled the natural final boss, we still need to defeat the real, boss, and of course, dare I say, the c o m p l e x Boss.
    It's actually quite intriguing to think about sin^z(z)/z^z, kind of a weird entity, although I think deviating much from the n case, since we're dealing with variables not constants, and that it's value depends on the parameter we choose, also sin^z(x)/x^z is even more intruiging, I'm no expert in choosing contours but I think we can use a rectangular contour? hmm... actually this looks alot like an exponential so I think it can be takled with a semi-circle integral... sorry it takes time for me to brainstorm lol!
    but awesome video anyways, it felt like an anniversery since I discovered this haven in the begining of chrismass vacation, and now it's about to end...
    btw I starting collecting all the goofy integrals inside a note book of mine called "the Big Book of Integrals" maybe one day we'll collect em all :D
    11:54 I spy a missing minus sign.
    after watching this video I'm really wondering how Matrix transformation and Diagonalization would help in cases like simplifing reduced formulas like this.... idk Just a weird Idea that came about.

  • @ianmi4i727
    @ianmi4i727 Рік тому +1

    This is first class!!!! 🤩

  • @AJ-et3vf
    @AJ-et3vf Рік тому

    Great video. Thank you

  • @SonnyBubba
    @SonnyBubba Рік тому +1

    What hardware and software do you use to make your videos?

    • @maths_505
      @maths_505  Рік тому

      It's the default Samsung notes app.

  • @nolsen1937
    @nolsen1937 Рік тому +1

    This is a NOICE homework?

  • @robertsandy3794
    @robertsandy3794 Рік тому +1

    Awesome

  • @Sky11631
    @Sky11631 Рік тому

    Late comment but could you not reduce the last 2 integrals (for general n) to sums and products using contour integration?

  • @ΑλέξανδροςΖεγγ

    I wonder if contour integration can deal with this integral.

  • @TheoH54
    @TheoH54 Рік тому +1

    Nice! However, after unleashing some of my own tricks on this integral I got a closed form solution, or rather two, one for the n=odd case and one for the n=even case. Both are simple finite sums.
    ADDED: both cases can be combined to obtain a single sum valid for all positive integer n.

    • @erictrefeu5041
      @erictrefeu5041 Рік тому

      ua-cam.com/users/shortsOadiTfmwjTI

    • @TheoH54
      @TheoH54 Рік тому

      @@erictrefeu5041 yes - see my latest comment.

    • @erictrefeu5041
      @erictrefeu5041 Рік тому +1

      @@TheoH54 oui j'ai vu. Nous sommes d'accord Theo. Cette formule est valable pour m pair ou impair. il m'apparaissait utile de compléter cette vidéo avec une formule générale.

    • @TheoH54
      @TheoH54 Рік тому

      ​​​@@erictrefeu5041, bon soir - I would like to ask you, do you have an expression for INT{0;inf} sin(x)^n/x dx ? If yes, maybe we can do something together. I've got a compact expression, apparently not known, with no sums for n=odd.

    • @erictrefeu5041
      @erictrefeu5041 Рік тому +1

      @@TheoH54 Hello Théo, for (sin(x)/x)^m dx ? , yes of course, follow this link : ua-cam.com/users/shortsOadiTfmwjTI
      Or for sin(x)^n/x dx ?... je vais réfléchir, je n'ai pas encore cherché (ca doit etre faisable)

  • @mangeshhebbalkar1715
    @mangeshhebbalkar1715 Рік тому

    I got a form which only involves a finite sum,is there a way I can show it to u

  • @txikitofandango
    @txikitofandango Рік тому +1

    Nice

  • @anestismoutafidis4575
    @anestismoutafidis4575 Рік тому

    => sin1 - (sin-1) •dx=0,0349

  • @trelosyiaellinika
    @trelosyiaellinika 9 місяців тому

    For the case of n being an odd integer, are you sure that the nominator is also n!? It seems to me it should be n!(n-1) instead.

  • @GreenMeansGOF
    @GreenMeansGOF 6 місяців тому +1

    Is there a closed form answer for this integral?

    • @erictrefeu5041
      @erictrefeu5041 6 місяців тому

      ua-cam.com/users/shortsDvT_6yRrSOk

    • @erictrefeu5041
      @erictrefeu5041 6 місяців тому

      ua-cam.com/users/shortsOadiTfmwjTI

  • @aaronrobertcattell8859
    @aaronrobertcattell8859 Рік тому +1

    interesting

  • @nicolascamargo8339
    @nicolascamargo8339 Рік тому

    Wow

  • @Charlesmadeit
    @Charlesmadeit Рік тому +1

    Swag

  • @yoav613
    @yoav613 Рік тому

    The solution for the Hw is: 👅

  • @Anonymous-Indian..2003
    @Anonymous-Indian..2003 Рік тому +7

    Integral 0 to infinity
    n=1 , then π/2
    n=2 , then π/2
    n=3 , then 3π/8
    n=4 , then π/3
    n=5 , then 115π/384
    n=6 , then 11π/40
    n=7 , then 5887π/23040
    I already solved this for any positive n by using only complex analysis and binomial theorem, and I've my own general formula. You've to only give value of n and then you'll get it
    Note: i solved it during my 3rd semester of B.Tech when i was 19 years old.........

  • @rahult1518
    @rahult1518 Рік тому

    when you do IBP , there is 1 mistake in the 2nd column 3rd integral. it should be e^-xt/t^2 ..

  • @erictrefeu5041
    @erictrefeu5041 Рік тому

    j'ai mieux

  • @Outcast_Raj
    @Outcast_Raj Рік тому +1

    Good one bhai

  • @illumexhisoka6181
    @illumexhisoka6181 Рік тому

    I turned the integrals into sums where you just need to substitute
    I used a limit so I don't need to change the value of n everywhere in the sum
    The sum doesn't work when n is 1 or 2
    Maybe if you can simplify the sum it would be the key to a more general formula
    This is the sum when n is odd
    DLimit((-1)^((t-1)/2)*t*Pi/2*(1/DProduct((2*y+1)^2-1,y,1,(t-1)/2)+t^(t-2)/DProduct((2*y+1)^2-t^2,y,0,(t-3)/2)+DSum((2*x+1)^(t-2)/(DProduct((2*y+1)^2-(2*x+1)^2,y,0,x-1)*DProduct((2*y+1)^2-(2*x+1)^2,y,x+1,(t-1)/2)),x,1,(t-3)/2)),t,n)
    And this is the sum when n is even
    DLimit((-1)^((t-2)/2)*t*Pi/2*(2^(t-3)/DProduct((2*y)^2-2^2,y,2,t/2)+t^(t-3)/DProduct((2*y)^2-t^2,y,1,(t-2)/2)+DSum((2*x)^(t-3)/(DProduct((2*y)^2-(2*x)^2,y,1,x-1)*DProduct((2*y)^2-(2*x)^2,y,x+1,t/2)),x,2,(t-2)/2)),t,n)
    If you need the sums in a different way of writing or photo of them written tell me
    I wouldn't have done this without your video and without blackpenredpen's way to do partial fraction (the cover up method)