Homogeneous Functions

Before solving homogeneous differential equations, we must understand what a homogeneous function is.

A function F(x,y)F(x, y) is said to be a homogeneous function of degree nn if, upon replacing xx with λx\lambda x and yy with λy\lambda y (where λ\lambda is any non-zero constant), the factor λ\lambda comes out completely as λn\lambda^n: F(λx,λy)=λnF(x,y)F(\lambda x, \lambda y) = \lambda^n F(x, y)

This means that every term of the function scales in the same way.

Examples:

  • F(x,y)=x2+2xy+y2F(x, y) = x^2 + 2xy + y^2 F(λx,λy)=(λx)2+2(λx)(λy)+(λy)2=λ2(x2+2xy+y2)=λ2F(x,y)F(\lambda x, \lambda y) = (\lambda x)^2 + 2(\lambda x)(\lambda y) + (\lambda y)^2 = \lambda^2(x^2 + 2xy + y^2) = \lambda^2 F(x, y) So this is homogeneous of degree 2.

  • G(x,y)=sin(yx)G(x, y) = \sin\left(\frac{y}{x}\right) G(λx,λy)=sin(λyλx)=sin(yx)=λ0G(x,y)G(\lambda x, \lambda y) = \sin\left(\frac{\lambda y}{\lambda x}\right) = \sin\left(\frac{y}{x}\right) = \lambda^0 G(x, y) So this is homogeneous of degree 0.

  • H(x,y)=sinx+cosyH(x, y) = \sin x + \cos y H(λx,λy)=sin(λx)+cos(λy)H(\lambda x, \lambda y) = \sin(\lambda x) + \cos(\lambda y) Here no single power of λ\lambda can be factored out. So this is not homogeneous.

Homogeneous Differential Equations

A first-order, first-degree differential equation of the form dydx=F(x,y)\frac{dy}{dx} = F(x, y) is called a homogeneous differential equation if F(x,y)F(x, y) is a homogeneous function of degree zero.

That means F(λx,λy)=F(x,y)F(\lambda x, \lambda y) = F(x, y) for every non-zero constant λ\lambda.

Equivalently, such an equation can be rewritten in the form dydx=f(yx)\frac{dy}{dx} = f\left(\frac{y}{x}\right) or, in some cases, more conveniently as dxdy=g(xy).\frac{dx}{dy} = g\left(\frac{x}{y}\right).

This is the key observation that allows us to solve homogeneous equations.

The Solving Algorithm

Homogeneous differential equations are generally not directly separable, but they become separable after a standard substitution.

Method 1: Substituting y=vxy = vx This is used when the equation is written as dydx=f(yx).\frac{dy}{dx} = f\left(\frac{y}{x}\right).

Steps:

  1. Put y=vx,y = vx, where vv is a function of xx.
  2. Differentiate using the product rule: dydx=v+xdvdx.\frac{dy}{dx} = v + x\frac{dv}{dx}.
  3. Substitute y=vxy = vx and dydx=v+xdvdx\frac{dy}{dx} = v + x\frac{dv}{dx} into the equation.
  4. The equation will reduce to a separable equation in vv and xx.
  5. Integrate.
  6. Replace vv by yx\frac{y}{x} to get the final answer.

Method 2: Substituting x=vyx = vy This is useful when it is easier to write the equation as dxdy=g(xy).\frac{dx}{dy} = g\left(\frac{x}{y}\right).

Steps:

  1. Put x=vy,x = vy, where vv is a function of yy.
  2. Differentiate with respect to yy: dxdy=v+ydvdy.\frac{dx}{dy} = v + y\frac{dv}{dy}.
  3. Substitute and simplify.
  4. Separate variables and integrate.
  5. Replace vv by xy\frac{x}{y}.

Important idea: The substitution works because a homogeneous function of degree zero depends only on the ratio yx\frac{y}{x} or xy\frac{x}{y}.

Example 1: Basic Homogeneous Equation

Solve the differential equation: (x+y)dy+(xy)dx=0(x+y)dy + (x-y)dx = 0.

