ode:=diff(diff(phi(x),x),x) = 4*Pi*n*c/(v__0^2+2*e/m*(phi(x)-V__0))^(1/2); dsolve(ode,phi(x), singsol=all);
Maple trace
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 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 symmetries using: way = 3 -> Calling odsolve with the ODE, diff(_b(_a),_a)*_b(_a)-4*Pi*n*c/(-(-m*v__0^2+2*V__0*e-2*_a*e)/m)^(1/2) = 0, _b(_a), HINT = [[-2/3*(-m*v__0^2+2*V__0*e-2*_a*e)/e, 1/3*_b]] *** Sublevel 2 *** symmetry methods on request 1st order, trying reduction of order with given symmetries: [-2/3*(-m*v__0^2+2*V__0*e-2*_a*e)/e, 1/3*_b] 1st order, trying the canonical coordinates of the invariance group -> Calling odsolve with the ODE, diff(y(x),x) = 1/2*y(x)/(2*(x-V__0)*e+v__0^2*m)*e, y(x) *** Sublevel 3 *** Methods for first order ODEs: --- Trying classification methods --- trying a quadrature trying 1st order linear <- 1st order linear successful <- 1st order, canonical coordinates successful <- differential order: 2; canonical coordinates successful <- differential order 2; missing variables successful