Explicit solution of a Hamiltonian system












5












$begingroup$


It is well-known that the following Hamiltonian system
begin{eqnarray}
left{begin{array}{rcl}
frac{dx}{dt}&=&y,\
frac{dy}{dt}&=&x(-1+x^2),
end{array}right.
end{eqnarray}

with
$$ H(x,y)=frac{x^2}2+frac{y^2}{2}-frac{x^4}{4}$$
has the solution
$$ x=tanhleft(frac{t}{sqrt2}right),qquad y=frac1{sqrt2}text{sech}^2left(frac{t}{sqrt2}right) $$
such that $H(x,y)=frac14$.



Now consider the following Hamiltonian system
begin{eqnarray}
left{begin{array}{rcl}
frac{dx}{dt}&=&y,\
frac{dy}{dt}&=&x(-1+x^4),
end{array}right.tag{1}
end{eqnarray}

with
$$ H(x,y)=frac{x^2}2+frac{y^2}{2}-frac{x^6}{6}$$
and I want to find the explicit solution such that $H(x,y)=frac{1}{3}$. I tried many ways to get the solution but failed. Any help is appreciated.










share|cite|improve this question









New contributor




xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.







$endgroup$

















    5












    $begingroup$


    It is well-known that the following Hamiltonian system
    begin{eqnarray}
    left{begin{array}{rcl}
    frac{dx}{dt}&=&y,\
    frac{dy}{dt}&=&x(-1+x^2),
    end{array}right.
    end{eqnarray}

    with
    $$ H(x,y)=frac{x^2}2+frac{y^2}{2}-frac{x^4}{4}$$
    has the solution
    $$ x=tanhleft(frac{t}{sqrt2}right),qquad y=frac1{sqrt2}text{sech}^2left(frac{t}{sqrt2}right) $$
    such that $H(x,y)=frac14$.



    Now consider the following Hamiltonian system
    begin{eqnarray}
    left{begin{array}{rcl}
    frac{dx}{dt}&=&y,\
    frac{dy}{dt}&=&x(-1+x^4),
    end{array}right.tag{1}
    end{eqnarray}

    with
    $$ H(x,y)=frac{x^2}2+frac{y^2}{2}-frac{x^6}{6}$$
    and I want to find the explicit solution such that $H(x,y)=frac{1}{3}$. I tried many ways to get the solution but failed. Any help is appreciated.










    share|cite|improve this question









    New contributor




    xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
    Check out our Code of Conduct.







    $endgroup$















      5












      5








      5


      1



      $begingroup$


      It is well-known that the following Hamiltonian system
      begin{eqnarray}
      left{begin{array}{rcl}
      frac{dx}{dt}&=&y,\
      frac{dy}{dt}&=&x(-1+x^2),
      end{array}right.
      end{eqnarray}

      with
      $$ H(x,y)=frac{x^2}2+frac{y^2}{2}-frac{x^4}{4}$$
      has the solution
      $$ x=tanhleft(frac{t}{sqrt2}right),qquad y=frac1{sqrt2}text{sech}^2left(frac{t}{sqrt2}right) $$
      such that $H(x,y)=frac14$.



      Now consider the following Hamiltonian system
      begin{eqnarray}
      left{begin{array}{rcl}
      frac{dx}{dt}&=&y,\
      frac{dy}{dt}&=&x(-1+x^4),
      end{array}right.tag{1}
      end{eqnarray}

      with
      $$ H(x,y)=frac{x^2}2+frac{y^2}{2}-frac{x^6}{6}$$
      and I want to find the explicit solution such that $H(x,y)=frac{1}{3}$. I tried many ways to get the solution but failed. Any help is appreciated.










      share|cite|improve this question









      New contributor




      xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.







      $endgroup$




      It is well-known that the following Hamiltonian system
      begin{eqnarray}
      left{begin{array}{rcl}
      frac{dx}{dt}&=&y,\
      frac{dy}{dt}&=&x(-1+x^2),
      end{array}right.
      end{eqnarray}

      with
      $$ H(x,y)=frac{x^2}2+frac{y^2}{2}-frac{x^4}{4}$$
      has the solution
      $$ x=tanhleft(frac{t}{sqrt2}right),qquad y=frac1{sqrt2}text{sech}^2left(frac{t}{sqrt2}right) $$
      such that $H(x,y)=frac14$.



      Now consider the following Hamiltonian system
      begin{eqnarray}
      left{begin{array}{rcl}
      frac{dx}{dt}&=&y,\
      frac{dy}{dt}&=&x(-1+x^4),
      end{array}right.tag{1}
      end{eqnarray}

      with
      $$ H(x,y)=frac{x^2}2+frac{y^2}{2}-frac{x^6}{6}$$
      and I want to find the explicit solution such that $H(x,y)=frac{1}{3}$. I tried many ways to get the solution but failed. Any help is appreciated.







      dg.differential-geometry ca.classical-analysis-and-odes differential-equations






      share|cite|improve this question









      New contributor




      xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.











      share|cite|improve this question









      New contributor




      xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.









      share|cite|improve this question




      share|cite|improve this question








      edited 2 days ago









      Ivo Terek

      442211




      442211






      New contributor




      xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.









      asked Mar 28 at 17:00









      xpaulxpaul

      1283




      1283




      New contributor




      xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.





      New contributor





      xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.






      xpaul is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
      Check out our Code of Conduct.






















          1 Answer
          1






          active

          oldest

          votes


















          10












          $begingroup$

          Edited.



          Surprisingly, there is an elementary solution, if I made no mistake in
          the following computation.



          Your equation is equivalent to
          $$left(frac{dx}{dt}right)^2=frac{1}{3}(x^6-3x^2+2),$$
          (I just plugged $y=dx/dt$ to your Hamiltonian, and used its value $H=1/3$.)
          This equation is separable,
          $$frac{t}{sqrt{3}}=intfrac{dx}{sqrt{x^6-3x^2+2}}=:I$$
          and requires inversion of the integral.



          To reduce it
          to a standard integral, change $x^2=1/(u+1), ; dx=-(1/2)(u+1)^{-3/2}du,$ and the integral becomes
          $$I=-frac{1}{sqrt{2}}intfrac{du}{usqrt{2u+3}}.$$
          The inverse function to this integral is elementary. Integrating and returning to the original variables, I obtained the general solution
          $$x(t)=sqrt{frac{1+cosh(2t+c)}{2+cosh(2t+c)}}.$$
          Can you check this computation?



          Of course, your value $1/3$ for the Hamiltonian is crucial here. With
          some other value you obtain a much more complicated integral which cannot be expressed in elementary functions. Where did this $1/3$ come from?






          share|cite|improve this answer











          $endgroup$









          • 2




            $begingroup$
            Thanks a lot. This is what I wanted.
            $endgroup$
            – xpaul
            Mar 28 at 21:29










          • $begingroup$
            Note it is much clearer than before. I really appreciate your help.
            $endgroup$
            – xpaul
            2 days ago






          • 1




            $begingroup$
            @xpaul: tell me where your $1/3$ comes from.
            $endgroup$
            – Alexandre Eremenko
            2 days ago






          • 1




            $begingroup$
            1/3 is from fact that when $H(x,y)=alpha$ passes through the separatrices $(pm1,0)$, $alpha=1/3$.
            $endgroup$
            – xpaul
            2 days ago












          • $begingroup$
            If I plug in your $x(t)$, I get $y=x'(t)=frac{sinh (c+2 t)}{sqrt{frac{cosh (c+2 t)+1}{cosh (c+2 t)+2}} (cosh (c+2 t)+2)^2}$ but I find $frac{x^2}2+frac{y^2}{2}-frac{x^6}6neqfrac13$.
            $endgroup$
            – xpaul
            2 days ago












          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: "504"
          };
          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
          });


          }
          });






          xpaul is a new contributor. Be nice, and check out our Code of Conduct.










          draft saved

          draft discarded


















          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathoverflow.net%2fquestions%2f326601%2fexplicit-solution-of-a-hamiltonian-system%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









          10












          $begingroup$

          Edited.



          Surprisingly, there is an elementary solution, if I made no mistake in
          the following computation.



          Your equation is equivalent to
          $$left(frac{dx}{dt}right)^2=frac{1}{3}(x^6-3x^2+2),$$
          (I just plugged $y=dx/dt$ to your Hamiltonian, and used its value $H=1/3$.)
          This equation is separable,
          $$frac{t}{sqrt{3}}=intfrac{dx}{sqrt{x^6-3x^2+2}}=:I$$
          and requires inversion of the integral.



          To reduce it
          to a standard integral, change $x^2=1/(u+1), ; dx=-(1/2)(u+1)^{-3/2}du,$ and the integral becomes
          $$I=-frac{1}{sqrt{2}}intfrac{du}{usqrt{2u+3}}.$$
          The inverse function to this integral is elementary. Integrating and returning to the original variables, I obtained the general solution
          $$x(t)=sqrt{frac{1+cosh(2t+c)}{2+cosh(2t+c)}}.$$
          Can you check this computation?



          Of course, your value $1/3$ for the Hamiltonian is crucial here. With
          some other value you obtain a much more complicated integral which cannot be expressed in elementary functions. Where did this $1/3$ come from?






          share|cite|improve this answer











          $endgroup$









          • 2




            $begingroup$
            Thanks a lot. This is what I wanted.
            $endgroup$
            – xpaul
            Mar 28 at 21:29










          • $begingroup$
            Note it is much clearer than before. I really appreciate your help.
            $endgroup$
            – xpaul
            2 days ago






          • 1




            $begingroup$
            @xpaul: tell me where your $1/3$ comes from.
            $endgroup$
            – Alexandre Eremenko
            2 days ago






          • 1




            $begingroup$
            1/3 is from fact that when $H(x,y)=alpha$ passes through the separatrices $(pm1,0)$, $alpha=1/3$.
            $endgroup$
            – xpaul
            2 days ago












          • $begingroup$
            If I plug in your $x(t)$, I get $y=x'(t)=frac{sinh (c+2 t)}{sqrt{frac{cosh (c+2 t)+1}{cosh (c+2 t)+2}} (cosh (c+2 t)+2)^2}$ but I find $frac{x^2}2+frac{y^2}{2}-frac{x^6}6neqfrac13$.
            $endgroup$
            – xpaul
            2 days ago
















          10












          $begingroup$

          Edited.



          Surprisingly, there is an elementary solution, if I made no mistake in
          the following computation.



          Your equation is equivalent to
          $$left(frac{dx}{dt}right)^2=frac{1}{3}(x^6-3x^2+2),$$
          (I just plugged $y=dx/dt$ to your Hamiltonian, and used its value $H=1/3$.)
          This equation is separable,
          $$frac{t}{sqrt{3}}=intfrac{dx}{sqrt{x^6-3x^2+2}}=:I$$
          and requires inversion of the integral.



          To reduce it
          to a standard integral, change $x^2=1/(u+1), ; dx=-(1/2)(u+1)^{-3/2}du,$ and the integral becomes
          $$I=-frac{1}{sqrt{2}}intfrac{du}{usqrt{2u+3}}.$$
          The inverse function to this integral is elementary. Integrating and returning to the original variables, I obtained the general solution
          $$x(t)=sqrt{frac{1+cosh(2t+c)}{2+cosh(2t+c)}}.$$
          Can you check this computation?



          Of course, your value $1/3$ for the Hamiltonian is crucial here. With
          some other value you obtain a much more complicated integral which cannot be expressed in elementary functions. Where did this $1/3$ come from?






          share|cite|improve this answer











          $endgroup$









          • 2




            $begingroup$
            Thanks a lot. This is what I wanted.
            $endgroup$
            – xpaul
            Mar 28 at 21:29










          • $begingroup$
            Note it is much clearer than before. I really appreciate your help.
            $endgroup$
            – xpaul
            2 days ago






          • 1




            $begingroup$
            @xpaul: tell me where your $1/3$ comes from.
            $endgroup$
            – Alexandre Eremenko
            2 days ago






          • 1




            $begingroup$
            1/3 is from fact that when $H(x,y)=alpha$ passes through the separatrices $(pm1,0)$, $alpha=1/3$.
            $endgroup$
            – xpaul
            2 days ago












          • $begingroup$
            If I plug in your $x(t)$, I get $y=x'(t)=frac{sinh (c+2 t)}{sqrt{frac{cosh (c+2 t)+1}{cosh (c+2 t)+2}} (cosh (c+2 t)+2)^2}$ but I find $frac{x^2}2+frac{y^2}{2}-frac{x^6}6neqfrac13$.
            $endgroup$
            – xpaul
            2 days ago














          10












          10








          10





          $begingroup$

          Edited.



          Surprisingly, there is an elementary solution, if I made no mistake in
          the following computation.



          Your equation is equivalent to
          $$left(frac{dx}{dt}right)^2=frac{1}{3}(x^6-3x^2+2),$$
          (I just plugged $y=dx/dt$ to your Hamiltonian, and used its value $H=1/3$.)
          This equation is separable,
          $$frac{t}{sqrt{3}}=intfrac{dx}{sqrt{x^6-3x^2+2}}=:I$$
          and requires inversion of the integral.



          To reduce it
          to a standard integral, change $x^2=1/(u+1), ; dx=-(1/2)(u+1)^{-3/2}du,$ and the integral becomes
          $$I=-frac{1}{sqrt{2}}intfrac{du}{usqrt{2u+3}}.$$
          The inverse function to this integral is elementary. Integrating and returning to the original variables, I obtained the general solution
          $$x(t)=sqrt{frac{1+cosh(2t+c)}{2+cosh(2t+c)}}.$$
          Can you check this computation?



          Of course, your value $1/3$ for the Hamiltonian is crucial here. With
          some other value you obtain a much more complicated integral which cannot be expressed in elementary functions. Where did this $1/3$ come from?






          share|cite|improve this answer











          $endgroup$



          Edited.



          Surprisingly, there is an elementary solution, if I made no mistake in
          the following computation.



          Your equation is equivalent to
          $$left(frac{dx}{dt}right)^2=frac{1}{3}(x^6-3x^2+2),$$
          (I just plugged $y=dx/dt$ to your Hamiltonian, and used its value $H=1/3$.)
          This equation is separable,
          $$frac{t}{sqrt{3}}=intfrac{dx}{sqrt{x^6-3x^2+2}}=:I$$
          and requires inversion of the integral.



          To reduce it
          to a standard integral, change $x^2=1/(u+1), ; dx=-(1/2)(u+1)^{-3/2}du,$ and the integral becomes
          $$I=-frac{1}{sqrt{2}}intfrac{du}{usqrt{2u+3}}.$$
          The inverse function to this integral is elementary. Integrating and returning to the original variables, I obtained the general solution
          $$x(t)=sqrt{frac{1+cosh(2t+c)}{2+cosh(2t+c)}}.$$
          Can you check this computation?



          Of course, your value $1/3$ for the Hamiltonian is crucial here. With
          some other value you obtain a much more complicated integral which cannot be expressed in elementary functions. Where did this $1/3$ come from?







          share|cite|improve this answer














          share|cite|improve this answer



          share|cite|improve this answer








          edited Mar 29 at 1:33

























          answered Mar 28 at 18:42









          Alexandre EremenkoAlexandre Eremenko

          51k6141261




          51k6141261








          • 2




            $begingroup$
            Thanks a lot. This is what I wanted.
            $endgroup$
            – xpaul
            Mar 28 at 21:29










          • $begingroup$
            Note it is much clearer than before. I really appreciate your help.
            $endgroup$
            – xpaul
            2 days ago






          • 1




            $begingroup$
            @xpaul: tell me where your $1/3$ comes from.
            $endgroup$
            – Alexandre Eremenko
            2 days ago






          • 1




            $begingroup$
            1/3 is from fact that when $H(x,y)=alpha$ passes through the separatrices $(pm1,0)$, $alpha=1/3$.
            $endgroup$
            – xpaul
            2 days ago












          • $begingroup$
            If I plug in your $x(t)$, I get $y=x'(t)=frac{sinh (c+2 t)}{sqrt{frac{cosh (c+2 t)+1}{cosh (c+2 t)+2}} (cosh (c+2 t)+2)^2}$ but I find $frac{x^2}2+frac{y^2}{2}-frac{x^6}6neqfrac13$.
            $endgroup$
            – xpaul
            2 days ago














          • 2




            $begingroup$
            Thanks a lot. This is what I wanted.
            $endgroup$
            – xpaul
            Mar 28 at 21:29










          • $begingroup$
            Note it is much clearer than before. I really appreciate your help.
            $endgroup$
            – xpaul
            2 days ago






          • 1




            $begingroup$
            @xpaul: tell me where your $1/3$ comes from.
            $endgroup$
            – Alexandre Eremenko
            2 days ago






          • 1




            $begingroup$
            1/3 is from fact that when $H(x,y)=alpha$ passes through the separatrices $(pm1,0)$, $alpha=1/3$.
            $endgroup$
            – xpaul
            2 days ago












          • $begingroup$
            If I plug in your $x(t)$, I get $y=x'(t)=frac{sinh (c+2 t)}{sqrt{frac{cosh (c+2 t)+1}{cosh (c+2 t)+2}} (cosh (c+2 t)+2)^2}$ but I find $frac{x^2}2+frac{y^2}{2}-frac{x^6}6neqfrac13$.
            $endgroup$
            – xpaul
            2 days ago








          2




          2




          $begingroup$
          Thanks a lot. This is what I wanted.
          $endgroup$
          – xpaul
          Mar 28 at 21:29




          $begingroup$
          Thanks a lot. This is what I wanted.
          $endgroup$
          – xpaul
          Mar 28 at 21:29












          $begingroup$
          Note it is much clearer than before. I really appreciate your help.
          $endgroup$
          – xpaul
          2 days ago




          $begingroup$
          Note it is much clearer than before. I really appreciate your help.
          $endgroup$
          – xpaul
          2 days ago




          1




          1




          $begingroup$
          @xpaul: tell me where your $1/3$ comes from.
          $endgroup$
          – Alexandre Eremenko
          2 days ago




          $begingroup$
          @xpaul: tell me where your $1/3$ comes from.
          $endgroup$
          – Alexandre Eremenko
          2 days ago




          1




          1




          $begingroup$
          1/3 is from fact that when $H(x,y)=alpha$ passes through the separatrices $(pm1,0)$, $alpha=1/3$.
          $endgroup$
          – xpaul
          2 days ago






          $begingroup$
          1/3 is from fact that when $H(x,y)=alpha$ passes through the separatrices $(pm1,0)$, $alpha=1/3$.
          $endgroup$
          – xpaul
          2 days ago














          $begingroup$
          If I plug in your $x(t)$, I get $y=x'(t)=frac{sinh (c+2 t)}{sqrt{frac{cosh (c+2 t)+1}{cosh (c+2 t)+2}} (cosh (c+2 t)+2)^2}$ but I find $frac{x^2}2+frac{y^2}{2}-frac{x^6}6neqfrac13$.
          $endgroup$
          – xpaul
          2 days ago




          $begingroup$
          If I plug in your $x(t)$, I get $y=x'(t)=frac{sinh (c+2 t)}{sqrt{frac{cosh (c+2 t)+1}{cosh (c+2 t)+2}} (cosh (c+2 t)+2)^2}$ but I find $frac{x^2}2+frac{y^2}{2}-frac{x^6}6neqfrac13$.
          $endgroup$
          – xpaul
          2 days ago










          xpaul is a new contributor. Be nice, and check out our Code of Conduct.










          draft saved

          draft discarded


















          xpaul is a new contributor. Be nice, and check out our Code of Conduct.













          xpaul is a new contributor. Be nice, and check out our Code of Conduct.












          xpaul is a new contributor. Be nice, and check out our Code of Conduct.
















          Thanks for contributing an answer to MathOverflow!


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




          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathoverflow.net%2fquestions%2f326601%2fexplicit-solution-of-a-hamiltonian-system%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

          數位音樂下載

          格利澤436b

          When can things happen in Etherscan, such as the picture below?