Runge-Kutta method for higher-order differential equations











up vote
0
down vote

favorite












I am studying Numerical Analysis with the book of Richard L.Burden.
A question which I'm struggling with right now is following.



Transform the second-order initial-value problem



$y'' - 2y' + 2y = e^{2t}sint$ for $0 leq t leq 1, $ with $y(0) = -0.4, y'(0) = -0.6, h=0.1$



into a system of ifrst order initial-value problems, and use the Runge-Kutta method ith h=0.1 to approximate the solution.



Then,
$$u_1(t) = y(t), u_2(t) = y'(t)$$
$$u_1'(t) = u_2(t)$$
$$u_2'(t) = e^{2t}sint - 2u_1(t) + u_2(t)$$
$$u_1(0) = -0.4, u_2(0) = -0.6$$



This initial conditions give $w_{1,0} = -0.4, w_{2,0}=-0.6$



I can understand that $k_{1,1} = hf_1(t_0, w_{1,0}, w_{2,0}) = hw_{2,0}$



$f_1 = u_1'= u_2(t)$,
So $f_1(t_0, w_{1,0}, w_{2,0}) = u_2(t_0, w_{1,0}, w_{2,0}) = w_{2,0}$ (By definition of $w_{i,j}$)



However, I can't understand the following.
$$k_{2,1} = hf_1(t_0 + frac{h}{2}, w_{1,0} + frac{1}{2}k_{1,1}, w_{2,0} + frac{1}{2}k_{1,2}) = hleft[w_{2,0} + frac{1}{2}k_{1,2}right]$$



Why does $f_1(t_0 + frac{h}{2}, w_{1,0} + frac{1}{2}k_{1,1}, w_{2,0} + frac{1}{2}k_{1,2})$ equal to $w_{2,0} + frac{1}{2}k_{1,2}$? It seems that third argument in the function comes out, but there is no detailed explanation in this book.










share|cite|improve this question







