Green function for solving ODE: continuity and discontinuity of derivatives












1












$begingroup$


I am studying Green functions and I want to find the Green function $G(x,x_0)$ of the linear operator
begin{equation}
mathcal{L}=left(frac{d^n}{dx^n}+a_1(x)frac{d^{n-1}}{dx^{n-1}}+dots+a_n(x) right)
end{equation}

This means solving the following problem
begin{equation}
frac{d^n G}{dx^n}+a_1(x)frac{d^{n-1}G}{dx^{n-1}}+dots+a_n(x)G=delta(x-x_0)
tag{1}
end{equation}



I understand that I have to impose continuity in $x_0$ for all the derivatives, up to $n-2$, and impose the jump discontinuity in $x_0$ for the $(n-1)$-th derivative. I don't exactly get why this is the case. I have searched in some lecture notes and found the following clues:



First approach



They proceed by integrating equation $(1)$, and they say "we obtain$G^{(n-1)}=H(x-x_0)$ + some continuous functions. The $(n-1)$-th derivative is not continuous, but suffers a discontinuous jump there. Integrating again shows that $G^{(n−2)}$ is continuous".



My problem here is the following:
how can I say that the result of integrating equation $(1)$ gets me "some continuous functions" if I don't know anything about the continuity of $G(x, x_0)$ and its derivatives?



Second approach



I found it for a $2^{nd}$ order ODE but I assume it can be generalized. It proceeds by saying "Suppose first that $G(x,x_0)$ was discontinuous at $x=x_0$, with the discontinuity modelled by a step function. Then $G' propto delta(x-x_0)$ and consequently $G''propto delta(x−x_0)$. However, the form of the ODE shows that $mathcal{L}G$ involves no generalized functions beyond $delta(x-x_0)$, thus $G(x, ξ)$ must be continuous." He THEN proceeds in integrating the ODE in a small neighborhood of $x_0$, and he asserts that "Also, since $G$ is continuous, $G'$ must be bounded so the term $G'$ also cannot contribute as the integration region shrinks to zero size". He then concludes that it's only the second derivative that contributes to the magnitude 1 jump.



My problems here:



1) we started by modeling the discontinuity of $G$ by a step function, noticed that there was a contradiction and consequently stated that "so the function must be continuous". But isn't it a bit restrictive, to just assume a jump discontinuity?



2) Can this be generalized for a generic $n$-order equation?



My questions



To resume, I would like to know:



1) What is the correct way to proceed to demonstrate that $G$ and its derivatives up to the $n-2$-th are continuous in $x_0$ and its $(n-1)$-th derivative has a jump discontinuity



2) As a side question: if the $(n-1)$-th derivative of the function is discontinuous in $x_0$ how can $G^{n}$ exist in $x_0$?



I'm a physicist not a mathematician, so please be kind =)










share|cite|improve this question











