Can we refine this asymptotic for Laguerre polynomials?












4














I just found an interesting and useful limit for Laguerre polynomials:



$$lim_{n to infty} L_n left( frac{2r}{n+1/2} right)=J_0(2 sqrt{2r})$$



I'm using specifically this form of the argument because it's the one I'm working with in the application. Of course, we can set any fixed number instead of $1/2$ in the denominator.



I found this limit in a paper, which references G. Szego. Orthogonal Polynomials. Amer. Math. Soc. Colloq. Publ. 23, Amer. Math. Soc. Providence,
RI, 1975. Fourth Edition.
, , Theorem 8.1.3.



While the limit is useful for very large orders and smallish $r$, I would really like to know if there's a refinement that could be applied to derive an asymptotic expansion, which would depend on $n$.



Here's an illustration which shows that the limit is not that good for larger $r$ (though it does approximate the roots better than the magnitude):



enter image description here



Not sure how we could refine this asymptotic or how the original limit was derived (as I don't have the linked book).



One way is considering the differential equations for both functions.



There's also an interesting result from Gradshteyn-Ryzhik:



$$L_n(z)= frac{2}{n!} e^z int_0^infty e^{-t^2} t^{2n+1} J_0(2t sqrt{z}) dt$$



Which may or may not be related to the limit above.










share|cite|improve this question



























    4














    I just found an interesting and useful limit for Laguerre polynomials:



    $$lim_{n to infty} L_n left( frac{2r}{n+1/2} right)=J_0(2 sqrt{2r})$$



    I'm using specifically this form of the argument because it's the one I'm working with in the application. Of course, we can set any fixed number instead of $1/2$ in the denominator.



    I found this limit in a paper, which references G. Szego. Orthogonal Polynomials. Amer. Math. Soc. Colloq. Publ. 23, Amer. Math. Soc. Providence,
    RI, 1975. Fourth Edition.
    , , Theorem 8.1.3.



    While the limit is useful for very large orders and smallish $r$, I would really like to know if there's a refinement that could be applied to derive an asymptotic expansion, which would depend on $n$.



    Here's an illustration which shows that the limit is not that good for larger $r$ (though it does approximate the roots better than the magnitude):



    enter image description here



    Not sure how we could refine this asymptotic or how the original limit was derived (as I don't have the linked book).



    One way is considering the differential equations for both functions.



    There's also an interesting result from Gradshteyn-Ryzhik:



    $$L_n(z)= frac{2}{n!} e^z int_0^infty e^{-t^2} t^{2n+1} J_0(2t sqrt{z}) dt$$



    Which may or may not be related to the limit above.










    share|cite|improve this question

























      4












      4








      4


      2





      I just found an interesting and useful limit for Laguerre polynomials:



      $$lim_{n to infty} L_n left( frac{2r}{n+1/2} right)=J_0(2 sqrt{2r})$$



      I'm using specifically this form of the argument because it's the one I'm working with in the application. Of course, we can set any fixed number instead of $1/2$ in the denominator.



      I found this limit in a paper, which references G. Szego. Orthogonal Polynomials. Amer. Math. Soc. Colloq. Publ. 23, Amer. Math. Soc. Providence,
      RI, 1975. Fourth Edition.
      , , Theorem 8.1.3.



      While the limit is useful for very large orders and smallish $r$, I would really like to know if there's a refinement that could be applied to derive an asymptotic expansion, which would depend on $n$.



      Here's an illustration which shows that the limit is not that good for larger $r$ (though it does approximate the roots better than the magnitude):



      enter image description here



      Not sure how we could refine this asymptotic or how the original limit was derived (as I don't have the linked book).



      One way is considering the differential equations for both functions.



      There's also an interesting result from Gradshteyn-Ryzhik:



      $$L_n(z)= frac{2}{n!} e^z int_0^infty e^{-t^2} t^{2n+1} J_0(2t sqrt{z}) dt$$



      Which may or may not be related to the limit above.










      share|cite|improve this question













      I just found an interesting and useful limit for Laguerre polynomials:



      $$lim_{n to infty} L_n left( frac{2r}{n+1/2} right)=J_0(2 sqrt{2r})$$



      I'm using specifically this form of the argument because it's the one I'm working with in the application. Of course, we can set any fixed number instead of $1/2$ in the denominator.



      I found this limit in a paper, which references G. Szego. Orthogonal Polynomials. Amer. Math. Soc. Colloq. Publ. 23, Amer. Math. Soc. Providence,
      RI, 1975. Fourth Edition.
      , , Theorem 8.1.3.



      While the limit is useful for very large orders and smallish $r$, I would really like to know if there's a refinement that could be applied to derive an asymptotic expansion, which would depend on $n$.



      Here's an illustration which shows that the limit is not that good for larger $r$ (though it does approximate the roots better than the magnitude):



      enter image description here



      Not sure how we could refine this asymptotic or how the original limit was derived (as I don't have the linked book).



      One way is considering the differential equations for both functions.



      There's also an interesting result from Gradshteyn-Ryzhik:



      $$L_n(z)= frac{2}{n!} e^z int_0^infty e^{-t^2} t^{2n+1} J_0(2t sqrt{z}) dt$$



      Which may or may not be related to the limit above.







      limits asymptotics approximation bessel-functions orthogonal-polynomials






      share|cite|improve this question













      share|cite|improve this question











      share|cite|improve this question




      share|cite|improve this question










      asked Nov 26 at 23:05









      Yuriy S

      15.7k433117




      15.7k433117






















          1 Answer
          1






          active

          oldest

          votes


















          2














          From the given representation, we can express
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{2}{n!} e^{frac{2r}{n+1/2}} int_0^infty e^{-t^2} t^{2n+1} J_0left(2t sqrt{frac{2r}{n+1/2}}right) ,dt
          end{equation}

          changing $ t=sqrt{uleft( n+1/2 right)}$,
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-uleft( n+1/2 right)}u^{n+1/2}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          or
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-left( n+1/2 right)left( u-ln u right)}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          The argument of the exponential is minimum at $u=1$. We can use the Laplace method to derive an asymptotic approximation for the integral. Near $u=1$,
          $u-ln(u)sim 1+(u-1)^2/2$ and $u^{-1/2}J_0left( 2sqrt{2ru}right)sim J_0left( 2sqrt{2r}right)$, then
          begin{equation}
          L_n(frac{2r}{n+1/2})sim sqrt{2pi}frac{(n+1/2)^{n+1/2}}{n!} e^{frac{2r}{n+1/2}}
          e^{-(n+1/2)}J_0left( 2sqrt{2r}right)
          end{equation}

          Now, plugging the Stirling approximation for $n!$ and an expansion for the exponential term, we obtain the expected result $ L_n(frac{2r}{n+1/2})sim J_0left( 2sqrt{2r}right)$. To improve the approximation, the above formula and/or higher orders in the Laplace expansion can be used.






          share|cite|improve this answer























          • Thank you! I will wait a little while before accepting the answer, but this is exactly what I wanted
            – Yuriy S
            Nov 30 at 1:02











          Your Answer





          StackExchange.ifUsing("editor", function () {
          return StackExchange.using("mathjaxEditing", function () {
          StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
          StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
          });
          });
          }, "mathjax-editing");

          StackExchange.ready(function() {
          var channelOptions = {
          tags: "".split(" "),
          id: "69"
          };
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function() {
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled) {
          StackExchange.using("snippets", function() {
          createEditor();
          });
          }
          else {
          createEditor();
          }
          });

          function createEditor() {
          StackExchange.prepareEditor({
          heartbeatType: 'answer',
          autoActivateHeartbeat: false,
          convertImagesToLinks: true,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          imageUploader: {
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          },
          noCode: true, onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          });


          }
          });














          draft saved

          draft discarded


















          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3015077%2fcan-we-refine-this-asymptotic-for-laguerre-polynomials%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown

























          1 Answer
          1






          active

          oldest

          votes








          1 Answer
          1






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes









          2














          From the given representation, we can express
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{2}{n!} e^{frac{2r}{n+1/2}} int_0^infty e^{-t^2} t^{2n+1} J_0left(2t sqrt{frac{2r}{n+1/2}}right) ,dt
          end{equation}

          changing $ t=sqrt{uleft( n+1/2 right)}$,
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-uleft( n+1/2 right)}u^{n+1/2}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          or
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-left( n+1/2 right)left( u-ln u right)}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          The argument of the exponential is minimum at $u=1$. We can use the Laplace method to derive an asymptotic approximation for the integral. Near $u=1$,
          $u-ln(u)sim 1+(u-1)^2/2$ and $u^{-1/2}J_0left( 2sqrt{2ru}right)sim J_0left( 2sqrt{2r}right)$, then
          begin{equation}
          L_n(frac{2r}{n+1/2})sim sqrt{2pi}frac{(n+1/2)^{n+1/2}}{n!} e^{frac{2r}{n+1/2}}
          e^{-(n+1/2)}J_0left( 2sqrt{2r}right)
          end{equation}

          Now, plugging the Stirling approximation for $n!$ and an expansion for the exponential term, we obtain the expected result $ L_n(frac{2r}{n+1/2})sim J_0left( 2sqrt{2r}right)$. To improve the approximation, the above formula and/or higher orders in the Laplace expansion can be used.






          share|cite|improve this answer























          • Thank you! I will wait a little while before accepting the answer, but this is exactly what I wanted
            – Yuriy S
            Nov 30 at 1:02
















          2














          From the given representation, we can express
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{2}{n!} e^{frac{2r}{n+1/2}} int_0^infty e^{-t^2} t^{2n+1} J_0left(2t sqrt{frac{2r}{n+1/2}}right) ,dt
          end{equation}

          changing $ t=sqrt{uleft( n+1/2 right)}$,
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-uleft( n+1/2 right)}u^{n+1/2}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          or
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-left( n+1/2 right)left( u-ln u right)}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          The argument of the exponential is minimum at $u=1$. We can use the Laplace method to derive an asymptotic approximation for the integral. Near $u=1$,
          $u-ln(u)sim 1+(u-1)^2/2$ and $u^{-1/2}J_0left( 2sqrt{2ru}right)sim J_0left( 2sqrt{2r}right)$, then
          begin{equation}
          L_n(frac{2r}{n+1/2})sim sqrt{2pi}frac{(n+1/2)^{n+1/2}}{n!} e^{frac{2r}{n+1/2}}
          e^{-(n+1/2)}J_0left( 2sqrt{2r}right)
          end{equation}

          Now, plugging the Stirling approximation for $n!$ and an expansion for the exponential term, we obtain the expected result $ L_n(frac{2r}{n+1/2})sim J_0left( 2sqrt{2r}right)$. To improve the approximation, the above formula and/or higher orders in the Laplace expansion can be used.






          share|cite|improve this answer























          • Thank you! I will wait a little while before accepting the answer, but this is exactly what I wanted
            – Yuriy S
            Nov 30 at 1:02














          2












          2








          2






          From the given representation, we can express
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{2}{n!} e^{frac{2r}{n+1/2}} int_0^infty e^{-t^2} t^{2n+1} J_0left(2t sqrt{frac{2r}{n+1/2}}right) ,dt
          end{equation}

          changing $ t=sqrt{uleft( n+1/2 right)}$,
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-uleft( n+1/2 right)}u^{n+1/2}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          or
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-left( n+1/2 right)left( u-ln u right)}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          The argument of the exponential is minimum at $u=1$. We can use the Laplace method to derive an asymptotic approximation for the integral. Near $u=1$,
          $u-ln(u)sim 1+(u-1)^2/2$ and $u^{-1/2}J_0left( 2sqrt{2ru}right)sim J_0left( 2sqrt{2r}right)$, then
          begin{equation}
          L_n(frac{2r}{n+1/2})sim sqrt{2pi}frac{(n+1/2)^{n+1/2}}{n!} e^{frac{2r}{n+1/2}}
          e^{-(n+1/2)}J_0left( 2sqrt{2r}right)
          end{equation}

          Now, plugging the Stirling approximation for $n!$ and an expansion for the exponential term, we obtain the expected result $ L_n(frac{2r}{n+1/2})sim J_0left( 2sqrt{2r}right)$. To improve the approximation, the above formula and/or higher orders in the Laplace expansion can be used.






          share|cite|improve this answer














          From the given representation, we can express
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{2}{n!} e^{frac{2r}{n+1/2}} int_0^infty e^{-t^2} t^{2n+1} J_0left(2t sqrt{frac{2r}{n+1/2}}right) ,dt
          end{equation}

          changing $ t=sqrt{uleft( n+1/2 right)}$,
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-uleft( n+1/2 right)}u^{n+1/2}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          or
          begin{equation}
          L_n(frac{2r}{n+1/2})= frac{(n+1/2)^{n+1}}{n!} e^{frac{2r}{n+1/2}}int_0^infty e^{-left( n+1/2 right)left( u-ln u right)}J_0left( 2sqrt{2ru} right)frac{du}{sqrt{u}}
          end{equation}

          The argument of the exponential is minimum at $u=1$. We can use the Laplace method to derive an asymptotic approximation for the integral. Near $u=1$,
          $u-ln(u)sim 1+(u-1)^2/2$ and $u^{-1/2}J_0left( 2sqrt{2ru}right)sim J_0left( 2sqrt{2r}right)$, then
          begin{equation}
          L_n(frac{2r}{n+1/2})sim sqrt{2pi}frac{(n+1/2)^{n+1/2}}{n!} e^{frac{2r}{n+1/2}}
          e^{-(n+1/2)}J_0left( 2sqrt{2r}right)
          end{equation}

          Now, plugging the Stirling approximation for $n!$ and an expansion for the exponential term, we obtain the expected result $ L_n(frac{2r}{n+1/2})sim J_0left( 2sqrt{2r}right)$. To improve the approximation, the above formula and/or higher orders in the Laplace expansion can be used.







          share|cite|improve this answer














          share|cite|improve this answer



          share|cite|improve this answer








          edited Nov 29 at 23:07

























          answered Nov 29 at 22:47









          Paul Enta

          4,15311129




          4,15311129












          • Thank you! I will wait a little while before accepting the answer, but this is exactly what I wanted
            – Yuriy S
            Nov 30 at 1:02


















          • Thank you! I will wait a little while before accepting the answer, but this is exactly what I wanted
            – Yuriy S
            Nov 30 at 1:02
















          Thank you! I will wait a little while before accepting the answer, but this is exactly what I wanted
          – Yuriy S
          Nov 30 at 1:02




          Thank you! I will wait a little while before accepting the answer, but this is exactly what I wanted
          – Yuriy S
          Nov 30 at 1:02


















          draft saved

          draft discarded




















































          Thanks for contributing an answer to Mathematics Stack Exchange!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          Use MathJax to format equations. MathJax reference.


          To learn more, see our tips on writing great answers.





          Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


          Please pay close attention to the following guidance:


          • Please be sure to answer the question. Provide details and share your research!

          But avoid



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3015077%2fcan-we-refine-this-asymptotic-for-laguerre-polynomials%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          Bundesstraße 106

          Le Mesnil-Réaume

          Ida-Boy-Ed-Garten