Differential equations with a parametric form
$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
calculus integration ordinary-differential-equations parametric
$endgroup$
add a comment |
$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
calculus integration ordinary-differential-equations parametric
$endgroup$
add a comment |
$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
calculus integration ordinary-differential-equations parametric
$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
calculus integration ordinary-differential-equations parametric
asked Jan 9 at 17:27
nipohc88nipohc88
163
163
add a comment |
add a comment |
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
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
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.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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