Solving an DE involving a logarithm












2












$begingroup$


I've the following DE, describing a physical phenomenon. And the prupose is to solve that DE:



$$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$



The intial conditon is equal to $x(0)=x_0$.



For the constants (because that is maybe important for an approximation):





  • $r$ can be very large;


  • $l$ can be very large;


  • $a$ is round about $0.02526$;


  • $b$ is very small, round about $300cdot10^{-6}$;


  • $x_0$ can be very large


I've no idea where to start what so ever. Thanks for any help or ideas










share|cite|improve this question









$endgroup$












  • $begingroup$
    I think that numerical methods would be required.
    $endgroup$
    – Claude Leibovici
    Dec 30 '18 at 16:36










  • $begingroup$
    @ClaudeLeibovici I've no idea how or why
    $endgroup$
    – Klopjas
    Dec 30 '18 at 16:42










  • $begingroup$
    It's a separable equation, so it's solvable by quadratures in principle, but probably not explicitly in practice since you get a rather nasty integral...
    $endgroup$
    – Hans Lundmark
    Dec 30 '18 at 17:47


















2












$begingroup$


I've the following DE, describing a physical phenomenon. And the prupose is to solve that DE:



$$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$



The intial conditon is equal to $x(0)=x_0$.



For the constants (because that is maybe important for an approximation):





  • $r$ can be very large;


  • $l$ can be very large;


  • $a$ is round about $0.02526$;


  • $b$ is very small, round about $300cdot10^{-6}$;


  • $x_0$ can be very large


I've no idea where to start what so ever. Thanks for any help or ideas










share|cite|improve this question









$endgroup$












  • $begingroup$
    I think that numerical methods would be required.
    $endgroup$
    – Claude Leibovici
    Dec 30 '18 at 16:36










  • $begingroup$
    @ClaudeLeibovici I've no idea how or why
    $endgroup$
    – Klopjas
    Dec 30 '18 at 16:42










  • $begingroup$
    It's a separable equation, so it's solvable by quadratures in principle, but probably not explicitly in practice since you get a rather nasty integral...
    $endgroup$
    – Hans Lundmark
    Dec 30 '18 at 17:47
















2












2








2


3



$begingroup$


I've the following DE, describing a physical phenomenon. And the prupose is to solve that DE:



$$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$



The intial conditon is equal to $x(0)=x_0$.



For the constants (because that is maybe important for an approximation):





  • $r$ can be very large;


  • $l$ can be very large;


  • $a$ is round about $0.02526$;


  • $b$ is very small, round about $300cdot10^{-6}$;


  • $x_0$ can be very large


I've no idea where to start what so ever. Thanks for any help or ideas










share|cite|improve this question









$endgroup$




I've the following DE, describing a physical phenomenon. And the prupose is to solve that DE:



$$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$



The intial conditon is equal to $x(0)=x_0$.



For the constants (because that is maybe important for an approximation):





  • $r$ can be very large;


  • $l$ can be very large;


  • $a$ is round about $0.02526$;


  • $b$ is very small, round about $300cdot10^{-6}$;


  • $x_0$ can be very large


I've no idea where to start what so ever. Thanks for any help or ideas







ordinary-differential-equations functions logarithms physics mathematical-physics






share|cite|improve this question













share|cite|improve this question











share|cite|improve this question




share|cite|improve this question










asked Dec 30 '18 at 14:33









KlopjasKlopjas

714




714












  • $begingroup$
    I think that numerical methods would be required.
    $endgroup$
    – Claude Leibovici
    Dec 30 '18 at 16:36










  • $begingroup$
    @ClaudeLeibovici I've no idea how or why
    $endgroup$
    – Klopjas
    Dec 30 '18 at 16:42










  • $begingroup$
    It's a separable equation, so it's solvable by quadratures in principle, but probably not explicitly in practice since you get a rather nasty integral...
    $endgroup$
    – Hans Lundmark
    Dec 30 '18 at 17:47




















  • $begingroup$
    I think that numerical methods would be required.
    $endgroup$
    – Claude Leibovici
    Dec 30 '18 at 16:36










  • $begingroup$
    @ClaudeLeibovici I've no idea how or why
    $endgroup$
    – Klopjas
    Dec 30 '18 at 16:42










  • $begingroup$
    It's a separable equation, so it's solvable by quadratures in principle, but probably not explicitly in practice since you get a rather nasty integral...
    $endgroup$
    – Hans Lundmark
    Dec 30 '18 at 17:47


