Solution: Step 1: Write the equation in the form dydx=F(x,y)\frac{dy}{dx} = F(x,y). (x+y)dy=(xy)dx(x+y)dy = -(x-y)dx dydx=xyx+y=yxx+y\frac{dy}{dx} = -\frac{x-y}{x+y} = \frac{y-x}{x+y} Now the numerator and denominator are both homogeneous expressions of degree 1, so their ratio is homogeneous of degree 0. Hence the equation is homogeneous.

Step 2: Use the substitution y=vx    dydx=v+xdvdx.y = vx \implies \frac{dy}{dx} = v + x\frac{dv}{dx}. Substitute into the equation: v+xdvdx=vxxvx+x=v1v+1.v + x\frac{dv}{dx} = \frac{vx - x}{vx + x} = \frac{v-1}{v+1}.

Step 3: Isolate the derivative term. xdvdx=v1v+1vx\frac{dv}{dx} = \frac{v-1}{v+1} - v =v1v(v+1)v+1= \frac{v-1-v(v+1)}{v+1} =v1v2vv+1=1+v2v+1.= \frac{v-1-v^2-v}{v+1} = -\frac{1+v^2}{v+1}.

Step 4: Separate variables. v+11+v2dv=dxx.\frac{v+1}{1+v^2}dv = -\frac{dx}{x}.

Step 5: Integrate both sides. Split the left side: v1+v2dv+11+v2dv=dxx.\int \frac{v}{1+v^2}dv + \int \frac{1}{1+v^2}dv = -\int \frac{dx}{x}. So, 12ln(1+v2)+tan1v=lnx+C.\frac{1}{2}\ln(1+v^2) + \tan^{-1}v = -\ln|x| + C.

Step 6: Put v=yxv = \frac{y}{x}. 12ln(1+y2x2)+tan1(yx)=lnx+C.\frac{1}{2}\ln\left(1+\frac{y^2}{x^2}\right) + \tan^{-1}\left(\frac{y}{x}\right) = -\ln|x| + C. Now, 12ln(x2+y2x2)=12ln(x2+y2)lnx,\frac{1}{2}\ln\left(\frac{x^2+y^2}{x^2}\right) = \frac{1}{2}\ln(x^2+y^2) - \ln|x|, so the lnx-\ln|x| on the right combines neatly, giving 12ln(x2+y2)+tan1(yx)=C.\frac{1}{2}\ln(x^2+y^2) + \tan^{-1}\left(\frac{y}{x}\right) = C.

Answer: 12ln(x2+y2)+tan1(yx)=C\frac{1}{2}\ln(x^2 + y^2) + \tan^{-1}\left(\frac{y}{x}\right) = C

Example 2: Verifying and Solving

Show that the differential equation x2dy+y(x+y)dx=0x^2 dy + y(x+y) dx = 0 is homogeneous and solve it.

Solution: Step 1: Rewrite in differential form. x2dydx=y(x+y)x^2\frac{dy}{dx} = -y(x+y) dydx=xy+y2x2.\frac{dy}{dx} = -\frac{xy+y^2}{x^2}. Let F(x,y)=xy+y2x2.F(x,y) = -\frac{xy+y^2}{x^2}. Then, F(λx,λy)=(λx)(λy)+(λy)2(λx)2=λ2(xy+y2)λ2x2=F(x,y).F(\lambda x, \lambda y) = -\frac{(\lambda x)(\lambda y) + (\lambda y)^2}{(\lambda x)^2} = -\frac{\lambda^2(xy+y^2)}{\lambda^2x^2} = F(x,y). So FF is homogeneous of degree 0. Hence the differential equation is homogeneous.

Step 2: Put y=vx    dydx=v+xdvdx.y = vx \implies \frac{dy}{dx} = v + x\frac{dv}{dx}. Substitute: v+xdvdx=(v+v2).v + x\frac{dv}{dx} = -(v+v^2).

Step 3: Simplify. xdvdx=vv2v=(v2+2v).x\frac{dv}{dx} = -v-v^2-v = -(v^2+2v). Thus, dvv(v+2)=dxx.\frac{dv}{v(v+2)} = -\frac{dx}{x}.

