Differential equations with a parametric form












0












$begingroup$


I have been exploring mathematical epidemiology for a school project and have decided to model epidemics using the SIR model (Susceptibles, Infected and Recovered), so $x(t), y(t), z(t)$ are functions describing the number of individuals belonging to each group. The model is described by a system of differential equations:
begin{equation}
frac{dx}{dt}=-beta x(t)y(t)end{equation}

begin{equation}
frac{dy}{dt}=beta x(t)y(t) - gamma y(t)
end{equation}

begin{equation}
frac{dz}{dt} = gamma y(t)
end{equation}

Where, $ beta, gamma $ are infection and recovery rates.
I have been reading a research paper showing how the system can be solved https://arxiv.org/abs/1403.2160, however, I have some difficulty with it.
(I'm only in the $10^{th}$ grade and have little understanding of differential equations}
This system can be reduced to an equation(I understand how):
begin{equation}
z''=-x_{0}beta z'e^{-frac{beta}{gamma}z}-gamma z
end{equation}

Where $x_{0}$ is the susceptible population at $t=0$. Then in the paper they introduce a function
begin{equation}
u(t)=e^{-frac{beta}{gamma}z(t)}
end{equation}

Which at $t_0$ takes the value
begin{equation}
u_{0}=e^{-frac{beta}{gamma}z}
end{equation}

Substitution of this equation into the previous one leads to:
begin{equation}
ufrac{d^{2}u}{dt^{2}}-bigg(frac{du}{dt}bigg)^{2}+bigg(gamma - x_{0}beta ubigg)u frac{du}{dt} = 0
end{equation}

Then they introduce an equation:
begin{equation}
phi=frac{dt}{du}
end{equation}

By substituting it into the previous equation they get:
begin{equation}
frac{dphi}{du}+frac{1}{u} phi = (gamma-x_{0}beta u)phi^{2}
end{equation}

The solution of this equation is:
begin{equation}
phi = frac{1}{u(C_{1}-gamma ln u +x_{0}beta u)}
end{equation}

Where $C_{1}$ is an arbitrary integration constant. Up to this point I have no problems, however, then they proceed with saying in view of equation $phi=frac{dt}{du}$ and the previous equation, they get:
begin{equation}
t-t_{0} = int_{u_0}^{u}frac{dxi}{xi(C_{1}-gamma lnxi + x_{0}beta xi}
end{equation}

where one can choose $t_{0}=0$ without loss of generality. Then they proceed by saying hence we have obtained the complete exact solution of the original system, given in parametric form by:
begin{equation}
x=x_{0}u,
end{equation}

begin{equation}
y=frac{gamma}{beta}ln u -x_{0}u - frac{C_{1}}{beta},
end{equation}

begin{equation}
z = -frac{gamma}{beta}ln u
end{equation}

and the next equation (which is derived from the fact that $x+y+z=N$ where $N$ is the whole population which stays constant during the epidemic) understand is:
begin{equation}
C_{1} = -beta N
end{equation}

Could someone please explain how they acquired the last few equations and what they mean, moreover, if I wanted to graph them how could I do it? Thank you










share|cite|improve this question









$endgroup$

















    0












    $begingroup$


    I have been exploring mathematical epidemiology for a school project and have decided to model epidemics using the SIR model (Susceptibles, Infected and Recovered), so $x(t), y(t), z(t)$ are functions describing the number of individuals belonging to each group. The model is described by a system of differential equations:
    begin{equation}
    frac{dx}{dt}=-beta x(t)y(t)end{equation}

    begin{equation}
    frac{dy}{dt}=beta x(t)y(t) - gamma y(t)
    end{equation}

    begin{equation}
    frac{dz}{dt} = gamma y(t)
    end{equation}

    Where, $ beta, gamma $ are infection and recovery rates.
    I have been reading a research paper showing how the system can be solved https://arxiv.org/abs/1403.2160, however, I have some difficulty with it.
    (I'm only in the $10^{th}$ grade and have little understanding of differential equations}
    This system can be reduced to an equation(I understand how):
    begin{equation}
    z''=-x_{0}beta z'e^{-frac{beta}{gamma}z}-gamma z
    end{equation}

    Where $x_{0}$ is the susceptible population at $t=0$. Then in the paper they introduce a function
    begin{equation}
    u(t)=e^{-frac{beta}{gamma}z(t)}
    end{equation}

    Which at $t_0$ takes the value
    begin{equation}
    u_{0}=e^{-frac{beta}{gamma}z}
    end{equation}

    Substitution of this equation into the previous one leads to:
    begin{equation}
    ufrac{d^{2}u}{dt^{2}}-bigg(frac{du}{dt}bigg)^{2}+bigg(gamma - x_{0}beta ubigg)u frac{du}{dt} = 0
    end{equation}

    Then they introduce an equation:
    begin{equation}
    phi=frac{dt}{du}
    end{equation}

    By substituting it into the previous equation they get:
    begin{equation}
    frac{dphi}{du}+frac{1}{u} phi = (gamma-x_{0}beta u)phi^{2}
    end{equation}

    The solution of this equation is:
    begin{equation}
    phi = frac{1}{u(C_{1}-gamma ln u +x_{0}beta u)}
    end{equation}

    Where $C_{1}$ is an arbitrary integration constant. Up to this point I have no problems, however, then they proceed with saying in view of equation $phi=frac{dt}{du}$ and the previous equation, they get:
    begin{equation}
    t-t_{0} = int_{u_0}^{u}frac{dxi}{xi(C_{1}-gamma lnxi + x_{0}beta xi}
    end{equation}

    where one can choose $t_{0}=0$ without loss of generality. Then they proceed by saying hence we have obtained the complete exact solution of the original system, given in parametric form by:
    begin{equation}
    x=x_{0}u,
    end{equation}

    begin{equation}
    y=frac{gamma}{beta}ln u -x_{0}u - frac{C_{1}}{beta},
    end{equation}

    begin{equation}
    z = -frac{gamma}{beta}ln u
    end{equation}

    and the next equation (which is derived from the fact that $x+y+z=N$ where $N$ is the whole population which stays constant during the epidemic) understand is:
    begin{equation}
    C_{1} = -beta N
    end{equation}

    Could someone please explain how they acquired the last few equations and what they mean, moreover, if I wanted to graph them how could I do it? Thank you










    share|cite|improve this question









    $endgroup$















      0












      0








      0





      $begingroup$


      I have been exploring mathematical epidemiology for a school project and have decided to model epidemics using the SIR model (Susceptibles, Infected and Recovered), so $x(t), y(t), z(t)$ are functions describing the number of individuals belonging to each group. The model is described by a system of differential equations:
      begin{equation}
      frac{dx}{dt}=-beta x(t)y(t)end{equation}

      begin{equation}
      frac{dy}{dt}=beta x(t)y(t) - gamma y(t)
      end{equation}

      begin{equation}
      frac{dz}{dt} = gamma y(t)
      end{equation}

      Where, $ beta, gamma $ are infection and recovery rates.
      I have been reading a research paper showing how the system can be solved https://arxiv.org/abs/1403.2160, however, I have some difficulty with it.
      (I'm only in the $10^{th}$ grade and have little understanding of differential equations}
      This system can be reduced to an equation(I understand how):
      begin{equation}
      z''=-x_{0}beta z'e^{-frac{beta}{gamma}z}-gamma z
      end{equation}

      Where $x_{0}$ is the susceptible population at $t=0$. Then in the paper they introduce a function
      begin{equation}
      u(t)=e^{-frac{beta}{gamma}z(t)}
      end{equation}

      Which at $t_0$ takes the value
      begin{equation}
      u_{0}=e^{-frac{beta}{gamma}z}
      end{equation}

      Substitution of this equation into the previous one leads to:
      begin{equation}
      ufrac{d^{2}u}{dt^{2}}-bigg(frac{du}{dt}bigg)^{2}+bigg(gamma - x_{0}beta ubigg)u frac{du}{dt} = 0
      end{equation}

      Then they introduce an equation:
      begin{equation}
      phi=frac{dt}{du}
      end{equation}

      By substituting it into the previous equation they get:
      begin{equation}
      frac{dphi}{du}+frac{1}{u} phi = (gamma-x_{0}beta u)phi^{2}
      end{equation}

      The solution of this equation is:
      begin{equation}
      phi = frac{1}{u(C_{1}-gamma ln u +x_{0}beta u)}
      end{equation}

      Where $C_{1}$ is an arbitrary integration constant. Up to this point I have no problems, however, then they proceed with saying in view of equation $phi=frac{dt}{du}$ and the previous equation, they get:
      begin{equation}
      t-t_{0} = int_{u_0}^{u}frac{dxi}{xi(C_{1}-gamma lnxi + x_{0}beta xi}
      end{equation}

      where one can choose $t_{0}=0$ without loss of generality. Then they proceed by saying hence we have obtained the complete exact solution of the original system, given in parametric form by:
      begin{equation}
      x=x_{0}u,
      end{equation}

      begin{equation}
      y=frac{gamma}{beta}ln u -x_{0}u - frac{C_{1}}{beta},
      end{equation}

      begin{equation}
      z = -frac{gamma}{beta}ln u
      end{equation}

      and the next equation (which is derived from the fact that $x+y+z=N$ where $N$ is the whole population which stays constant during the epidemic) understand is:
      begin{equation}
      C_{1} = -beta N
      end{equation}

      Could someone please explain how they acquired the last few equations and what they mean, moreover, if I wanted to graph them how could I do it? Thank you










      share|cite|improve this question









      $endgroup$




      I have been exploring mathematical epidemiology for a school project and have decided to model epidemics using the SIR model (Susceptibles, Infected and Recovered), so $x(t), y(t), z(t)$ are functions describing the number of individuals belonging to each group. The model is described by a system of differential equations:
      begin{equation}
      frac{dx}{dt}=-beta x(t)y(t)end{equation}

      begin{equation}
      frac{dy}{dt}=beta x(t)y(t) - gamma y(t)
      end{equation}

      begin{equation}
      frac{dz}{dt} = gamma y(t)
      end{equation}

      Where, $ beta, gamma $ are infection and recovery rates.
      I have been reading a research paper showing how the system can be solved https://arxiv.org/abs/1403.2160, however, I have some difficulty with it.
      (I'm only in the $10^{th}$ grade and have little understanding of differential equations}
      This system can be reduced to an equation(I understand how):
      begin{equation}
      z''=-x_{0}beta z'e^{-frac{beta}{gamma}z}-gamma z
      end{equation}

      Where $x_{0}$ is the susceptible population at $t=0$. Then in the paper they introduce a function
      begin{equation}
      u(t)=e^{-frac{beta}{gamma}z(t)}
      end{equation}

      Which at $t_0$ takes the value
      begin{equation}
      u_{0}=e^{-frac{beta}{gamma}z}
      end{equation}

      Substitution of this equation into the previous one leads to:
      begin{equation}
      ufrac{d^{2}u}{dt^{2}}-bigg(frac{du}{dt}bigg)^{2}+bigg(gamma - x_{0}beta ubigg)u frac{du}{dt} = 0
      end{equation}

      Then they introduce an equation:
      begin{equation}
      phi=frac{dt}{du}
      end{equation}

      By substituting it into the previous equation they get:
      begin{equation}
      frac{dphi}{du}+frac{1}{u} phi = (gamma-x_{0}beta u)phi^{2}
      end{equation}

      The solution of this equation is:
      begin{equation}
      phi = frac{1}{u(C_{1}-gamma ln u +x_{0}beta u)}
      end{equation}

      Where $C_{1}$ is an arbitrary integration constant. Up to this point I have no problems, however, then they proceed with saying in view of equation $phi=frac{dt}{du}$ and the previous equation, they get:
      begin{equation}
      t-t_{0} = int_{u_0}^{u}frac{dxi}{xi(C_{1}-gamma lnxi + x_{0}beta xi}
      end{equation}

      where one can choose $t_{0}=0$ without loss of generality. Then they proceed by saying hence we have obtained the complete exact solution of the original system, given in parametric form by:
      begin{equation}
      x=x_{0}u,
      end{equation}

      begin{equation}
      y=frac{gamma}{beta}ln u -x_{0}u - frac{C_{1}}{beta},
      end{equation}

      begin{equation}
      z = -frac{gamma}{beta}ln u
      end{equation}

      and the next equation (which is derived from the fact that $x+y+z=N$ where $N$ is the whole population which stays constant during the epidemic) understand is:
      begin{equation}
      C_{1} = -beta N
      end{equation}

      Could someone please explain how they acquired the last few equations and what they mean, moreover, if I wanted to graph them how could I do it? Thank you







      calculus integration ordinary-differential-equations parametric






      share|cite|improve this question













      share|cite|improve this question











      share|cite|improve this question




      share|cite|improve this question










      asked Jan 9 at 17:27









      nipohc88nipohc88

      163




      163






















          0






          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',
          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%2f3067740%2fdifferential-equations-with-a-parametric-form%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown

























          0






          active

          oldest

          votes








          0






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes
















          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.




          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3067740%2fdifferential-equations-with-a-parametric-form%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

          MongoDB - Not Authorized To Execute Command

          How to fix TextFormField cause rebuild widget in Flutter

          in spring boot 2.1 many test slices are not allowed anymore due to multiple @BootstrapWith