How to solve the coupled differential equation












0












$begingroup$


I have the following differential equations.



$$ tau Vfrac{partial r}{partial x} + Wfrac{partial ^2r}{partial x^2} + [r - Cu(1-r^2)][1-r^2]=0,$$



$$ Vfrac{partial u}{partial x} + Dfrac{partial^2 u}{partial x^2} - frac{V}{2}frac{partial r}{partial x}=0,$$ with far field boundary condition as $u_{x=infty} = 0.6.$ and $u_{x=0} = 0$. Here, $tau, V, W, D, C$ are dimensionless constants.



Now, I would like to solve the second equation subject to boundary conditions. But the last term in the second equation confuses me on how to proceed. Can anyone help me out with this ?



thanks










share|cite|improve this question











$endgroup$












  • $begingroup$
    I would rather say that the last term in the first equation is the most difficult to deal with. Is $x$ the only independent variable? Why are those partial differential equations instead of ordinary equations? There are a handful numerical methods that would help you to solve these equations if they are ordinary.
    $endgroup$
    – rafa11111
    Dec 18 '18 at 12:51










  • $begingroup$
    @rafa11111 Ok. let us say these are ordinary equations. I would like to solve it analytically.
    $endgroup$
    – newstudent
    Dec 18 '18 at 12:52










  • $begingroup$
    If you are sure these are ordinary equations, please edit the question to address so (change the equation, the tag, etc.). I don't think there is an analytical solution to these equations, mainly due to the nonlinear term in the first equation. If $C$ is a small number, you can use a perturbation method, that provides an approximate analytical solution.
    $endgroup$
    – rafa11111
    Dec 18 '18 at 13:01










  • $begingroup$
    You can integrate the second equation once then substitute for $r(x)$ in the first to give an equation in $u(x)$ alone. I think your journey will end there!
    $endgroup$
    – user121049
    Dec 19 '18 at 8:00










  • $begingroup$
    You need some more boundary conditions on $r$. The only hope I see is to assume $r$ stays small so you can set $1-r^2approx 1$ and linearize the first equation. The two equations are then easy to solve. Then check $r$ really stays small.
    $endgroup$
    – user121049
    Dec 20 '18 at 9:00
















0












$begingroup$


I have the following differential equations.



$$ tau Vfrac{partial r}{partial x} + Wfrac{partial ^2r}{partial x^2} + [r - Cu(1-r^2)][1-r^2]=0,$$



$$ Vfrac{partial u}{partial x} + Dfrac{partial^2 u}{partial x^2} - frac{V}{2}frac{partial r}{partial x}=0,$$ with far field boundary condition as $u_{x=infty} = 0.6.$ and $u_{x=0} = 0$. Here, $tau, V, W, D, C$ are dimensionless constants.



Now, I would like to solve the second equation subject to boundary conditions. But the last term in the second equation confuses me on how to proceed. Can anyone help me out with this ?



thanks










share|cite|improve this question











$endgroup$












  • $begingroup$
    I would rather say that the last term in the first equation is the most difficult to deal with. Is $x$ the only independent variable? Why are those partial differential equations instead of ordinary equations? There are a handful numerical methods that would help you to solve these equations if they are ordinary.
    $endgroup$
    – rafa11111
    Dec 18 '18 at 12:51










  • $begingroup$
    @rafa11111 Ok. let us say these are ordinary equations. I would like to solve it analytically.
    $endgroup$
    – newstudent
    Dec 18 '18 at 12:52










  • $begingroup$
    If you are sure these are ordinary equations, please edit the question to address so (change the equation, the tag, etc.). I don't think there is an analytical solution to these equations, mainly due to the nonlinear term in the first equation. If $C$ is a small number, you can use a perturbation method, that provides an approximate analytical solution.
    $endgroup$
    – rafa11111
    Dec 18 '18 at 13:01










  • $begingroup$
    You can integrate the second equation once then substitute for $r(x)$ in the first to give an equation in $u(x)$ alone. I think your journey will end there!
    $endgroup$
    – user121049
    Dec 19 '18 at 8:00










  • $begingroup$
    You need some more boundary conditions on $r$. The only hope I see is to assume $r$ stays small so you can set $1-r^2approx 1$ and linearize the first equation. The two equations are then easy to solve. Then check $r$ really stays small.
    $endgroup$
    – user121049
    Dec 20 '18 at 9:00














0












0








0





$begingroup$


I have the following differential equations.