Step 4: Integrate. Using partial fractions, 1v(v+2)=12(1v1v+2).\frac{1}{v(v+2)} = \frac{1}{2}\left(\frac{1}{v} - \frac{1}{v+2}\right). Therefore, 12lnvv+2=lnx+C.\frac{1}{2}\ln\left|\frac{v}{v+2}\right| = -\ln|x| + C. Multiply by 2: lnvv+2=2lnx+C1=lnKx2.\ln\left|\frac{v}{v+2}\right| = -2\ln|x| + C_1 = \ln\left|\frac{K}{x^2}\right|. So, vv+2=Kx2.\frac{v}{v+2} = \frac{K}{x^2}.

Step 5: Replace v=yxv=\frac{y}{x}. y/xy/x+2=Kx2\frac{y/x}{y/x+2} = \frac{K}{x^2} yy+2x=Kx2.\frac{y}{y+2x} = \frac{K}{x^2}. Cross-multiplying, x2y=K(y+2x).x^2y = K(y+2x).

Answer: x2y=K(y+2x)x^2 y = K(y + 2x)

Example 3: When to use x=vyx = vy

Solve the differential equation: 2yex/ydx+(y2xex/y)dy=02ye^{x/y} dx + (y - 2xe^{x/y}) dy = 0.

Solution: Step 1: Since the term ex/ye^{x/y} appears, it is more natural to write the equation in terms of dxdy\frac{dx}{dy} and use the substitution x=vyx=vy.

Rearrange: 2yex/ydx=(2xex/yy)dy2ye^{x/y}dx = (2xe^{x/y} - y)dy dxdy=2xex/yy2yex/y.\frac{dx}{dy} = \frac{2xe^{x/y} - y}{2ye^{x/y}}.

Step 2: Put x=vy    dxdy=v+ydvdy.x = vy \implies \frac{dx}{dy} = v + y\frac{dv}{dy}. Substitute: v+ydvdy=2(vy)evy2yev=2vev12ev.v + y\frac{dv}{dy} = \frac{2(vy)e^v - y}{2ye^v} = \frac{2ve^v - 1}{2e^v}.

Step 3: Isolate the derivative term. ydvdy=2vev12evv=12ev.y\frac{dv}{dy} = \frac{2ve^v - 1}{2e^v} - v = -\frac{1}{2e^v}.

Step 4: Separate variables. 2evdv=dyy.2e^v dv = -\frac{dy}{y}.

Step 5: Integrate. 2evdv=dyy\int 2e^v dv = -\int \frac{dy}{y} 2ev=lny+C.2e^v = -\ln|y| + C. Rearranging, 2ev+lny=C.2e^v + \ln|y| = C.

Step 6: Replace v=xyv = \frac{x}{y}. 2ex/y+lny=C.2e^{x/y} + \ln|y| = C.

Answer: 2ex/y+lny=C2e^{x/y} + \ln|y| = C

Example 4: Trigonometric Homogeneous Equation

Solve: xdydx=y+xtan(yx)x \frac{dy}{dx} = y + x \tan\left(\frac{y}{x}\right).

Solution: Step 1: Divide by xx. dydx=yx+tan(yx).\frac{dy}{dx} = \frac{y}{x} + \tan\left(\frac{y}{x}\right). This is of the form f(y/x)f(y/x), so it is homogeneous.

Step 2: Put y=vx    dydx=v+xdvdx.y = vx \implies \frac{dy}{dx} = v + x\frac{dv}{dx}. Substitute: v+xdvdx=v+tanv.v + x\frac{dv}{dx} = v + \tan v.

Step 3: Simplify. xdvdx=tanvx\frac{dv}{dx} = \tan v cotvdv=dxx.\cot v \, dv = \frac{dx}{x}.

Step 4: Integrate. cotvdv=dxx\int \cot v \, dv = \int \frac{dx}{x} lnsinv=lnx+C=lnCx.\ln|\sin v| = \ln|x| + C = \ln|Cx|. So, sinv=Cx.\sin v = Cx.

