Substitution in a differential equation


By Hal 2001 (P3046) on Sunday, October 22, 2000 - 11:58 pm :

I am doing some work on integration and differential equations. I got stumped by the following problems.

Q1: Show that the substitution y = vx, where v is a function of x, transforms the differential equation
xy(dy/dx) = x 2 + y2
into the differential equation vx(dv/dx) = 1.
Then it says, "Hence solve the original differential eqn given that y=2 at x=1" - this part is fine, I can do, no probs. But the part that I stated before this is the problem.

Q2: The normal at every point on a given curve passes through the point (2,3). Prove that the eqn of the curve satisfies the differential eqn:
(x-2) + (dy/dx)(y-3) = 0

Hence find the cartesian eqn of the family of curves which satisfy the condition above. (I think I can do this part)

Find the eqn of the perpendicular curve which passes through the point (4,2).

Hope someone can help.
Thanks in advance.
HAL2001
By Kerwin Hui (Kwkh2) on Monday, October 23, 2000 - 09:27 am :

Hal,

If y=vx, we have dy/dx=v+x(dv/dx). Substituting, we have

x(vx)(v+x(dv/dx))=x2 +(vx)2

Now cancel the common factor x2 and a little bit more rearranging gives you the required form.

Now, for the second question. Do you know the general formula for the equation of normal? The equation of normal is

(dy/dx)(x0,y0) (y-y0 )+(x-x0 )=0

And for the very last bit, two curves y=y1 (x) and y=y2 (x) are perpendicular if and only if

(dy1 /dx)(dy2 /dx)=-1

at the point(s) of intersection.

Hope this helps.

Kerwin


By Hal 2001 (P3046) on Monday, October 23, 2000 - 01:30 pm :

Hi Kerwin,

Thanks for the quick response.

The problem for me to solving Q1 was that I was not sure how exactly to tackle y=vx. I see that you have done this : (dy/dx)= v + x(dv/dx). I however don't completely understand how you have implicitly differentiated this. Why does the 'vx' part become 'v + x(dv/dx)'. Is this becaise v is a function of x. I think this part would be better understood if I fully comprehended the meaning of the line 'v is a function of x'. Am I right in thinking that x is a function of y, thats why you get the LHS of y ---> (dy/dx), and v is a function of x, so vx are a product so it becomes x/(dv/dx). I am I on the right line of thinking?

The second part of the Q1 is posing me more of a problem than I initially thought. As the original differential eqn is xy(dy/dx)= x 2 + y 2 , I am not fully aware of what to do with the with any of the sides to get the x's and y's to either side of the eqn. It got very messy when I tried rearranging it. Here is how far I got:

Originally we have: xy(dy/dx) = x 2 + y2
Then lets try to: y(dy/x) = x + (y 2 /x)
And divide through by y 2 : y -1 (dy/dx) = (x/y) + (1/x)
Then (1/y)(1-x)(dy/dx) = (1/x)

So: ò(1/y)(dy)=ò((1-x)/(x) (dx)
Hence: lny = x-1 - 1 + C


I am nearly sure that the above steps are wrong.
Any ideas?

Thanks again.
HAL2001
By Hal 2001 (P3046) on Monday, October 23, 2000 - 01:26 pm :

Regarding question 2:
The general form that you showed I think I know in the form, (note m=gradient), so (y-y1 ) = (1/m)(x-x1 ), then rearrange you get m1 (y-y1 ) -(x-x1 ).

But how does that answer the question in proving it passes the point (2,3). Even though the point initially takes was (2,0). How do you put it into writing as an answer.

I did not get the right answer in the cartesian form either.

Regards
HAL2001


By Kerwin Hui (Kwkh2) on Monday, October 23, 2000 - 02:00 pm :
Hal,

For the first part, we differentiate y=v x wrt x. We have not made the assumption that v is a constant, so we have

dy/dx=d(v x)/dx=v+x(dv/dx)

by the product rule.

The second part follows from the first part. Having got to the equation

v x(dv/dx)=1

we can separate the variables and integrate, giving

òv dv=òdx/x

v2 = ln(Ax2)

giving you v in terms of x when you evaluate A. Substituting v=y/x gives you y in terms of x.

Hopefully this is clear.

Kerwin


By Hal 2001 (P3046) on Monday, October 23, 2000 - 02:32 pm :

Thanks Kerwin!

Its much clearer now.
The trick was not to assume that v was a constant and diff wrt x. And treat the vx with the product rule.

The second part of the question, the reason I got thrown off the path was that the question said to use the original differential eqn, thus my goofed up attempt. The eqn you used works.

Thanks once again.
HAL2001


By Kerwin Hui (Kwkh2) on Monday, October 23, 2000 - 02:57 pm :

Hal,

Now for the second question, the normal to the curve has gradient -(dx/dy) evaluated at x=x0 and y=y0 ., where (x0 ,y0 ) is the point on the curve.

so we have, assuming y=y1 (x), gradient of normal is -1/(dy1 /dx)(x0,y0) passing through the point (x0 ,y0 ) on the curve. So we have

(y-y0 )=-1/(dy1 /dx)(x0,y0) (x-x0 )

Rearranging gives (dy1 /dx)(x0,y0) (y-y0 )+(x-x0 )=0.

Since the normal passes through (2,3), we have

(dy1 /dx)x0,y0 (3-y0 )+(2-x0 )=0

Now, since this properites is true for all point on the curve, we can replace x0 by x and y0 by y, giving the required answer as y1 =y.

Kerwin


By Hal 2001 (P3046) on Monday, October 23, 2000 - 05:41 pm :

Thanks for that further explanation Kerwin.

I don't quite understand why the last line of the post turned out like (3-y 0 )+(2-x 0 ), I thought that it would be the other way round like, m 2 (y-3) + (x-2)=0

I also got confused with the notation in using subscript 1s and 0s generally. What does the line 'assuming y=y1 (x)' mean?

I think I get the general drift of the solution that you proposed. Though it has not quite completly hit me to be instantly obvious, yet.

Regards
HAL2001


By Kerwin Hui (Kwkh2) on Monday, October 23, 2000 - 05:55 pm :

Hal,

Well, the last line should read

(dy1 /dx)(x0,y0) (3-y0 )+(2-x0 )=0,

which when multiplied both side by -1 and replacing x0 by x, y0 and y1 by y gives the required solution.

The reason I used y1 (x) is completely arbitrary. You can use simply y(x) if you wanted. Using y1 (x) has the advantage on the second part, when I can use y2 (x) for the curve normal to y1 (x).

Hope this explains it.

Kerwin


By Hal 2001 (P3046) on Monday, October 23, 2000 - 06:53 pm :

Thanks Kerwin,

I see what you mean now.
The m1 x m2=-1.

Once again, thanks for your help.
HAL2001