$$ tau Vfrac{partial r}{partial x} + Wfrac{partial ^2r}{partial x^2} + [r - Cu(1-r^2)][1-r^2]=0,$$



$$ Vfrac{partial u}{partial x} + Dfrac{partial^2 u}{partial x^2} - frac{V}{2}frac{partial r}{partial x}=0,$$ with far field boundary condition as $u_{x=infty} = 0.6.$ and $u_{x=0} = 0$. Here, $tau, V, W, D, C$ are dimensionless constants.



Now, I would like to solve the second equation subject to boundary conditions. But the last term in the second equation confuses me on how to proceed. Can anyone help me out with this ?



thanks










share|cite|improve this question











$endgroup$




I have the following differential equations.



$$ tau Vfrac{partial r}{partial x} + Wfrac{partial ^2r}{partial x^2} + [r - Cu(1-r^2)][1-r^2]=0,$$



$$ Vfrac{partial u}{partial x} + Dfrac{partial^2 u}{partial x^2} - frac{V}{2}frac{partial r}{partial x}=0,$$ with far field boundary condition as $u_{x=infty} = 0.6.$ and $u_{x=0} = 0$. Here, $tau, V, W, D, C$ are dimensionless constants.



Now, I would like to solve the second equation subject to boundary conditions. But the last term in the second equation confuses me on how to proceed. Can anyone help me out with this ?



thanks







ordinary-differential-equations numerical-methods






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Dec 18 '18 at 13:03







newstudent

















asked Dec 18 '18 at 12:38









newstudentnewstudent

366




366












  • $begingroup$
    I would rather say that the last term in the first equation is the most difficult to deal with. Is $x$ the only independent variable? Why are those partial differential equations instead of ordinary equations? There are a handful numerical methods that would help you to solve these equations if they are ordinary.
    $endgroup$
    – rafa11111
    Dec 18 '18 at 12:51










  • $begingroup$
    @rafa11111 Ok. let us say these are ordinary equations. I would like to solve it analytically.
    $endgroup$
    – newstudent
    Dec 18 '18 at 12:52










  • $begingroup$
    If you are sure these are ordinary equations, please edit the question to address so (change the equation, the tag, etc.). I don't think there is an analytical solution to these equations, mainly due to the nonlinear term in the first equation. If $C$ is a small number, you can use a perturbation method, that provides an approximate analytical solution.
    $endgroup$
    – rafa11111
    Dec 18 '18 at 13:01










  • $begingroup$
    You can integrate the second equation once then substitute for $r(x)$ in the first to give an equation in $u(x)$ alone. I think your journey will end there!
    $endgroup$
    – user121049
    Dec 19 '18 at 8:00










  • $begingroup$
    You need some more boundary conditions on $r$. The only hope I see is to assume $r$ stays small so you can set $1-r^2approx 1$ and linearize the first equation. The two equations are then easy to solve. Then check $r$ really stays small.
    $endgroup$
    – user121049
    Dec 20 '18 at 9:00


















  • $begingroup$
    I would rather say that the last term in the first equation is the most difficult to deal with. Is $x$ the only independent variable? Why are those partial differential equations instead of ordinary equations? There are a handful numerical methods that would help you to solve these equations if they are ordinary.
    $endgroup$
    – rafa11111
    Dec 18 '18 at 12:51










  • $begingroup$
    @rafa11111 Ok. let us say these are ordinary equations. I would like to solve it analytically.
    $endgroup$
    – newstudent
    Dec 18 '18 at 12:52










  • $begingroup$
    If you are sure these are ordinary equations, please edit the question to address so (change the equation, the tag, etc.). I don't think there is an analytical solution to these equations, mainly due to the nonlinear term in the first equation. If $C$ is a small number, you can use a perturbation method, that provides an approximate analytical solution.
    $endgroup$
    – rafa11111
    Dec 18 '18 at 13:01










  • $begingroup$
    You can integrate the second equation once then substitute for $r(x)$ in the first to give an equation in $u(x)$ alone. I think your journey will end there!
    $endgroup$
    – user121049
    Dec 19 '18 at 8:00










  • $begingroup$
    You need some more boundary conditions on $r$. The only hope I see is to assume $r$ stays small so you can set $1-r^2approx 1$ and linearize the first equation. The two equations are then easy to solve. Then check $r$ really stays small.
    $endgroup$
    – user121049
    Dec 20 '18 at 9:00
















