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

Поделиться
HTML-код
  • Опубликовано: 15 дек 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 Год назад +13

    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)

  • @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 Год назад

      ruclips.net/user/shortsOadiTfmwjTI

  • @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

  • @vladimirlucic1276
    @vladimirlucic1276 Год назад +12

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

  • @ichwillfrieden1635
    @ichwillfrieden1635 Год назад +3

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

  • @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 Год назад

      ruclips.net/user/shortsOadiTfmwjTI

    • @svencollister2355
      @svencollister2355 9 месяцев назад

      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 9 месяцев назад

      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

  • @alankuo5579
    @alankuo5579 Год назад +1

    You have made my day.

  • @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!

  • @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 RUclips) 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 8 месяцев назад

      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 8 месяцев назад

      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 ruclips.net/user/shortsDvT_6yRrSOk

  • @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.

  • @jayaadithya
    @jayaadithya Год назад +4

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

  • @VerSalieri
    @VerSalieri Год назад

    I really enjoyed this.. thank you.

  • @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 Год назад

      ruclips.net/user/shortsOadiTfmwjTI

    • @joaomatos6598
      @joaomatos6598 Год назад

      And I wold be happy to read it

    • @ShanBojack
      @ShanBojack Год назад +1

      Share it bro share it

  • @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)

  • @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!

  • @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.

  • @erivaldolopes632
    @erivaldolopes632 4 месяца назад

    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

  • @shanmugasundaram9688
    @shanmugasundaram9688 Год назад

    Excellent calculation.

  • @امینظاهرزاده
    @امینظاهرزاده Год назад +1

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

  • @ianmi4i727
    @ianmi4i727 Год назад +1

    This is first class!!!! 🤩

  • @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.

  • @shubhamkumar-vx4ld
    @shubhamkumar-vx4ld Год назад +2

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

  • @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.

  • @nicholaselias9312
    @nicholaselias9312 Месяц назад

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

  • @AJ-et3vf
    @AJ-et3vf Год назад

    Great video. Thank you

  • @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 Год назад

      ruclips.net/user/shortsOadiTfmwjTI

  • @robertsandy3794
    @robertsandy3794 Год назад +1

    Awesome

  • @GreenMeansGOF
    @GreenMeansGOF 7 месяцев назад +1

    Is there a closed form answer for this integral?

    • @erictrefeu5041
      @erictrefeu5041 6 месяцев назад

      ruclips.net/user/shortsDvT_6yRrSOk

    • @erictrefeu5041
      @erictrefeu5041 6 месяцев назад

      ruclips.net/user/shortsOadiTfmwjTI

  • @Sky11631
    @Sky11631 Год назад

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

  • @mangeshhebbalkar1715
    @mangeshhebbalkar1715 Год назад

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

  • @ΑλέξανδροςΖεγγ
    @ΑλέξανδροςΖεγγ Год назад

    I wonder if contour integration can deal with this integral.

  • @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.

  • @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 Год назад

      ruclips.net/user/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 : ruclips.net/user/shortsOadiTfmwjTI
      Or for sin(x)^n/x dx ?... je vais réfléchir, je n'ai pas encore cherché (ca doit etre faisable)

  • @txikitofandango
    @txikitofandango Год назад +1

    Nice

  • @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.........

  • @nolsen1937
    @nolsen1937 Год назад +1

    This is a NOICE homework?

  • @aaronrobertcattell8859
    @aaronrobertcattell8859 Год назад +1

    interesting

  • @nicolascamargo8339
    @nicolascamargo8339 Год назад

    Wow

  • @anestismoutafidis4575
    @anestismoutafidis4575 Год назад

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

  • @Charlesmadeit
    @Charlesmadeit Год назад +1

    Swag

  • @yoav613
    @yoav613 Год назад

    The solution for the Hw is: 👅

  • @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)