Step 5: Replace v=yxv = \frac{y}{x}. sin(yx)=Cx.\sin\left(\frac{y}{x}\right) = Cx.

Answer: sin(yx)=Cx\sin\left(\frac{y}{x}\right) = Cx

Example 5: Initial Value Problem (IVP)

Solve the initial value problem: xdydxy+xsin(yx)=0x \frac{dy}{dx} - y + x \sin\left(\frac{y}{x}\right) = 0, given that y=πy = \pi when x=2x = 2.

Solution: Step 1: Rewrite the equation. xdydx=yxsin(yx)x\frac{dy}{dx} = y - x\sin\left(\frac{y}{x}\right) dydx=yxsin(yx).\frac{dy}{dx} = \frac{y}{x} - \sin\left(\frac{y}{x}\right).

Step 2: Put y=vx    dydx=v+xdvdx.y = vx \implies \frac{dy}{dx} = v + x\frac{dv}{dx}. Substitute: v+xdvdx=vsinv.v + x\frac{dv}{dx} = v - \sin v.

Step 3: Simplify and separate. xdvdx=sinvx\frac{dv}{dx} = -\sin v cscvdv=dxx.\csc v \, dv = -\frac{dx}{x}.

Step 4: Integrate. Using cscvdv=lncscvcotv\int \csc v \, dv = \ln|\csc v - \cot v|, lncscvcotv=lnx+C=lnCx.\ln|\csc v - \cot v| = -\ln|x| + C = \ln\left|\frac{C}{x}\right|. Therefore, cscvcotv=Cx.\csc v - \cot v = \frac{C}{x}.

Step 5: Replace v=yxv = \frac{y}{x}. csc(yx)cot(yx)=Cx.\csc\left(\frac{y}{x}\right) - \cot\left(\frac{y}{x}\right) = \frac{C}{x}.

Step 6: Apply the condition x=2x=2, y=πy=\pi. Then yx=π2\frac{y}{x} = \frac{\pi}{2}, so csc(π2)cot(π2)=C2\csc\left(\frac{\pi}{2}\right) - \cot\left(\frac{\pi}{2}\right) = \frac{C}{2} 10=C21 - 0 = \frac{C}{2} C=2.C = 2.

Hence the particular solution is csc(yx)cot(yx)=2x.\csc\left(\frac{y}{x}\right) - \cot\left(\frac{y}{x}\right) = \frac{2}{x}.

Answer: csc(yx)cot(yx)=2x\csc\left(\frac{y}{x}\right) - \cot\left(\frac{y}{x}\right) = \frac{2}{x}

Example 6: Radical Homogeneous Equation

Solve: xdydx=y+x2+y2x \frac{dy}{dx} = y + \sqrt{x^2 + y^2}.

Solution: Step 1: Divide by xx: dydx=yx+x2+y2x.\frac{dy}{dx} = \frac{y}{x} + \frac{\sqrt{x^2+y^2}}{x}. Now, x2+y2x=1+(yx)2\frac{\sqrt{x^2+y^2}}{x} = \sqrt{1+\left(\frac{y}{x}\right)^2} when working in the standard homogeneous setting. Thus, dydx=yx+1+(yx)2.\frac{dy}{dx} = \frac{y}{x} + \sqrt{1+\left(\frac{y}{x}\right)^2}.

Step 2: Put y=vx    dydx=v+xdvdx.y = vx \implies \frac{dy}{dx} = v + x\frac{dv}{dx}. Substitute: v+xdvdx=v+1+v2.v + x\frac{dv}{dx} = v + \sqrt{1+v^2}.

Step 3: Simplify. xdvdx=1+v2.x\frac{dv}{dx} = \sqrt{1+v^2}. So, dv1+v2=dxx.\frac{dv}{\sqrt{1+v^2}} = \frac{dx}{x}.

Step 4: Integrate. Using dv1+v2=lnv+1+v2,\int \frac{dv}{\sqrt{1+v^2}} = \ln|v+\sqrt{1+v^2}|, we get lnv+1+v2=lnx+C=lnCx.\ln|v+\sqrt{1+v^2}| = \ln|x| + C = \ln|Cx|. Therefore, v+1+v2=Cx.v+\sqrt{1+v^2} = Cx.