$endgroup$

















    1












    $begingroup$


    I am studying Green functions and I want to find the Green function $G(x,x_0)$ of the linear operator
    begin{equation}
    mathcal{L}=left(frac{d^n}{dx^n}+a_1(x)frac{d^{n-1}}{dx^{n-1}}+dots+a_n(x) right)
    end{equation}

    This means solving the following problem
    begin{equation}
    frac{d^n G}{dx^n}+a_1(x)frac{d^{n-1}G}{dx^{n-1}}+dots+a_n(x)G=delta(x-x_0)
    tag{1}
    end{equation}



    I understand that I have to impose continuity in $x_0$ for all the derivatives, up to $n-2$, and impose the jump discontinuity in $x_0$ for the $(n-1)$-th derivative. I don't exactly get why this is the case. I have searched in some lecture notes and found the following clues:



    First approach



    They proceed by integrating equation $(1)$, and they say "we obtain$G^{(n-1)}=H(x-x_0)$ + some continuous functions. The $(n-1)$-th derivative is not continuous, but suffers a discontinuous jump there. Integrating again shows that $G^{(n−2)}$ is continuous".



    My problem here is the following:
    how can I say that the result of integrating equation $(1)$ gets me "some continuous functions" if I don't know anything about the continuity of $G(x, x_0)$ and its derivatives?



    Second approach



    I found it for a $2^{nd}$ order ODE but I assume it can be generalized. It proceeds by saying "Suppose first that $G(x,x_0)$ was discontinuous at $x=x_0$, with the discontinuity modelled by a step function. Then $G' propto delta(x-x_0)$ and consequently $G''propto delta(x−x_0)$. However, the form of the ODE shows that $mathcal{L}G$ involves no generalized functions beyond $delta(x-x_0)$, thus $G(x, ξ)$ must be continuous." He THEN proceeds in integrating the ODE in a small neighborhood of $x_0$, and he asserts that "Also, since $G$ is continuous, $G'$ must be bounded so the term $G'$ also cannot contribute as the integration region shrinks to zero size". He then concludes that it's only the second derivative that contributes to the magnitude 1 jump.



    My problems here:



    1) we started by modeling the discontinuity of $G$ by a step function, noticed that there was a contradiction and consequently stated that "so the function must be continuous". But isn't it a bit restrictive, to just assume a jump discontinuity?



    2) Can this be generalized for a generic $n$-order equation?



    My questions



    To resume, I would like to know:



    1) What is the correct way to proceed to demonstrate that $G$ and its derivatives up to the $n-2$-th are continuous in $x_0$ and its $(n-1)$-th derivative has a jump discontinuity



    2) As a side question: if the $(n-1)$-th derivative of the function is discontinuous in $x_0$ how can $G^{n}$ exist in $x_0$?



    I'm a physicist not a mathematician, so please be kind =)










    share|cite|improve this question











    $endgroup$















      1












      1








      1





      $begingroup$


      I am studying Green functions and I want to find the Green function $G(x,x_0)$ of the linear operator
      begin{equation}
      mathcal{L}=left(frac{d^n}{dx^n}+a_1(x)frac{d^{n-1}}{dx^{n-1}}+dots+a_n(x) right)
      end{equation}

      This means solving the following problem
      begin{equation}
      frac{d^n G}{dx^n}+a_1(x)frac{d^{n-1}G}{dx^{n-1}}+dots+a_n(x)G=delta(x-x_0)
      tag{1}
      end{equation}



      I understand that I have to impose continuity in $x_0$ for all the derivatives, up to $n-2$, and impose the jump discontinuity in $x_0$ for the $(n-1)$-th derivative. I don't exactly get why this is the case. I have searched in some lecture notes and found the following clues:



      First approach



      They proceed by integrating equation $(1)$, and they say "we obtain$G^{(n-1)}=H(x-x_0)$ + some continuous functions. The $(n-1)$-th derivative is not continuous, but suffers a discontinuous jump there. Integrating again shows that $G^{(n−2)}$ is continuous".



      My problem here is the following:
      how can I say that the result of integrating equation $(1)$ gets me "some continuous functions" if I don't know anything about the continuity of $G(x, x_0)$ and its derivatives?



      Second approach



      I found it for a $2^{nd}$ order ODE but I assume it can be generalized. It proceeds by saying "Suppose first that $G(x,x_0)$ was discontinuous at $x=x_0$, with the discontinuity modelled by a step function. Then $G' propto delta(x-x_0)$ and consequently $G''propto delta(x−x_0)$. However, the form of the ODE shows that $mathcal{L}G$ involves no generalized functions beyond $delta(x-x_0)$, thus $G(x, ξ)$ must be continuous." He THEN proceeds in integrating the ODE in a small neighborhood of $x_0$, and he asserts that "Also, since $G$ is continuous, $G'$ must be bounded so the term $G'$ also cannot contribute as the integration region shrinks to zero size". He then concludes that it's only the second derivative that contributes to the magnitude 1 jump.



      My problems here:



      1) we started by modeling the discontinuity of $G$ by a step function, noticed that there was a contradiction and consequently stated that "so the function must be continuous". But isn't it a bit restrictive, to just assume a jump discontinuity?



      2) Can this be generalized for a generic $n$-order equation?



      My questions



      To resume, I would like to know:



      1) What is the correct way to proceed to demonstrate that $G$ and its derivatives up to the $n-2$-th are continuous in $x_0$ and its $(n-1)$-th derivative has a jump discontinuity



      2) As a side question: if the $(n-1)$-th derivative of the function is discontinuous in $x_0$ how can $G^{n}$ exist in $x_0$?



      I'm a physicist not a mathematician, so please be kind =)










      share|cite|improve this question











      $endgroup$




      I am studying Green functions and I want to find the Green function $G(x,x_0)$ of the linear operator
      begin{equation}
      mathcal{L}=left(frac{d^n}{dx^n}+a_1(x)frac{d^{n-1}}{dx^{n-1}}+dots+a_n(x) right)
      end{equation}

      This means solving the following problem
      begin{equation}
      frac{d^n G}{dx^n}+a_1(x)frac{d^{n-1}G}{dx^{n-1}}+dots+a_n(x)G=delta(x-x_0)
      tag{1}
      end{equation}



      I understand that I have to impose continuity in $x_0$ for all the derivatives, up to $n-2$, and impose the jump discontinuity in $x_0$ for the $(n-1)$-th derivative. I don't exactly get why this is the case. I have searched in some lecture notes and found the following clues:



      First approach



      They proceed by integrating equation $(1)$, and they say "we obtain$G^{(n-1)}=H(x-x_0)$ + some continuous functions. The $(n-1)$-th derivative is not continuous, but suffers a discontinuous jump there. Integrating again shows that $G^{(n−2)}$ is continuous".



      My problem here is the following:
      how can I say that the result of integrating equation $(1)$ gets me "some continuous functions" if I don't know anything about the continuity of $G(x, x_0)$ and its derivatives?



      Second approach



      I found it for a $2^{nd}$ order ODE but I assume it can be generalized. It proceeds by saying "Suppose first that $G(x,x_0)$ was discontinuous at $x=x_0$, with the discontinuity modelled by a step function. Then $G' propto delta(x-x_0)$ and consequently $G''propto delta(x−x_0)$. However, the form of the ODE shows that $mathcal{L}G$ involves no generalized functions beyond $delta(x-x_0)$, thus $G(x, ξ)$ must be continuous." He THEN proceeds in integrating the ODE in a small neighborhood of $x_0$, and he asserts that "Also, since $G$ is continuous, $G'$ must be bounded so the term $G'$ also cannot contribute as the integration region shrinks to zero size". He then concludes that it's only the second derivative that contributes to the magnitude 1 jump.



      My problems here:



      1) we started by modeling the discontinuity of $G$ by a step function, noticed that there was a contradiction and consequently stated that "so the function must be continuous". But isn't it a bit restrictive, to just assume a jump discontinuity?



      2) Can this be generalized for a generic $n$-order equation?



      My questions



      To resume, I would like to know:



      1) What is the correct way to proceed to demonstrate that $G$ and its derivatives up to the $n-2$-th are continuous in $x_0$ and its $(n-1)$-th derivative has a jump discontinuity



      2) As a side question: if the $(n-1)$-th derivative of the function is discontinuous in $x_0$ how can $G^{n}$ exist in $x_0$?



      I'm a physicist not a mathematician, so please be kind =)







      ordinary-differential-equations greens-function






      share|cite|improve this question















      share|cite|improve this question













      share|cite|improve this question




      share|cite|improve this question








      edited Jan 13 at 2:12







      Luthien

















      asked Jan 13 at 2:00









      LuthienLuthien

      1448




      1448






















          1 Answer
          1






          active

          oldest

          votes


















          1












          $begingroup$

          1) One can start other way round. Namely to present a Green's function with required properties. Let for simplicity sake the equation's coefficients be constant. Let $v$ be the solution of the Cauchy problem ${cal L}v=0$, $v(x_0)=v'(x_0)=ldots=v^{(n-2)}(x_0)=0$, $v^{(n-1)}(x_0)=1$ and put $G(x,x_0)=v(x)$ for $x>x_0$ and $G(x,x_0)=0$ for $x<x_0$. It's straightforward to check that $G(x,x_0)$ is a Green's function for ${cal L}$.



          Btw in the PDE theory such functions are called also fundamental solutions and the term Green's function is usually reserved for fundamental solutions with some homogeneous boundary conditions (e.g. zero Dirichlet condition).



          By construction function $G(x,x_0)$ has continuous derivatives up to order $n-2$ and the derivative of order $n-1$ has a jump discontinuity in $x_0$. Now if $tilde G(x,x_0)$ is another Green's function then their difference $G-tilde G$ satisfy homogeneous equation ${cal L}u=0$ and therefore is smooth. In particular it has a continuous derivative of order $n-1$. So derivative $tilde G^{(n-1)}(x,x_0)$ also has a jump discontinuity...



          2) It doesn't exist, no problem here. The Dirac function $delta(x-x_0)$ also "does not exist" in $x_0$. The equation with the delta function in the RHS is understood in the sense of distributions.



          ADDED:



          I don't see the need to follow their argument to the point.



          To see that the required parts are "some continuous functions" one can start other way round without assuming any smoothness of $G$ but keeping in mind the theorems about smoothness of ODE solutions (with regular enough coefficients). Fix $x_0$ and define function $v$ as $v(x)=(x-x_0)^{n-1}/(n-1)!$ for $xge x_0$ and $v(x)=0$ for $x< x_0$. Since $partial_x^n v(x)=delta(x-x_0)$ for function $u(x)=G(x,x_0)-v(x)$ it follows that
          $$
          cal Lu=-(L-partial_x^n)v=f(x),
          $$

          where the rhs $f$ is bounded in some segment containing $x_0$, say $I=[x_0-1,x_0+1]$. In the (linear) ODE theory there are theorems stating that smoothness of the solution's senior derivative $partial_x^nu$ has the same smoothness as the rhs. In particular if $fin L_infty(I)$ then $partial_x^nuin L_infty(I)$. And this means that $partial_x^{n-1}u$ is continuous (even Lipschitz continuous). Returning to the Green function one has
          $$
          partial_x^{n-1}G(x,x_0)=partial_x^{n-1}u(x)+partial_x^{n-1}v(x)=
          partial_x^{n-1}u(x)+H(x-x_0)
          $$

          which is a sum of a continuous function and the step function.






          share|cite|improve this answer











          $endgroup$













          • $begingroup$
            Thank you for your answer! Can I also ask you what do you think of the two approaches I mentioned in my question? Are they "wrong" in some part?
            $endgroup$
            – Luthien
            Jan 15 at 15:36






          • 1




            $begingroup$
            @Luthien I'e edited the answer rewriting the arguments from your points 1) and 2) in somewhat other form. But essentially it's all the same.
            $endgroup$
            – Andrew
            Jan 16 at 12:59












          • $begingroup$
            Thank you, great!
            $endgroup$
            – Luthien
            Jan 16 at 13:40











          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%2f3071624%2fgreen-function-for-solving-ode-continuity-and-discontinuity-of-derivatives%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$

          1) One can start other way round. Namely to present a Green's function with required properties. Let for simplicity sake the equation's coefficients be constant. Let $v$ be the solution of the Cauchy problem ${cal L}v=0$, $v(x_0)=v'(x_0)=ldots=v^{(n-2)}(x_0)=0$, $v^{(n-1)}(x_0)=1$ and put $G(x,x_0)=v(x)$ for $x>x_0$ and $G(x,x_0)=0$ for $x<x_0$. It's straightforward to check that $G(x,x_0)$ is a Green's function for ${cal L}$.



          Btw in the PDE theory such functions are called also fundamental solutions and the term Green's function is usually reserved for fundamental solutions with some homogeneous boundary conditions (e.g. zero Dirichlet condition).



          By construction function $G(x,x_0)$ has continuous derivatives up to order $n-2$ and the derivative of order $n-1$ has a jump discontinuity in $x_0$. Now if $tilde G(x,x_0)$ is another Green's function then their difference $G-tilde G$ satisfy homogeneous equation ${cal L}u=0$ and therefore is smooth. In particular it has a continuous derivative of order $n-1$. So derivative $tilde G^{(n-1)}(x,x_0)$ also has a jump discontinuity...



          2) It doesn't exist, no problem here. The Dirac function $delta(x-x_0)$ also "does not exist" in $x_0$. The equation with the delta function in the RHS is understood in the sense of distributions.



          ADDED:



          I don't see the need to follow their argument to the point.



          To see that the required parts are "some continuous functions" one can start other way round without assuming any smoothness of $G$ but keeping in mind the theorems about smoothness of ODE solutions (with regular enough coefficients). Fix $x_0$ and define function $v$ as $v(x)=(x-x_0)^{n-1}/(n-1)!$ for $xge x_0$ and $v(x)=0$ for $x< x_0$. Since $partial_x^n v(x)=delta(x-x_0)$ for function $u(x)=G(x,x_0)-v(x)$ it follows that
          $$
          cal Lu=-(L-partial_x^n)v=f(x),
          $$

          where the rhs $f$ is bounded in some segment containing $x_0$, say $I=[x_0-1,x_0+1]$. In the (linear) ODE theory there are theorems stating that smoothness of the solution's senior derivative $partial_x^nu$ has the same smoothness as the rhs. In particular if $fin L_infty(I)$ then $partial_x^nuin L_infty(I)$. And this means that $partial_x^{n-1}u$ is continuous (even Lipschitz continuous). Returning to the Green function one has
          $$
          partial_x^{n-1}G(x,x_0)=partial_x^{n-1}u(x)+partial_x^{n-1}v(x)=
          partial_x^{n-1}u(x)+H(x-x_0)
          $$

          which is a sum of a continuous function and the step function.






          share|cite|improve this answer











          $endgroup$













          • $begingroup$
            Thank you for your answer! Can I also ask you what do you think of the two approaches I mentioned in my question? Are they "wrong" in some part?
            $endgroup$
            – Luthien
            Jan 15 at 15:36






          • 1




            $begingroup$
            @Luthien I'e edited the answer rewriting the arguments from your points 1) and 2) in somewhat other form. But essentially it's all the same.
            $endgroup$
            – Andrew
            Jan 16 at 12:59












          • $begingroup$
            Thank you, great!
            $endgroup$
            – Luthien
            Jan 16 at 13:40
















          1












          $begingroup$

          1) One can start other way round. Namely to present a Green's function with required properties. Let for simplicity sake the equation's coefficients be constant. Let $v$ be the solution of the Cauchy problem ${cal L}v=0$, $v(x_0)=v'(x_0)=ldots=v^{(n-2)}(x_0)=0$, $v^{(n-1)}(x_0)=1$ and put $G(x,x_0)=v(x)$ for $x>x_0$ and $G(x,x_0)=0$ for $x<x_0$. It's straightforward to check that $G(x,x_0)$ is a Green's function for ${cal L}$.



          Btw in the PDE theory such functions are called also fundamental solutions and the term Green's function is usually reserved for fundamental solutions with some homogeneous boundary conditions (e.g. zero Dirichlet condition).



          By construction function $G(x,x_0)$ has continuous derivatives up to order $n-2$ and the derivative of order $n-1$ has a jump discontinuity in $x_0$. Now if $tilde G(x,x_0)$ is another Green's function then their difference $G-tilde G$ satisfy homogeneous equation ${cal L}u=0$ and therefore is smooth. In particular it has a continuous derivative of order $n-1$. So derivative $tilde G^{(n-1)}(x,x_0)$ also has a jump discontinuity...



          2) It doesn't exist, no problem here. The Dirac function $delta(x-x_0)$ also "does not exist" in $x_0$. The equation with the delta function in the RHS is understood in the sense of distributions.



          ADDED:



          I don't see the need to follow their argument to the point.



          To see that the required parts are "some continuous functions" one can start other way round without assuming any smoothness of $G$ but keeping in mind the theorems about smoothness of ODE solutions (with regular enough coefficients). Fix $x_0$ and define function $v$ as $v(x)=(x-x_0)^{n-1}/(n-1)!$ for $xge x_0$ and $v(x)=0$ for $x< x_0$. Since $partial_x^n v(x)=delta(x-x_0)$ for function $u(x)=G(x,x_0)-v(x)$ it follows that
          $$
          cal Lu=-(L-partial_x^n)v=f(x),
          $$

          where the rhs $f$ is bounded in some segment containing $x_0$, say $I=[x_0-1,x_0+1]$. In the (linear) ODE theory there are theorems stating that smoothness of the solution's senior derivative $partial_x^nu$ has the same smoothness as the rhs. In particular if $fin L_infty(I)$ then $partial_x^nuin L_infty(I)$. And this means that $partial_x^{n-1}u$ is continuous (even Lipschitz continuous). Returning to the Green function one has
          $$
          partial_x^{n-1}G(x,x_0)=partial_x^{n-1}u(x)+partial_x^{n-1}v(x)=
          partial_x^{n-1}u(x)+H(x-x_0)
          $$

          which is a sum of a continuous function and the step function.






          share|cite|improve this answer











          $endgroup$













          • $begingroup$
            Thank you for your answer! Can I also ask you what do you think of the two approaches I mentioned in my question? Are they "wrong" in some part?
            $endgroup$
            – Luthien
            Jan 15 at 15:36






          • 1




            $begingroup$
            @Luthien I'e edited the answer rewriting the arguments from your points 1) and 2) in somewhat other form. But essentially it's all the same.
            $endgroup$
            – Andrew
            Jan 16 at 12:59












          • $begingroup$
            Thank you, great!
            $endgroup$
            – Luthien
            Jan 16 at 13:40














          1












          1








          1





          $begingroup$

          1) One can start other way round. Namely to present a Green's function with required properties. Let for simplicity sake the equation's coefficients be constant. Let $v$ be the solution of the Cauchy problem ${cal L}v=0$, $v(x_0)=v'(x_0)=ldots=v^{(n-2)}(x_0)=0$, $v^{(n-1)}(x_0)=1$ and put $G(x,x_0)=v(x)$ for $x>x_0$ and $G(x,x_0)=0$ for $x<x_0$. It's straightforward to check that $G(x,x_0)$ is a Green's function for ${cal L}$.



          Btw in the PDE theory such functions are called also fundamental solutions and the term Green's function is usually reserved for fundamental solutions with some homogeneous boundary conditions (e.g. zero Dirichlet condition).



          By construction function $G(x,x_0)$ has continuous derivatives up to order $n-2$ and the derivative of order $n-1$ has a jump discontinuity in $x_0$. Now if $tilde G(x,x_0)$ is another Green's function then their difference $G-tilde G$ satisfy homogeneous equation ${cal L}u=0$ and therefore is smooth. In particular it has a continuous derivative of order $n-1$. So derivative $tilde G^{(n-1)}(x,x_0)$ also has a jump discontinuity...



          2) It doesn't exist, no problem here. The Dirac function $delta(x-x_0)$ also "does not exist" in $x_0$. The equation with the delta function in the RHS is understood in the sense of distributions.



          ADDED:



          I don't see the need to follow their argument to the point.



          To see that the required parts are "some continuous functions" one can start other way round without assuming any smoothness of $G$ but keeping in mind the theorems about smoothness of ODE solutions (with regular enough coefficients). Fix $x_0$ and define function $v$ as $v(x)=(x-x_0)^{n-1}/(n-1)!$ for $xge x_0$ and $v(x)=0$ for $x< x_0$. Since $partial_x^n v(x)=delta(x-x_0)$ for function $u(x)=G(x,x_0)-v(x)$ it follows that
          $$
          cal Lu=-(L-partial_x^n)v=f(x),
          $$

          where the rhs $f$ is bounded in some segment containing $x_0$, say $I=[x_0-1,x_0+1]$. In the (linear) ODE theory there are theorems stating that smoothness of the solution's senior derivative $partial_x^nu$ has the same smoothness as the rhs. In particular if $fin L_infty(I)$ then $partial_x^nuin L_infty(I)$. And this means that $partial_x^{n-1}u$ is continuous (even Lipschitz continuous). Returning to the Green function one has
          $$
          partial_x^{n-1}G(x,x_0)=partial_x^{n-1}u(x)+partial_x^{n-1}v(x)=
          partial_x^{n-1}u(x)+H(x-x_0)
          $$

          which is a sum of a continuous function and the step function.






          share|cite|improve this answer











          $endgroup$



          1) One can start other way round. Namely to present a Green's function with required properties. Let for simplicity sake the equation's coefficients be constant. Let $v$ be the solution of the Cauchy problem ${cal L}v=0$, $v(x_0)=v'(x_0)=ldots=v^{(n-2)}(x_0)=0$, $v^{(n-1)}(x_0)=1$ and put $G(x,x_0)=v(x)$ for $x>x_0$ and $G(x,x_0)=0$ for $x<x_0$. It's straightforward to check that $G(x,x_0)$ is a Green's function for ${cal L}$.



          Btw in the PDE theory such functions are called also fundamental solutions and the term Green's function is usually reserved for fundamental solutions with some homogeneous boundary conditions (e.g. zero Dirichlet condition).



          By construction function $G(x,x_0)$ has continuous derivatives up to order $n-2$ and the derivative of order $n-1$ has a jump discontinuity in $x_0$. Now if $tilde G(x,x_0)$ is another Green's function then their difference $G-tilde G$ satisfy homogeneous equation ${cal L}u=0$ and therefore is smooth. In particular it has a continuous derivative of order $n-1$. So derivative $tilde G^{(n-1)}(x,x_0)$ also has a jump discontinuity...



          2) It doesn't exist, no problem here. The Dirac function $delta(x-x_0)$ also "does not exist" in $x_0$. The equation with the delta function in the RHS is understood in the sense of distributions.



          ADDED:



          I don't see the need to follow their argument to the point.



          To see that the required parts are "some continuous functions" one can start other way round without assuming any smoothness of $G$ but keeping in mind the theorems about smoothness of ODE solutions (with regular enough coefficients). Fix $x_0$ and define function $v$ as $v(x)=(x-x_0)^{n-1}/(n-1)!$ for $xge x_0$ and $v(x)=0$ for $x< x_0$. Since $partial_x^n v(x)=delta(x-x_0)$ for function $u(x)=G(x,x_0)-v(x)$ it follows that
          $$
          cal Lu=-(L-partial_x^n)v=f(x),
          $$

          where the rhs $f$ is bounded in some segment containing $x_0$, say $I=[x_0-1,x_0+1]$. In the (linear) ODE theory there are theorems stating that smoothness of the solution's senior derivative $partial_x^nu$ has the same smoothness as the rhs. In particular if $fin L_infty(I)$ then $partial_x^nuin L_infty(I)$. And this means that $partial_x^{n-1}u$ is continuous (even Lipschitz continuous). Returning to the Green function one has
          $$
          partial_x^{n-1}G(x,x_0)=partial_x^{n-1}u(x)+partial_x^{n-1}v(x)=
          partial_x^{n-1}u(x)+H(x-x_0)
          $$

          which is a sum of a continuous function and the step function.







          share|cite|improve this answer














          share|cite|improve this answer



          share|cite|improve this answer








          edited Jan 16 at 12:57

























          answered Jan 13 at 3:12









          AndrewAndrew

          9,55211946




          9,55211946












          • $begingroup$
            Thank you for your answer! Can I also ask you what do you think of the two approaches I mentioned in my question? Are they "wrong" in some part?
            $endgroup$
            – Luthien
            Jan 15 at 15:36






          • 1




            $begingroup$
            @Luthien I'e edited the answer rewriting the arguments from your points 1) and 2) in somewhat other form. But essentially it's all the same.
            $endgroup$
            – Andrew
            Jan 16 at 12:59












          • $begingroup$
            Thank you, great!
            $endgroup$
            – Luthien
            Jan 16 at 13:40


















          • $begingroup$
            Thank you for your answer! Can I also ask you what do you think of the two approaches I mentioned in my question? Are they "wrong" in some part?
            $endgroup$
            – Luthien
            Jan 15 at 15:36






          • 1




            $begingroup$
            @Luthien I'e edited the answer rewriting the arguments from your points 1) and 2) in somewhat other form. But essentially it's all the same.
            $endgroup$
            – Andrew
            Jan 16 at 12:59












          • $begingroup$
            Thank you, great!
            $endgroup$
            – Luthien
            Jan 16 at 13:40
















          $begingroup$
          Thank you for your answer! Can I also ask you what do you think of the two approaches I mentioned in my question? Are they "wrong" in some part?
          $endgroup$
          – Luthien
          Jan 15 at 15:36




          $begingroup$
          Thank you for your answer! Can I also ask you what do you think of the two approaches I mentioned in my question? Are they "wrong" in some part?
          $endgroup$
          – Luthien
          Jan 15 at 15:36




          1




          1




          $begingroup$
          @Luthien I'e edited the answer rewriting the arguments from your points 1) and 2) in somewhat other form. But essentially it's all the same.
          $endgroup$
          – Andrew
          Jan 16 at 12:59






          $begingroup$
          @Luthien I'e edited the answer rewriting the arguments from your points 1) and 2) in somewhat other form. But essentially it's all the same.
          $endgroup$
          – Andrew
          Jan 16 at 12:59














          $begingroup$
          Thank you, great!
          $endgroup$
          – Luthien
          Jan 16 at 13:40




          $begingroup$
          Thank you, great!
          $endgroup$
          – Luthien
          Jan 16 at 13:40


















          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%2f3071624%2fgreen-function-for-solving-ode-continuity-and-discontinuity-of-derivatives%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

          Can a sorcerer learn a 5th-level spell early by creating spell slots using the Font of Magic feature?

          Does disintegrating a polymorphed enemy still kill it after the 2018 errata?

          A Topological Invariant for $pi_3(U(n))$