2.1.65 Problem 65

Solved as second order missing x ode
Maple
Mathematica
Sympy

Internal problem ID [8777]
Book : Own collection of miscellaneous problems
Section : section 1.0
Problem number : 65
Date solved : Sunday, March 30, 2025 at 01:32:59 PM
CAS classification : [[_2nd_order, _quadrature]]

Solved as second order missing x ode

Time used: 0.125 (sec)

Solve

y2y=0

This is missing independent variable second order ode. Solved by reduction of order by using substitution which makes the dependent variable y an independent variable. Using

y=p

Then

y=dpdx=dpdydydx=pdpdy

Hence the ode becomes

y2p(y)(ddyp(y))=0

Which is now solved as first order ode for p(y).

Since the ode has the form p=f(y), then we only need to integrate f(y).

dp=0dy+c1p=c1

For solution (1) found earlier, since p=y then we now have a new first order ode to solve which is

y=c1

Since the ode has the form y=f(x), then we only need to integrate f(x).

dy=c1dxy=c1x+c2

Will add steps showing solving for IC soon.

Summary of solutions found

y=c1x+c2

Maple. Time used: 0.002 (sec). Leaf size: 13
ode:=y(x)^2*diff(diff(y(x),x),x) = 0; 
dsolve(ode,y(x), singsol=all);
 
y=0y=c1x+c2

Maple trace

Methods for second order ODEs: 
--- Trying classification methods --- 
trying a quadrature 
<- quadrature successful
 

Maple step by step

Let’s solvey(x)2(ddxddxy(x))=0Highest derivative means the order of the ODE is2ddxddxy(x)Isolate 2nd derivativeddxddxy(x)=0Characteristic polynomial of ODEr2=0Use quadratic formula to solve forrr=0±(0)2Roots of the characteristic polynomialr=01st solution of the ODEy1(x)=1Repeated root, multiplyy1(x)byxto ensure linear independencey2(x)=xGeneral solution of the ODEy(x)=C1y1(x)+C2y2(x)Substitute in solutionsy(x)=C2x+C1
Mathematica. Time used: 0.003 (sec). Leaf size: 17
ode=y[x]^2*D[y[x],{x,2}]==0; 
ic={}; 
DSolve[{ode,ic},y[x],x,IncludeSingularSolutions->True]
 
y(x)0y(x)c2x+c1
Sympy. Time used: 0.150 (sec). Leaf size: 7
from sympy import * 
x = symbols("x") 
y = Function("y") 
ode = Eq(y(x)**2*Derivative(y(x), (x, 2)),0) 
ics = {} 
dsolve(ode,func=y(x),ics=ics)
 
y(x)=C1+C2x