$begingroup$
I think that numerical methods would be required.
$endgroup$
– Claude Leibovici
Dec 30 '18 at 16:36




$begingroup$
I think that numerical methods would be required.
$endgroup$
– Claude Leibovici
Dec 30 '18 at 16:36












$begingroup$
@ClaudeLeibovici I've no idea how or why
$endgroup$
– Klopjas
Dec 30 '18 at 16:42




$begingroup$
@ClaudeLeibovici I've no idea how or why
$endgroup$
– Klopjas
Dec 30 '18 at 16:42












$begingroup$
It's a separable equation, so it's solvable by quadratures in principle, but probably not explicitly in practice since you get a rather nasty integral...
$endgroup$
– Hans Lundmark
Dec 30 '18 at 17:47






$begingroup$
It's a separable equation, so it's solvable by quadratures in principle, but probably not explicitly in practice since you get a rather nasty integral...
$endgroup$
– Hans Lundmark
Dec 30 '18 at 17:47












1 Answer
1






active

oldest

votes


















2












$begingroup$

$$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$
This is a separable ODE.
$$x'(t)=frac{dx}{dt}=-frac{1}{l}left(rx(t)+alnleft(1+frac{x(t)}{b}right)right)$$
$$dt=-frac{l}{rx+alnleft(1+frac{x}{b}right)}dx$$
$$t=-lintfrac{dx}{rx+alnleft(1+frac{x}{b}right)}$$
With the condition $x(0)=x_0$ :
$$t(x)=-lint_{x_0}^xfrac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$$
As far as I know there is no closed form of this integral in terms of a finite number of standard functions. The same for the inverse function $x(t)$.



So, the analytic solution is a function defined by an integral. This is very common in practice. One have to proceed with numerical integration.



It is very easy to draw $x(t)$ thanks to usual numerical integration : Draw $t(x)$ from the above integral. Plot the points $(t,x)$ instead of $(x,t)$ , i.e. with $t$ on horizontal axis and $x$ on vertical axis.



To compute the value of $x(t)$ at a given value $t$ , proceed to the numerical integration of $-lint_{x_0}frac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$ with $xi$ increasing up to reach the specified value $t$ . The value of $xi$ at this point gives the value of $x(t)$ .






share|cite|improve this answer