$begingroup$
I would rather say that the last term in the first equation is the most difficult to deal with. Is $x$ the only independent variable? Why are those partial differential equations instead of ordinary equations? There are a handful numerical methods that would help you to solve these equations if they are ordinary.
$endgroup$
– rafa11111
Dec 18 '18 at 12:51




$begingroup$
I would rather say that the last term in the first equation is the most difficult to deal with. Is $x$ the only independent variable? Why are those partial differential equations instead of ordinary equations? There are a handful numerical methods that would help you to solve these equations if they are ordinary.
$endgroup$
– rafa11111
Dec 18 '18 at 12:51












$begingroup$
@rafa11111 Ok. let us say these are ordinary equations. I would like to solve it analytically.
$endgroup$
– newstudent
Dec 18 '18 at 12:52




$begingroup$
@rafa11111 Ok. let us say these are ordinary equations. I would like to solve it analytically.
$endgroup$
– newstudent
Dec 18 '18 at 12:52












$begingroup$
If you are sure these are ordinary equations, please edit the question to address so (change the equation, the tag, etc.). I don't think there is an analytical solution to these equations, mainly due to the nonlinear term in the first equation. If $C$ is a small number, you can use a perturbation method, that provides an approximate analytical solution.
$endgroup$
– rafa11111
Dec 18 '18 at 13:01




$begingroup$
If you are sure these are ordinary equations, please edit the question to address so (change the equation, the tag, etc.). I don't think there is an analytical solution to these equations, mainly due to the nonlinear term in the first equation. If $C$ is a small number, you can use a perturbation method, that provides an approximate analytical solution.
$endgroup$
– rafa11111
Dec 18 '18 at 13:01












$begingroup$
You can integrate the second equation once then substitute for $r(x)$ in the first to give an equation in $u(x)$ alone. I think your journey will end there!
$endgroup$
– user121049
Dec 19 '18 at 8:00




$begingroup$
You can integrate the second equation once then substitute for $r(x)$ in the first to give an equation in $u(x)$ alone. I think your journey will end there!
$endgroup$
– user121049
Dec 19 '18 at 8:00












$begingroup$
You need some more boundary conditions on $r$. The only hope I see is to assume $r$ stays small so you can set $1-r^2approx 1$ and linearize the first equation. The two equations are then easy to solve. Then check $r$ really stays small.
$endgroup$
– user121049
Dec 20 '18 at 9:00




$begingroup$
You need some more boundary conditions on $r$. The only hope I see is to assume $r$ stays small so you can set $1-r^2approx 1$ and linearize the first equation. The two equations are then easy to solve. Then check $r$ really stays small.
$endgroup$
– user121049
Dec 20 '18 at 9:00










1 Answer
1






active

oldest

votes


















1












$begingroup$

Disclaimer: I've just realized you can't actually solve for $r(x)$ in the first equation, which makes my answer kind of useless. I'll leave it up anyway, to learn from my mistakes.



Suppose you've already solved the first equation, and $r(x)$ is known, then the second equation is linear and non-homogeneous, and can easily be solved using variations of parameter



Let $a = frac{V}{D}$



$$ frac{d^2u}{dx^2} + afrac{du}{dx} = frac{a}{2}frac{dr}{dx} $$



The homogeneous equation has solution $u_h = c_1 + c_2e^{-ax}$, so guess a particular solution of the form



$$ u_p(x) = p(x) + q(x)e^{-ax} $$



Then



begin{align} p'(x) + q'(x)e^{-ax} &= 0 \ -aq'(x)e^{-ax} &= frac{a}{2}r'(x) end{align}



or $p'(x) = frac{1}{2}r'(x)$ and $q'(x) = -frac{1}{2}r'(x)e^{-ax}$



You can integrate to find



begin{align} p(x) &= frac{1}{2}r(x) \ q(x) &=-frac12int r'(x) e^{ax} = -frac12 r(x)e^{ax} dx + frac{a}{2}int r(x)e^{ax} dx end{align}



Therefore



$$ u_p(x) = frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



Matching the boundary conditions, we find



$$ u(x) = u_{infty}(1-e^{-ax}) + frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



where $u_infty$ is some constant depending on the behavior of $r(x)$ at ${xtoinfty}$






share|cite|improve this answer