Step 5: Replace v=yxv = \frac{y}{x}. yx+1+y2x2=Cx.\frac{y}{x} + \sqrt{1+\frac{y^2}{x^2}} = Cx. Multiplying through by xx gives y+x2+y2=Cx2.y + \sqrt{x^2+y^2} = Cx^2.

Answer: y+x2+y2=Cx2y + \sqrt{x^2 + y^2} = Cx^2

Example 7: Implicit General Solution

Solve: (x2+y2)dx2xydy=0(x^2+y^2)dx - 2xy \, dy = 0.

Solution: Step 1: Rearrange. 2xydy=(x2+y2)dx2xy \, dy = (x^2+y^2)dx dydx=x2+y22xy.\frac{dy}{dx} = \frac{x^2+y^2}{2xy}. This is homogeneous because numerator and denominator are both of degree 2.

Step 2: Put y=vx    dydx=v+xdvdx.y = vx \implies \frac{dy}{dx} = v + x\frac{dv}{dx}. Substitute: v+xdvdx=1+v22v.v + x\frac{dv}{dx} = \frac{1+v^2}{2v}.

Step 3: Simplify. xdvdx=1+v22vv=1v22v.x\frac{dv}{dx} = \frac{1+v^2}{2v} - v = \frac{1-v^2}{2v}. Thus, 2v1v2dv=dxx.\frac{2v}{1-v^2}dv = \frac{dx}{x}.

Step 4: Integrate. Let u=1v2u = 1-v^2, so du=2vdvdu = -2v \, dv. Then 2v1v2dv=duu=ln1v2.\int \frac{2v}{1-v^2}dv = -\int \frac{du}{u} = -\ln|1-v^2|. Therefore, ln1v2=lnx+C.-\ln|1-v^2| = \ln|x| + C. So, ln11v2=lnCx\ln\left|\frac{1}{1-v^2}\right| = \ln|Cx| 11v2=Cx.\frac{1}{1-v^2} = Cx.

Step 5: Replace v=yxv = \frac{y}{x}. 11y2/x2=Cx\frac{1}{1-y^2/x^2} = Cx x2x2y2=Cx.\frac{x^2}{x^2-y^2} = Cx. This simplifies to x=C(x2y2).x = C(x^2-y^2). Equivalently, x2y2=Kx,x^2-y^2 = Kx, where K=1CK = \frac{1}{C}.

Answer: x2y2=Kxx^2 - y^2 = Kx

Example 8: Reverse Substitution Case (x=vyx = vy)

Solve: 2xydx+(x2+2y2)dy=02xy \, dx + (x^2 + 2y^2) dy = 0.

Solution: Step 1: It is easier to isolate dxdy\frac{dx}{dy}. 2xydx=(x2+2y2)dy2xy \, dx = -(x^2+2y^2)dy dxdy=x2+2y22xy.\frac{dx}{dy} = -\frac{x^2+2y^2}{2xy}. This is homogeneous, so use x=vyx=vy.

Step 2: Put x=vy    dxdy=v+ydvdy.x = vy \implies \frac{dx}{dy} = v + y\frac{dv}{dy}. Substitute: v+ydvdy=v2+22v.v + y\frac{dv}{dy} = -\frac{v^2+2}{2v}.

Step 3: Simplify. ydvdy=v2+22vv=3v2+22v.y\frac{dv}{dy} = -\frac{v^2+2}{2v} - v = -\frac{3v^2+2}{2v}. So, 2v3v2+2dv=dyy.\frac{2v}{3v^2+2}dv = -\frac{dy}{y}.

