# ODE Reduction in Order

**URL:** <https://boards.straightdope.com/t/ode-reduction-in-order/668402>\
**Category:** Factual Questions\
**Created:** [September 9, 2013, 11:46am UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402 "2013-09-09T11:46:04Z")\
**Posts on this page:** 13\
**Page:** 1

<div class="post-metadata">

**Author:** ![dgrdfd](https://avatars.discourse-cdn.com/v4/letter/d/eb8c5e/32.png) [@dgrdfd](https://boards.straightdope.com/u/dgrdfd)\
**Post date:** [September 9, 2013, 11:46am UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/1 "2013-09-09T11:46:04Z")

</div>

First, yes this is homework, but no I am not looking for someone to just give me the answers. Getting the answers is easy enough with Matlab/Maple/Wolfram Alpha, etc. I just can’t figure out how to solve these two problems by reducing them to a first order ODE. As I understand it,

1. The general solution y=y1+y2 where y1 and y2 are each solutions, given y1 is a solution then y2=u(x) \* y1

2. Find the first and second derivative of y2 and then plug into the original ODE

3. Now set w=u’ and the second order ODE has been reduced to a first order.

4. Solve for w. Then integrate to find u.

5. Now, take the original (given, or obvious by inspection solution) y1 and multiple by u to solve the ODE and find y

That’s a rough outline of my plan. However, I can’t find an “obvious” y1 solution and none is given. Here are the ODEs…

A: y_y’’=3_(y’)^2  
B: y’’=1+(y’)^2

I’ve tried solutions of the form y=x^m, y=e^x, etc. Any advice would be appreciated.

---

<div class="post-metadata">

**Author:** ![Thudlow\_Boink](https://sea3.discourse-cdn.com/straightdope/user_avatar/boards.straightdope.com/thudlow_boink/32/320_2.png) [@Thudlow\_Boink](https://boards.straightdope.com/u/Thudlow_Boink)\
**Post date:** [September 9, 2013, 2:00pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/2 "2013-09-09T14:00:11Z")

</div>

For B, if I understand your notation, it looks like there’s no y, just y’ and y’’. So couldn’t you reduce it to order 1 by simply setting w = y’ ?

---

<div class="post-metadata">

**Author:** ![MikeS](https://avatars.discourse-cdn.com/v4/letter/m/919ad9/32.png) [@MikeS](https://boards.straightdope.com/u/MikeS)\
**Post date:** [September 9, 2013, 2:08pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/3 "2013-09-09T14:08:38Z")

</div>

For the first one: I thought I had an idea how to do this, but I didn’t. Let me get back to you on this one.

I don’t think you need your whole procedure for the second one. There’s a way to reduce that one in order without “guessing” a solution. ETA: which **Thudlow Boink** described without being as cagey as I was being.

By the way:

> [@dgrdfd](#):
>
> 1. The general solution y=y1+y2 where y1 and y2 are each solutions…

The general solution of a _linear_ ODE is the arbitrary linear combination of two solutions y[sub]1[/sub] and y[sub]2[/sub]. The ODEs you have here aren’t linear.

Hope this helps! Feel free to ask follow-up questions.

---

<div class="post-metadata">

**Author:** ![MikeS](https://avatars.discourse-cdn.com/v4/letter/m/919ad9/32.png) [@MikeS](https://boards.straightdope.com/u/MikeS)\
**Post date:** [September 9, 2013, 2:17pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/4 "2013-09-09T14:17:49Z")

</div>

For the first one: try a power law. This is the form you should “guess” when the equation is _homogeneous_ in both the independent and the dependent variable. (Homogeneous means that if you rescale the variable by a factor lambda, you get the same equation back.)

---

<div class="post-metadata">

**Author:** ![MikeS](https://avatars.discourse-cdn.com/v4/letter/m/919ad9/32.png) [@MikeS](https://boards.straightdope.com/u/MikeS)\
**Post date:** [September 9, 2013, 2:58pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/5 "2013-09-09T14:58:44Z")

</div>

> [@MikeS](#):
>
> The general solution of a _linear_ ODE is the arbitrary linear combination of two solutions y[sub]1[/sub] and y[sub]2[/sub]. The ODEs you have here aren’t linear.

… and upon further reflection, I realized that the method described in the OP doesn’t work unless the ODE is linear. You’ll need a different technique to reduce the order of the ODE in question A. The fact that the ODE is homogeneous (as I mentioned in my previous post) is useful.

---

<div class="post-metadata">

**Author:** ![Chronos](https://sea3.discourse-cdn.com/straightdope/user_avatar/boards.straightdope.com/chronos/32/134_2.png) [@Chronos](https://boards.straightdope.com/u/Chronos)\
**Post date:** [September 9, 2013, 5:14pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/6 "2013-09-09T17:14:47Z")

</div>

While we’re at it, if you’re setting up a problem for an eventual numerical solution, you can _always_ reduce the order of a system of differential equations to 1, by increasing the number of equations in the system. Just define a new variable for y’, and another one for y’’, and so on as necessary. Though this doesn’t appear to be what the OP is looking for in this case.

---

<div class="post-metadata">

**Author:** ![dgrdfd](https://avatars.discourse-cdn.com/v4/letter/d/eb8c5e/32.png) [@dgrdfd](https://boards.straightdope.com/u/dgrdfd)\
**Post date:** [September 10, 2013, 11:44am UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/7 "2013-09-10T11:44:32Z")

</div>

First, sorry it has taken me so long to respond. My phone is giving me issues posting from mobile so I have to actually be at a computer so I couldn’t respond as quickly as i would have liked.

Second, thanks very much for the help guys. I think my fundamental misconception was that I was only thinking of the linear case and assumed I needed to either be given a solution y1 or obtain one via some other means such as inspection. Then I could do the y2=u(x)\*y1 and y=y1+y2.

I ended up following the advice I got here for

B: y’’=1+(y’)^2

by setting w=y’. Then I could solve via separation of variables to obtain y=c1 - ln(cos(x+c2))

where c1, c2 are arbitrary constants that could be obtained via initial conditions.

For A,

A: y\*(y’’)=3(y’)^2

I set w=y’, and, this is the part I kept messing up, I saw in my book where they then say that y’’=(dw/dy)_w (as opposed to simply y’’=w’). Once you plug this is and get y_(dw/dy)_w=3_w^2 you can use separation of variables to obtain y=(2x+C)^-0.5 where C is an arbritary constant obtained by using the initial conditions.

Thanks again for the help!

---

<div class="post-metadata">

**Author:** ![MikeS](https://avatars.discourse-cdn.com/v4/letter/m/919ad9/32.png) [@MikeS](https://boards.straightdope.com/u/MikeS)\
**Post date:** [September 10, 2013, 12:25pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/8 "2013-09-10T12:25:49Z")

</div>

> [@dgrdfd](#):
>
> A: y\*(y’‘)=3(y’)^2
> 
> I set w=y’, and, this is the part I kept messing up, I saw in my book where they then say that y’'=(dw/dy)_w (as opposed to simply y’‘=w’). Once you plug this is and get y_(dw/dy)_w=3_w^2 you can use separation of variables to obtain y=(2x+C)^-0.5 where C is an arbritary constant obtained by using the initial conditions.

That’s almost the same as the answer I got. Keep in mind, though, that any second-order ODE should have two arbitrary constants in its general solution, not just the one that you have here. So your answer isn’t quite complete.

---

<div class="post-metadata">

**Author:** ![leahcim](https://avatars.discourse-cdn.com/v4/letter/l/b4bc9f/32.png) [@leahcim](https://boards.straightdope.com/u/leahcim)\
**Post date:** [September 10, 2013, 1:51pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/9 "2013-09-10T13:51:55Z")

</div>

> [@dgrdfd](#):
>
> I saw in my book where they then say that y’‘=(dw/dy)\*w (as opposed to simply y’‘=w’).

I think this is referring to a trick you can do if the ODE is \*homogenous \*(i.e. x does not appear explicitly anywhere in the equation). In that case you can treat the equation as if y was the independent variable, and y’ and y’’ are functions of y.

So when the book is doing that transformation, instead of defining w as being equal to y’ as function of x, w is y’ as a function of y, so derivatives with respect to x have to be converted to derivatives with respect to y through the chain rule, so:

dw/dx = dw/dy \* dy/dx = dw/dy \* w

which explains the weirdness of the transformation.

---

<div class="post-metadata">

**Author:** ![Chronos](https://sea3.discourse-cdn.com/straightdope/user_avatar/boards.straightdope.com/chronos/32/134_2.png) [@Chronos](https://boards.straightdope.com/u/Chronos)\
**Post date:** [September 10, 2013, 2:06pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/10 "2013-09-10T14:06:02Z")

</div>

This is why I prefer to use Liebnitz notation all the way. When you put the "d"s in explicitly, everywhere, everything is so much clearer.

---

<div class="post-metadata">

**Author:** ![leahcim](https://avatars.discourse-cdn.com/v4/letter/l/b4bc9f/32.png) [@leahcim](https://boards.straightdope.com/u/leahcim)\
**Post date:** [September 10, 2013, 5:42pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/11 "2013-09-10T17:42:41Z")

</div>

> [@Chronos](#):
>
> This is why I prefer to use Liebnitz notation all the way. When you put the "d"s in explicitly, everywhere, everything is so much clearer.

But then you run into other problems, like if you are evaluating the function f at two different variables u, and v, do you name the derivative df/du, or df/dv? Lagrange notation resolves that ambiguity.

---

<div class="post-metadata">

**Author:** ![Thudlow\_Boink](https://sea3.discourse-cdn.com/straightdope/user_avatar/boards.straightdope.com/thudlow_boink/32/320_2.png) [@Thudlow\_Boink](https://boards.straightdope.com/u/Thudlow_Boink)\
**Post date:** [September 10, 2013, 6:09pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/12 "2013-09-10T18:09:14Z")

</div>

> [@Chronos](#):
>
> This is why I prefer to use Liebnitz notation all the way.

Spelling nitpick: All the modern sources I’ve seen refer to him as “Leibniz,” but I’ve seen older references to “Leibnitz” (with the t). Upon researching, I found the [New World Encyclopedia](http://www.newworldencyclopedia.org/entry/Gottfried_Leibniz)’s explanation that

> [@](#):
>
> Leibniz was born on July 1, 1646, the child of Friedrich Leibnütz and Catherina Schmuck. He began spelling his name “Leibniz” early in adult life, but others often referred to him as “Leibnitz,” a spelling which persisted until the twentieth century.

IIUC the “z” in German sounds like “tz,” so the pronunciation would be the same either way. However, “ei” and “ie” are quite distinct in German, and Leibniz’s name is unambiguously spelled with “ei.”

> [@leahcim](#):
>
> But then you run into other problems, like if you are evaluating the function f at two different variables u, and v

I’m not sure what kind of situation you’re describing. Can you provide some context where this would be an issue?

---

<div class="post-metadata">

**Author:** ![leahcim](https://avatars.discourse-cdn.com/v4/letter/l/b4bc9f/32.png) [@leahcim](https://boards.straightdope.com/u/leahcim)\
**Post date:** [September 10, 2013, 8:00pm UTC](https://boards.straightdope.com/t/ode-reduction-in-order/668402/13 "2013-09-10T20:00:54Z")

</div>

> [@Thudlow\_Boink](#):
>
> I’m not sure what kind of situation you’re describing. Can you provide some context where this would be an issue?

Say g(x) = f(v(x)) + f(u(x)). Find dg/dx

dg/dx(x) = df/dv (v(x)) \* dv/dx (x) + df/du (u(x)) \* du/dx (x)

But df/dv and df/du here are the same function (the derivative of f with respect to its parameter) but Leibniz notation gives it different names depending on what letter we assign to the parameter. Lagrange notation would call it just f’ in both locations.

In my experience people who want to think of the variables as representing actual physical quantities varying with each other prefer Leibniz notation, while people who want to think of the functions as first class objects and the parameters as arbitrary prefer Lagrange or other notations.
