Resolver un PVI con Sympy

 

https://relopezbriega.github.io/blog/2016/01/10/ecuaciones-diferenciales-con-python/

# Defino incognitas
x = sympy.symbols('x')
y = sympy.Function('y')

# Defino la función
f = y(x)**2 + x**2 -1

# Condición inicial
ics = {y(0): 0}

# Resolviendo la ecuación diferencial
edo_sol = sympy.dsolve(y(x).diff(x) - f, ics=ics)
edo_sol