2.6.4.4 Mathematica. Time used: 80.965 (sec). Leaf size: 2754
ode=D[phi[x],{x,2}]==4*Pi*n*c/Sqrt[v0^2+2*e/m*(phi[x]-V0)]; 
ic={}; 
DSolve[{ode,ic},phi[x],x,IncludeSingularSolutions->True]
 

Too large to display