Step 4: Integrate. Let u=3v2+2u = 3v^2+2, then du=6vdvdu = 6v \, dv, so 2vdv=du32v \, dv = \frac{du}{3}. Hence, 13duu=dyy\frac{1}{3}\int \frac{du}{u} = -\int \frac{dy}{y} 13ln3v2+2=lny+C.\frac{1}{3}\ln|3v^2+2| = -\ln|y| + C. Multiply by 3: ln3v2+2=3lny+C1=lnKy3.\ln|3v^2+2| = -3\ln|y| + C_1 = \ln\left|\frac{K}{y^3}\right|. Therefore, y3(3v2+2)=K.y^3(3v^2+2) = K.

Step 5: Replace v=xyv = \frac{x}{y}. y3(3x2y2+2)=Ky^3\left(3\frac{x^2}{y^2} + 2\right) = K 3x2y+2y3=K.3x^2y + 2y^3 = K.

Answer: 3x2y+2y3=K3x^2 y + 2y^3 = K

Example 9: Complex Logarithmic Homogeneous Equation

Solve: xdydx=y(lnylnx+1)x \frac{dy}{dx} = y(\ln y - \ln x + 1).

Solution: Step 1: Use log properties. lnylnx=ln(yx).\ln y - \ln x = \ln\left(\frac{y}{x}\right). So the equation becomes xdydx=y(ln(yx)+1)x\frac{dy}{dx} = y\left(\ln\left(\frac{y}{x}\right)+1\right) dydx=yx(ln(yx)+1).\frac{dy}{dx} = \frac{y}{x}\left(\ln\left(\frac{y}{x}\right)+1\right). This is homogeneous.

Step 2: Put y=vx    dydx=v+xdvdx.y = vx \implies \frac{dy}{dx} = v + x\frac{dv}{dx}. Substitute: v+xdvdx=v(lnv+1)=vlnv+v.v + x\frac{dv}{dx} = v(\ln v + 1) = v\ln v + v.

Step 3: Simplify. xdvdx=vlnv.x\frac{dv}{dx} = v\ln v. Thus, dvvlnv=dxx.\frac{dv}{v\ln v} = \frac{dx}{x}.

Step 4: Integrate. Let u=lnvu = \ln v, then du=1vdvdu = \frac{1}{v}dv. So, duu=dxx\int \frac{du}{u} = \int \frac{dx}{x} lnu=lnx+C=lnCx.\ln|u| = \ln|x| + C = \ln|Cx|. Hence, u=Cxu = Cx that is, lnv=Cx.\ln v = Cx.

Step 5: Replace v=yxv = \frac{y}{x}. ln(yx)=Cx.\ln\left(\frac{y}{x}\right) = Cx.

Answer: ln(yx)=Cx\ln\left(\frac{y}{x}\right) = Cx

Example 10: Special Form

Solve: dydx=y2xyx2\frac{dy}{dx} = \frac{y^2}{xy - x^2}.

Solution: Step 1: Rewrite in homogeneous form by dividing numerator and denominator by x2x^2. dydx=(y/x)2(y/x)1.\frac{dy}{dx} = \frac{(y/x)^2}{(y/x)-1}. So the equation is homogeneous.

Step 2: Put y=vx    dydx=v+xdvdx.y = vx \implies \frac{dy}{dx} = v + x\frac{dv}{dx}. Substitute: v+xdvdx=v2v1.v + x\frac{dv}{dx} = \frac{v^2}{v-1}.

Step 3: Simplify. xdvdx=v2v1vx\frac{dv}{dx} = \frac{v^2}{v-1} - v =v2v(v1)v1=vv1.= \frac{v^2 - v(v-1)}{v-1} = \frac{v}{v-1}. So, v1vdv=dxx.\frac{v-1}{v}dv = \frac{dx}{x}.

Step 4: Integrate. (11v)dv=dxx\int \left(1 - \frac{1}{v}\right) dv = \int \frac{dx}{x} vlnv=lnx+C.v - \ln|v| = \ln|x| + C. Now, lnx+lnv=lnvx,\ln|x| + \ln|v| = \ln|vx|, so we may write v=lnvx+C.v = \ln|vx| + C.

Step 5: Replace v=yxv = \frac{y}{x}. Since vx=yvx = y, yx=lny+C.\frac{y}{x} = \ln|y| + C.

Answer: yx=lny+C\frac{y}{x} = \ln|y| + C