How to solve this second-order nonlinear ordinary differential equation?

dimension10

I was trying to solve a calculus problem (about the equation of the curve traced by a dog chasing a postman) and I came across the following equation. I would like to know how to solve it.

$$\frac{w}{v}(1-\frac{y(x)y''(x)}{{y'(x)}^{2}})=\sqrt{1+{y(x)}^{2}}-\sqrt{1+{C}^{2}}$$

Thanks.

Related Differential Equations News on Phys.org

jackmell

I was trying to solve a calculus problem (about the equation of the curve traced by a dog chasing a postman) and I came across the following equation. I would like to know how to solve it.

$$\frac{w}{v}(1-\frac{y(x)y''(x)}{{y'(x)}^{2}})=\sqrt{1+{y(x)}^{2}}-\sqrt{1+{C}^{2}}$$

Thanks.
You have got to remove all that clutter. It just makes it way more confussing. And what's with all those extra letters. Need to absorb them into one another for now. How about start with writing it as:

$$k\left(1-\frac{yy''}{(y')^2}\right)=\sqrt{1+y^2}-c$$

See what I mean?
Now, itsn't that an equation with the independent variable missing? You know the kind where you let:

$$\frac{dy}{dx}=p$$
$$\frac{d^2y}{dx^2}=p\frac{dp}{dy}$$

Now, make all those substitutions, solve it for p in terms of y, and then integrate one more time.

Last edited:

dimension10

You have got to remove all that clutter. It just makes it way more confussing. And what's with all those extra letters. Need to absorb them into one another for now. How about start with writing it as:

$$k\left(1-\frac{yy''}{(y')^2}\right)=\sqrt{1+y^2}-c$$

See what I mean?
Now, itsn't that an equation with the independent variable missing? You know the kind where you let:

$$\frac{dy}{dx}=p$$
$$\frac{d^2y}{dx^2}=p\frac{dp}{dy}$$

Now, make all those substitutions, solve it for p in terms of y, and then integrate one more time.
But I tried that and got

$$p=\int {(\frac{dy}{dx})}^{2}(\frac{c-\sqrt{1+{y}^{2}}}{k}+1)dx$$

which is equivalent to

$$y=\int \int {(\frac{dy}{dx})}^{2}(\frac{c-\sqrt{1+{y}^{2}}}{k}+1)dx\; dx$$

but I have no idea how to solve that...

jackmell

Let's solve this one:

$$k\left(1-\frac{yy''}{(y')^2}\right)=\sqrt{1+y^2}-c,\quad y(0)=1,\quad y'(0)=1$$
Letting:

$$\frac{dy}{dt}=p,\quad \frac{d^2y}{dy^2}=p\frac{dp}{dy}$$

then:

$$k\left(1-\frac{ypp'}{p^2}\right)=\sqrt{1+y^2}-c$$

at this point, separate the variables first, then let $u=k+c$ and then integrate to obtain:

$$\int_{1}^{y'} k\frac{dp}{p}=\int_0^y \frac{u-\sqrt{1+y^2}}{y}dy$$

$$k\log(y')=\left(u\log(y)+\log(1+\sqrt{1+y^2})-\log(y)-\sqrt{1+y^2}\right)-k2$$

so that:

\begin{align*} \frac{dy}{dt}&=\exp\left\{\frac{1}{k}\biggr[u\log(y)+\log(1+\sqrt{1+y^2}-\log(y)-\sqrt{1+y^2}-k2\biggr]\right\}\\ &=e^{f(y)} \end{align*}

then we integrate again:

$$\int_{1}^{y} e^{-f(y)}dy=\int_0^t dt$$

to obtain the implicit answer:

$$h(y)=t$$

or:

$$y=h^{-1}(t)$$

Now, ain't that beautiful or what?

and we can check this numerically in Mathematica by first solving the IVP numerically, then integrate our analytic solution numerically and compare the two. Note in the Table command below, I implicitly "invert" the function h(y) by reversing the calculated values {t(y),y}:

Code:
w = 0.5;
v = 1.7;
k = w/v;
c2 = 1.5;
c = Sqrt[1 + c2^2];
y0 = 1;
y1 = 1;
u = k + c;
k2 = Log[1 + Sqrt] - Sqrt

mysol = NDSolve[{k*(1 - (y[t]*Derivative[y][t])/Derivative[y][t]^2) == Sqrt[1 + y[t]^2] - c,
y == y0, Derivative[y] == y1}, y, {t, 0, 5}]

p1 = Plot[y[t] /. mysol, {t, 0, 5}, PlotStyle -> Red]

f[y_] := Exp[(-k^(-1))*(u*Log[y] + Log[1 + Sqrt[1 + y^2]] - Log[y] - Sqrt[1 + y^2] - k2)]

myt[y_] := NIntegrate[f[s], {s, y0, y}]

mysol = Table[{myt[y], y}, {y, y0, 5, 0.1}];

p2 = ListPlot[mysol, Joined -> True, PlotStyle -> Blue]
Show[{p1, p2}]
There may however be a better way to solve this.

Last edited:

dimension10

Let's solve this one:

$$k\left(1-\frac{yy''}{(y')^2}\right)=\sqrt{1+y^2}-c,\quad y(0)=1,\quad y'(0)=1$$
Letting:

$$\frac{dy}{dt}=p,\quad \frac{d^2y}{dy^2}=p\frac{dp}{dy}$$

then:

$$k\left(1-\frac{ypp'}{p^2}\right)=\sqrt{1+y^2}-c$$

at this point, separate the variables first, then let $u=k+c$ and then integrate to obtain:

$$\int_{1}^{y'} k\frac{dp}{p}=\int_0^y \frac{u-\sqrt{1+y^2}}{y}dy$$

$$k\log(y')=\left(u\log(y)+\log(1+\sqrt{1+y^2})-\log(y)-\sqrt{1+y^2}\right)-k2$$

so that:

\begin{align*} \frac{dy}{dt}&=\exp\left\{\frac{1}{k}\biggr[u\log(y)+\log(1+\sqrt{1+y^2}-\log(y)-\sqrt{1+y^2}-k2\biggr]\right\}\\ &=e^{f(y)} \end{align*}

then we integrate again:

$$\int_{1}^{y} e^{-f(y)}dy=\int_0^t dt$$

to obtain the implicit answer:

$$h(y)=t$$

or:

$$y=h^{-1}(t)$$

Now, ain't that beautiful or what?

and we can check this numerically in Mathematica by first solving the IVP numerically, then integrate our analytic solution numerically and compare the two. Note in the Table command below, I implicitly "invert" the function h(y) by reversing the calculated values {t(y),y}:

Code:
w = 0.5;
v = 1.7;
k = w/v;
c2 = 1.5;
c = Sqrt[1 + c2^2];
y0 = 1;
y1 = 1;
u = k + c;
k2 = Log[1 + Sqrt] - Sqrt

mysol = NDSolve[{k*(1 - (y[t]*Derivative[y][t])/Derivative[y][t]^2) == Sqrt[1 + y[t]^2] - c,
y == y0, Derivative[y] == y1}, y, {t, 0, 5}]

p1 = Plot[y[t] /. mysol, {t, 0, 5}, PlotStyle -> Red]

f[y_] := Exp[(-k^(-1))*(u*Log[y] + Log[1 + Sqrt[1 + y^2]] - Log[y] - Sqrt[1 + y^2] - k2)]

myt[y_] := NIntegrate[f[s], {s, y0, y}]

mysol = Table[{myt[y], y}, {y, y0, 5, 0.1}];

p2 = ListPlot[mysol, Joined -> True, PlotStyle -> Blue]
Show[{p1, p2}]
There may however be a better way to solve this.
Thanks.

Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving