> 
ode:=diff(y(x),x$2)=(diff(y(x),x))^3(diff(y(x),x))^2;
IC:=y(0)=3,D(y)(0)=1;

Error, (in dsolve) numeric exception: division by zero
> 
sol:=y(x)=3+x;
odetest(sol,[ode,IC])

> 
DEtools:odeadvisor(ode);

Methods for second order ODEs:
 Trying classification methods 
trying 2nd order Liouville
trying 2nd order WeierstrassP
trying 2nd order JacobiSN
differential order: 2; trying a linearization to 3rd order
 trying a change of variables {x > y(x), y(x) > x}
differential order: 2; trying a linearization to 3rd order
trying 2nd order ODE linearizable_by_differentiation
trying 2nd order, 2 integrating factors of the form mu(x,y)
trying differential order: 2; missing variables
> Computing canonical coordinates for the symmetry [0, 1]
> Rewriting ODE in canonical coordinates by means of differential invariants
Try computing 1 more symmetries for ODE written in canonical coordinates
> Computing symmetries using: way = 3
Found another symmetry:
Found another symmetry:
Computing a convenient ordering to use the 3 symmetries available
< differential order: 2; canonical coordinates successful
< differential order 2; missing variables successful
Error, (in dsolve) numeric exception: division by zero
Methods for second order ODEs:
 Trying classification methods 
trying 2nd order Liouville
trying 2nd order WeierstrassP
trying 2nd order JacobiSN
differential order: 2; trying a linearization to 3rd order
 trying a change of variables {x > y(x), y(x) > x}
differential order: 2; trying a linearization to 3rd order
trying 2nd order ODE linearizable_by_differentiation
trying 2nd order, 2 integrating factors of the form mu(x,y)
trying differential order: 2; missing variables
> Computing canonical coordinates for the symmetry [0, 1]
> Rewriting ODE in canonical coordinates by means of differential invariants
Try computing 1 more symmetries for ODE written in canonical coordinates
> Computing symmetries using: way = 3
Found another symmetry:
Found another symmetry:
Computing a convenient ordering to use the 3 symmetries available
> Calling odsolve with the ODE diff(_b(_a) _a) = _b(_a)^3_b(_a)^2 _b(_a) HINT = [[1 0] [_a+y _b*(_b1)]]
*** Sublevel 2 ***
symmetry methods on request
1st order, trying reduction of order with given symmetries:
1st order, trying the canonical coordinates of the invariance group
< 1st order, canonical coordinates successful
< differential order: 2; canonical coordinates successful
< differential order 2; missing variables successful