$endgroup$














    Your Answer








    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%2f3056875%2fsolving-an-de-involving-a-logarithm%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









    2












    $begingroup$

    $$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$
    This is a separable ODE.
    $$x'(t)=frac{dx}{dt}=-frac{1}{l}left(rx(t)+alnleft(1+frac{x(t)}{b}right)right)$$
    $$dt=-frac{l}{rx+alnleft(1+frac{x}{b}right)}dx$$
    $$t=-lintfrac{dx}{rx+alnleft(1+frac{x}{b}right)}$$
    With the condition $x(0)=x_0$ :
    $$t(x)=-lint_{x_0}^xfrac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$$
    As far as I know there is no closed form of this integral in terms of a finite number of standard functions. The same for the inverse function $x(t)$.



    So, the analytic solution is a function defined by an integral. This is very common in practice. One have to proceed with numerical integration.



    It is very easy to draw $x(t)$ thanks to usual numerical integration : Draw $t(x)$ from the above integral. Plot the points $(t,x)$ instead of $(x,t)$ , i.e. with $t$ on horizontal axis and $x$ on vertical axis.



    To compute the value of $x(t)$ at a given value $t$ , proceed to the numerical integration of $-lint_{x_0}frac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$ with $xi$ increasing up to reach the specified value $t$ . The value of $xi$ at this point gives the value of $x(t)$ .






    share|cite|improve this answer









    $endgroup$


















      2












      $begingroup$

      $$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$
      This is a separable ODE.
      $$x'(t)=frac{dx}{dt}=-frac{1}{l}left(rx(t)+alnleft(1+frac{x(t)}{b}right)right)$$
      $$dt=-frac{l}{rx+alnleft(1+frac{x}{b}right)}dx$$
      $$t=-lintfrac{dx}{rx+alnleft(1+frac{x}{b}right)}$$
      With the condition $x(0)=x_0$ :
      $$t(x)=-lint_{x_0}^xfrac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$$
      As far as I know there is no closed form of this integral in terms of a finite number of standard functions. The same for the inverse function $x(t)$.



      So, the analytic solution is a function defined by an integral. This is very common in practice. One have to proceed with numerical integration.



      It is very easy to draw $x(t)$ thanks to usual numerical integration : Draw $t(x)$ from the above integral. Plot the points $(t,x)$ instead of $(x,t)$ , i.e. with $t$ on horizontal axis and $x$ on vertical axis.



      To compute the value of $x(t)$ at a given value $t$ , proceed to the numerical integration of $-lint_{x_0}frac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$ with $xi$ increasing up to reach the specified value $t$ . The value of $xi$ at this point gives the value of $x(t)$ .






      share|cite|improve this answer









      $endgroup$
















        2












        2








        2





        $begingroup$

        $$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$
        This is a separable ODE.
        $$x'(t)=frac{dx}{dt}=-frac{1}{l}left(rx(t)+alnleft(1+frac{x(t)}{b}right)right)$$
        $$dt=-frac{l}{rx+alnleft(1+frac{x}{b}right)}dx$$
        $$t=-lintfrac{dx}{rx+alnleft(1+frac{x}{b}right)}$$
        With the condition $x(0)=x_0$ :
        $$t(x)=-lint_{x_0}^xfrac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$$
        As far as I know there is no closed form of this integral in terms of a finite number of standard functions. The same for the inverse function $x(t)$.



        So, the analytic solution is a function defined by an integral. This is very common in practice. One have to proceed with numerical integration.



        It is very easy to draw $x(t)$ thanks to usual numerical integration : Draw $t(x)$ from the above integral. Plot the points $(t,x)$ instead of $(x,t)$ , i.e. with $t$ on horizontal axis and $x$ on vertical axis.



        To compute the value of $x(t)$ at a given value $t$ , proceed to the numerical integration of $-lint_{x_0}frac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$ with $xi$ increasing up to reach the specified value $t$ . The value of $xi$ at this point gives the value of $x(t)$ .






        share|cite|improve this answer









        $endgroup$



        $$x(t)cdot r+x'(t)cdot l+acdotlnleft(1+frac{x(t)}{b}right)=0spaceLongleftrightarrowspace x(t)=dots$$
        This is a separable ODE.
        $$x'(t)=frac{dx}{dt}=-frac{1}{l}left(rx(t)+alnleft(1+frac{x(t)}{b}right)right)$$
        $$dt=-frac{l}{rx+alnleft(1+frac{x}{b}right)}dx$$
        $$t=-lintfrac{dx}{rx+alnleft(1+frac{x}{b}right)}$$
        With the condition $x(0)=x_0$ :
        $$t(x)=-lint_{x_0}^xfrac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$$
        As far as I know there is no closed form of this integral in terms of a finite number of standard functions. The same for the inverse function $x(t)$.



        So, the analytic solution is a function defined by an integral. This is very common in practice. One have to proceed with numerical integration.



        It is very easy to draw $x(t)$ thanks to usual numerical integration : Draw $t(x)$ from the above integral. Plot the points $(t,x)$ instead of $(x,t)$ , i.e. with $t$ on horizontal axis and $x$ on vertical axis.



        To compute the value of $x(t)$ at a given value $t$ , proceed to the numerical integration of $-lint_{x_0}frac{dxi}{rxi+alnleft(1+frac{xi}{b}right)}$ with $xi$ increasing up to reach the specified value $t$ . The value of $xi$ at this point gives the value of $x(t)$ .







        share|cite|improve this answer












        share|cite|improve this answer



        share|cite|improve this answer










        answered Dec 30 '18 at 18:08









        JJacquelinJJacquelin

        45.8k21858




        45.8k21858






























            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%2f3056875%2fsolving-an-de-involving-a-logarithm%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

            How to change which sound is reproduced for terminal bell?

            Can I use Tabulator js library in my java Spring + Thymeleaf project?

            Title Spacing in Bjornstrup Chapter, Removing Chapter Number From Contents