$endgroup$













  • $begingroup$
    There is a $u(x)$ in the first equation so how do you solve it for $r(x)$?
    $endgroup$
    – user121049
    Dec 19 '18 at 20:37










  • $begingroup$
    @user121049 Wow I didn't even notice. Thanks for pointing it out.
    $endgroup$
    – Dylan
    Dec 20 '18 at 3:57











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%2f3045117%2fhow-to-solve-the-coupled-differential-equation%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









1












$begingroup$

Disclaimer: I've just realized you can't actually solve for $r(x)$ in the first equation, which makes my answer kind of useless. I'll leave it up anyway, to learn from my mistakes.



Suppose you've already solved the first equation, and $r(x)$ is known, then the second equation is linear and non-homogeneous, and can easily be solved using variations of parameter



Let $a = frac{V}{D}$



$$ frac{d^2u}{dx^2} + afrac{du}{dx} = frac{a}{2}frac{dr}{dx} $$



The homogeneous equation has solution $u_h = c_1 + c_2e^{-ax}$, so guess a particular solution of the form



$$ u_p(x) = p(x) + q(x)e^{-ax} $$



Then



begin{align} p'(x) + q'(x)e^{-ax} &= 0 \ -aq'(x)e^{-ax} &= frac{a}{2}r'(x) end{align}



or $p'(x) = frac{1}{2}r'(x)$ and $q'(x) = -frac{1}{2}r'(x)e^{-ax}$



You can integrate to find



begin{align} p(x) &= frac{1}{2}r(x) \ q(x) &=-frac12int r'(x) e^{ax} = -frac12 r(x)e^{ax} dx + frac{a}{2}int r(x)e^{ax} dx end{align}



Therefore



$$ u_p(x) = frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



Matching the boundary conditions, we find



$$ u(x) = u_{infty}(1-e^{-ax}) + frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



where $u_infty$ is some constant depending on the behavior of $r(x)$ at ${xtoinfty}$






share|cite|improve this answer











$endgroup$













  • $begingroup$
    There is a $u(x)$ in the first equation so how do you solve it for $r(x)$?
    $endgroup$
    – user121049
    Dec 19 '18 at 20:37










  • $begingroup$
    @user121049 Wow I didn't even notice. Thanks for pointing it out.
    $endgroup$
    – Dylan
    Dec 20 '18 at 3:57
















1












$begingroup$

Disclaimer: I've just realized you can't actually solve for $r(x)$ in the first equation, which makes my answer kind of useless. I'll leave it up anyway, to learn from my mistakes.



Suppose you've already solved the first equation, and $r(x)$ is known, then the second equation is linear and non-homogeneous, and can easily be solved using variations of parameter



Let $a = frac{V}{D}$



$$ frac{d^2u}{dx^2} + afrac{du}{dx} = frac{a}{2}frac{dr}{dx} $$



The homogeneous equation has solution $u_h = c_1 + c_2e^{-ax}$, so guess a particular solution of the form



$$ u_p(x) = p(x) + q(x)e^{-ax} $$



Then



begin{align} p'(x) + q'(x)e^{-ax} &= 0 \ -aq'(x)e^{-ax} &= frac{a}{2}r'(x) end{align}



or $p'(x) = frac{1}{2}r'(x)$ and $q'(x) = -frac{1}{2}r'(x)e^{-ax}$



You can integrate to find



begin{align} p(x) &= frac{1}{2}r(x) \ q(x) &=-frac12int r'(x) e^{ax} = -frac12 r(x)e^{ax} dx + frac{a}{2}int r(x)e^{ax} dx end{align}



Therefore



$$ u_p(x) = frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



Matching the boundary conditions, we find



$$ u(x) = u_{infty}(1-e^{-ax}) + frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



where $u_infty$ is some constant depending on the behavior of $r(x)$ at ${xtoinfty}$






share|cite|improve this answer











$endgroup$













  • $begingroup$
    There is a $u(x)$ in the first equation so how do you solve it for $r(x)$?
    $endgroup$
    – user121049
    Dec 19 '18 at 20:37










  • $begingroup$
    @user121049 Wow I didn't even notice. Thanks for pointing it out.
    $endgroup$
    – Dylan
    Dec 20 '18 at 3:57














1












1








1





$begingroup$

Disclaimer: I've just realized you can't actually solve for $r(x)$ in the first equation, which makes my answer kind of useless. I'll leave it up anyway, to learn from my mistakes.



Suppose you've already solved the first equation, and $r(x)$ is known, then the second equation is linear and non-homogeneous, and can easily be solved using variations of parameter



Let $a = frac{V}{D}$