New contributor




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
























    up vote
    0
    down vote

    favorite












    I am studying Numerical Analysis with the book of Richard L.Burden.
    A question which I'm struggling with right now is following.



    Transform the second-order initial-value problem



    $y'' - 2y' + 2y = e^{2t}sint$ for $0 leq t leq 1, $ with $y(0) = -0.4, y'(0) = -0.6, h=0.1$



    into a system of ifrst order initial-value problems, and use the Runge-Kutta method ith h=0.1 to approximate the solution.



    Then,
    $$u_1(t) = y(t), u_2(t) = y'(t)$$
    $$u_1'(t) = u_2(t)$$
    $$u_2'(t) = e^{2t}sint - 2u_1(t) + u_2(t)$$
    $$u_1(0) = -0.4, u_2(0) = -0.6$$



    This initial conditions give $w_{1,0} = -0.4, w_{2,0}=-0.6$



    I can understand that $k_{1,1} = hf_1(t_0, w_{1,0}, w_{2,0}) = hw_{2,0}$



    $f_1 = u_1'= u_2(t)$,
    So $f_1(t_0, w_{1,0}, w_{2,0}) = u_2(t_0, w_{1,0}, w_{2,0}) = w_{2,0}$ (By definition of $w_{i,j}$)



    However, I can't understand the following.
    $$k_{2,1} = hf_1(t_0 + frac{h}{2}, w_{1,0} + frac{1}{2}k_{1,1}, w_{2,0} + frac{1}{2}k_{1,2}) = hleft[w_{2,0} + frac{1}{2}k_{1,2}right]$$



    Why does $f_1(t_0 + frac{h}{2}, w_{1,0} + frac{1}{2}k_{1,1}, w_{2,0} + frac{1}{2}k_{1,2})$ equal to $w_{2,0} + frac{1}{2}k_{1,2}$? It seems that third argument in the function comes out, but there is no detailed explanation in this book.










    share|cite|improve this question







    New contributor




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






















      up vote
      0
      down vote

      favorite









      up vote
      0
      down vote

      favorite











      I am studying Numerical Analysis with the book of Richard L.Burden.
      A question which I'm struggling with right now is following.



      Transform the second-order initial-value problem



      $y'' - 2y' + 2y = e^{2t}sint$ for $0 leq t leq 1, $ with $y(0) = -0.4, y'(0) = -0.6, h=0.1$



      into a system of ifrst order initial-value problems, and use the Runge-Kutta method ith h=0.1 to approximate the solution.



      Then,
      $$u_1(t) = y(t), u_2(t) = y'(t)$$
      $$u_1'(t) = u_2(t)$$
      $$u_2'(t) = e^{2t}sint - 2u_1(t) + u_2(t)$$
      $$u_1(0) = -0.4, u_2(0) = -0.6$$



      This initial conditions give $w_{1,0} = -0.4, w_{2,0}=-0.6$



      I can understand that $k_{1,1} = hf_1(t_0, w_{1,0}, w_{2,0}) = hw_{2,0}$



      $f_1 = u_1'= u_2(t)$,
      So $f_1(t_0, w_{1,0}, w_{2,0}) = u_2(t_0, w_{1,0}, w_{2,0}) = w_{2,0}$ (By definition of $w_{i,j}$)



      However, I can't understand the following.
      $$k_{2,1} = hf_1(t_0 + frac{h}{2}, w_{1,0} + frac{1}{2}k_{1,1}, w_{2,0} + frac{1}{2}k_{1,2}) = hleft[w_{2,0} + frac{1}{2}k_{1,2}right]$$



      Why does $f_1(t_0 + frac{h}{2}, w_{1,0} + frac{1}{2}k_{1,1}, w_{2,0} + frac{1}{2}k_{1,2})$ equal to $w_{2,0} + frac{1}{2}k_{1,2}$? It seems that third argument in the function comes out, but there is no detailed explanation in this book.










      share|cite|improve this question







      New contributor




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











      I am studying Numerical Analysis with the book of Richard L.Burden.
      A question which I'm struggling with right now is following.



      Transform the second-order initial-value problem



      $y'' - 2y' + 2y = e^{2t}sint$ for $0 leq t leq 1, $ with $y(0) = -0.4, y'(0) = -0.6, h=0.1$



      into a system of ifrst order initial-value problems, and use the Runge-Kutta method ith h=0.1 to approximate the solution.



      Then,
      $$u_1(t) = y(t), u_2(t) = y'(t)$$
      $$u_1'(t) = u_2(t)$$
      $$u_2'(t) = e^{2t}sint - 2u_1(t) + u_2(t)$$
      $$u_1(0) = -0.4, u_2(0) = -0.6$$



      This initial conditions give $w_{1,0} = -0.4, w_{2,0}=-0.6$



      I can understand that $k_{1,1} = hf_1(t_0, w_{1,0}, w_{2,0}) = hw_{2,0}$



      $f_1 = u_1'= u_2(t)$,
      So $f_1(t_0, w_{1,0}, w_{2,0}) = u_2(t_0, w_{1,0}, w_{2,0}) = w_{2,0}$ (By definition of $w_{i,j}$)



      However, I can't understand the following.
      $$k_{2,1} = hf_1(t_0 + frac{h}{2}, w_{1,0} + frac{1}{2}k_{1,1}, w_{2,0} + frac{1}{2}k_{1,2}) = hleft[w_{2,0} + frac{1}{2}k_{1,2}right]$$



      Why does $f_1(t_0 + frac{h}{2}, w_{1,0} + frac{1}{2}k_{1,1}, w_{2,0} + frac{1}{2}k_{1,2})$ equal to $w_{2,0} + frac{1}{2}k_{1,2}$? It seems that third argument in the function comes out, but there is no detailed explanation in this book.







      runge-kutta-methods






      share|cite|improve this question







      New contributor




      James 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




      James 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






      New contributor




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









      asked 6 hours ago









      James

      1




      1




      New contributor




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





      New contributor





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






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



























          active

          oldest

          votes











          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',
          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
          });


          }
          });






          James 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%2fmath.stackexchange.com%2fquestions%2f3004621%2frunge-kutta-method-for-higher-order-differential-equations%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown






























          active

          oldest

          votes













          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes








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










           

          draft saved


          draft discarded


















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













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












          James 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%2fmath.stackexchange.com%2fquestions%2f3004621%2frunge-kutta-method-for-higher-order-differential-equations%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

          'app-layout' is not a known element: how to share Component with different Modules

          android studio warns about leanback feature tag usage required on manifest while using Unity exported app?

          WPF add header to Image with URL pettitions [duplicate]