Existence and Uniqueness Theorems for ODEs

  • Thread starter Thread starter twoflower
  • Start date Start date
  • Tags Tags
    Odes
twoflower
Messages
363
Reaction score
0
Hi all,
I'd be happy if someone could clarify these two things to me:

1. While solving linear first-order ODE, I first solve homogenous equation (with the right side equal to 0) and eventually I get to the point (just an example):

<br /> \log |y| = \log C(e^{x} - 1)<br />

Now, is it ok to compute C for given non-homogenous equation and then write there two solutions

<br /> y = C(e^{x} - 1)<br />

<br /> y = -C(e^{x} - 1)\mbox{ ?}<br />

Because, you know, both satisfy

<br /> \log |y| = \log C(e^{x} - 1)<br />

Anyway, sometimes it gave me correct results (both C(e^{x} - 1) and -C(e^{x} - 1) were solutions) while other times only plus-signed solution was ok.
2. Solving Bernoulli's ODE, let's say

<br /> y&#039; = 2y + 2x\sqrt{y}<br />

and substituting

<br /> z = \sqrt{y}<br />

<br /> y = z^2<br />

<br /> y&#039; = 2zz&#039;<br />

We get

<br /> 2zz&#039; = 2z^2 + 2xz<br />

and first thing I do is dividing with 2z and so getting the condition z \neq 0
After some computing, I get the result

<br /> z = Ce^{x} - x - 1<br />

and thus

<br /> y = (Ce^{x} - x - 1)^2<br />

a) first question
I'd suppose this expression must not get zeroed because at the beginning we divided the equation with 2z. Anyway, even if C=1 and x=0 (and thus z = 0 = y), the original equation holds true. How it comes?b) second question
Having the substitution above in mind, I know that z itself must be \geq 0. Anyway, if we take C \in (0, 1), then the equation

<br /> Ce^{x} -x - 1 = 0<br />

has two roots, x_0 and x_1. Then we know that
for x \in (x_0, x_1)

<br /> Ce^{x} -x - 1 &lt; 0<br />

Ok, still remembering that this is z and it must be \geq 0, we shouldn't accept this interval for x. However, what's the problem with this solution on this interval? When you put it in the original equation, it's ok I think.

I can't understand why the interval (x_0,x_1) is excluded from the solution (according to our professor).Thank you for any suggestions.
 
Last edited:
Physics news on Phys.org
twoflower said:
b) second question
Having the substitution above in mind, I know that z itself must be \geq 0. Anyway, if we take C \in (0, 1), then the equation

<br /> Ce^{x} -x - 1 = 0<br />

has two roots, x_0 and x_1. Then we know that
for x \in (x_0, x_1)

<br /> Ce^{x} -x - 1 &lt; 0<br />

Ok, still remembering that this is z and it must be \geq 0, we shouldn't accept this interval for x. However, what's the problem with this solution on this interval? When you put it in the original equation, it's ok I think.

I can't understand why the interval (x_0,x_1) is excluded from the solution (according to our professor).


Thank you for any suggestions.

Well, I would've lost money on that one Twoflower. Apparently when you assign:

z=y^{1/2}

You implicitly define z\geq 0.

Note that the solution:

y(t)=(ce^x-(1+x))^2

does NOT satisfiy the ODE in the interval where z is less than zero. For example, take the plot below for c=0.5 and:

y(t)=(ce^x-(1+x))^2

For (approx):

x\in (-0.77,1.68)

the curve does not satisfy the ODE. New for me:eek: . Thanks!
 

Attachments

  • Bernoulli.JPG
    Bernoulli.JPG
    5.1 KB · Views: 433
I wish to clear up something with this problem from the perspective of existence and uniqueness:

Consider:

y^{&#039;}=f(x,y)=2y+2x\sqrt{y},\quad y(0)=0.25

Now,since f(x,y) and the partial of f with respect to y exists in a bounded region about the point (0,0.25), we can expect a unique solution passing through this point.

The solution above was determined to be:

y(x)=(ce^x-1-x)^2\;\text{with}\;ce^x-1-x\geq 0

Solving for c, I obtain:

c=0.5\;\text{and}\;c=1.5

Now, initially this would suggest two solutions satisfy the equation. However, the value of c=0.5 would not satisify the restiction ce^x-1-x\geq 0[/tex] and so must be discarded. I&#039;m left then with the &quot;unique&quot; solution:<br /> <br /> y(x)=(1.5e^x-1-x)^2
 
Thank you Saltydog, now I'm clear about that, it helped me a lot.

saltydog said:
Now,since f(x,y) and the partial of f with respect to y exists in a bounded region about the point (0,0.25), we can expect a unique solution passing through this point.
What does this follow from? Is there some theorem saying this?

Could you clarify those two other questions to me too? Thank you very much.
 
Last edited:
twoflower said:
What does this follow from? Is there some theorem saying this?

This is based on Existence and Uniqueness Theorems for ordinary differential equations:

Consider:

y^{&#039;}=f(x,y)

Let T denote the rectangular region defined by:

|x-x_0|\leq a\quad |y-y_0|\leq b

with the point (x_0,y_0) at its center. Let f and \frac{\partial f}{\partial y} be continuous at each point in T, then there exists an interval|x-x_0|\leq h and a function y(x) such that:

1. y(x) is a solution to the ODE on the interval |x-x_0|\leq h.

2. y(x) is unique in the sense that if another equation h(x) satisfies the ODE, then y(x)=h(x).

Note that the interval |x-x_0|\leq h can be very small and nothing is said about uniqueness and existence if f and the partial do not exist at the point x0.

Thus for a particular IVP, I know if the above is satisfied, I can be guaranteed a unique solution. However, I'm not sure about the case of y=0. Note that for the IVP:

y^{&#039;}=2y+2x\sqrt{y};\quad y(0)=0

The partial is indeterminant at the initial condition so uniqueness and existence cannot be guaranteed.

Could you clarify those two other questions to me too? Thank you very much.

Not sure I can do better than the above. Perhaps someone in here more knowledgeable than I could elaborate on the matter.
 
Last edited:
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top