$$ frac{d^2u}{dx^2} + afrac{du}{dx} = frac{a}{2}frac{dr}{dx} $$



The homogeneous equation has solution $u_h = c_1 + c_2e^{-ax}$, so guess a particular solution of the form



$$ u_p(x) = p(x) + q(x)e^{-ax} $$



Then



begin{align} p'(x) + q'(x)e^{-ax} &= 0 \ -aq'(x)e^{-ax} &= frac{a}{2}r'(x) end{align}



or $p'(x) = frac{1}{2}r'(x)$ and $q'(x) = -frac{1}{2}r'(x)e^{-ax}$



You can integrate to find



begin{align} p(x) &= frac{1}{2}r(x) \ q(x) &=-frac12int r'(x) e^{ax} = -frac12 r(x)e^{ax} dx + frac{a}{2}int r(x)e^{ax} dx end{align}



Therefore



$$ u_p(x) = frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



Matching the boundary conditions, we find



$$ u(x) = u_{infty}(1-e^{-ax}) + frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



where $u_infty$ is some constant depending on the behavior of $r(x)$ at ${xtoinfty}$






share|cite|improve this answer











$endgroup$



Disclaimer: I've just realized you can't actually solve for $r(x)$ in the first equation, which makes my answer kind of useless. I'll leave it up anyway, to learn from my mistakes.



Suppose you've already solved the first equation, and $r(x)$ is known, then the second equation is linear and non-homogeneous, and can easily be solved using variations of parameter



Let $a = frac{V}{D}$



$$ frac{d^2u}{dx^2} + afrac{du}{dx} = frac{a}{2}frac{dr}{dx} $$



The homogeneous equation has solution $u_h = c_1 + c_2e^{-ax}$, so guess a particular solution of the form



$$ u_p(x) = p(x) + q(x)e^{-ax} $$



Then



begin{align} p'(x) + q'(x)e^{-ax} &= 0 \ -aq'(x)e^{-ax} &= frac{a}{2}r'(x) end{align}



or $p'(x) = frac{1}{2}r'(x)$ and $q'(x) = -frac{1}{2}r'(x)e^{-ax}$



You can integrate to find



begin{align} p(x) &= frac{1}{2}r(x) \ q(x) &=-frac12int r'(x) e^{ax} = -frac12 r(x)e^{ax} dx + frac{a}{2}int r(x)e^{ax} dx end{align}



Therefore



$$ u_p(x) = frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



Matching the boundary conditions, we find



$$ u(x) = u_{infty}(1-e^{-ax}) + frac{a}{2}e^{-ax}int_0^x r(t)e^{at} dt $$



where $u_infty$ is some constant depending on the behavior of $r(x)$ at ${xtoinfty}$







share|cite|improve this answer














share|cite|improve this answer



share|cite|improve this answer








edited Dec 20 '18 at 4:03

























answered Dec 19 '18 at 7:48









DylanDylan

14.1k31127




14.1k31127












  • $begingroup$
    There is a $u(x)$ in the first equation so how do you solve it for $r(x)$?
    $endgroup$
    – user121049
    Dec 19 '18 at 20:37










  • $begingroup$
    @user121049 Wow I didn't even notice. Thanks for pointing it out.
    $endgroup$
    – Dylan
    Dec 20 '18 at 3:57


















  • $begingroup$
    There is a $u(x)$ in the first equation so how do you solve it for $r(x)$?
    $endgroup$
    – user121049
    Dec 19 '18 at 20:37










  • $begingroup$
    @user121049 Wow I didn't even notice. Thanks for pointing it out.
    $endgroup$
    – Dylan
    Dec 20 '18 at 3:57
















$begingroup$
There is a $u(x)$ in the first equation so how do you solve it for $r(x)$?
$endgroup$
– user121049
Dec 19 '18 at 20:37




$begingroup$
There is a $u(x)$ in the first equation so how do you solve it for $r(x)$?
$endgroup$
– user121049
Dec 19 '18 at 20:37












$begingroup$
@user121049 Wow I didn't even notice. Thanks for pointing it out.
$endgroup$
– Dylan
Dec 20 '18 at 3:57




$begingroup$
@user121049 Wow I didn't even notice. Thanks for pointing it out.
$endgroup$
– Dylan
Dec 20 '18 at 3:57


















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%2f3045117%2fhow-to-solve-the-coupled-differential-equation%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

Plaza Victoria

Puebla de Zaragoza

Musa