2.6.4.5 Sympy
from sympy import * 
x = symbols("x") 
V__0 = symbols("V__0") 
c = symbols("c") 
e = symbols("e") 
m = symbols("m") 
n = symbols("n") 
v__0 = symbols("v__0") 
phi = Function("phi") 
ode = Eq(-4*pi*c*n/sqrt(2*e*(-V__0 + phi(x))/m + v__0**2) + Derivative(phi(x), (x, 2)),0) 
ics = {} 
dsolve(ode,func=phi(x),ics=ics)
 
Timed Out
 
Python version: 3.12.3 (main, Aug 14 2025, 17:47:21) [GCC 13.3.0] 
Sympy version 1.14.0