# definimos nuestras variables
x,y=var("x,y") 
# definimos las condiciones iniciales
x0=-10 
y0=-6
#creamos el punto asociado a la condición inicia
punto=point((x0,y0), color="red", size="26")
#R es el ancho del cuadrado del punto inicial, alrededor del cual graficaremos
R=10
#creamos el campo vectorial correspondiente a nuestra ecuación diferencial
#en este ejemplo, la ecuación es y'=2*y+1
#aquí debe sustituir "2*y+1" el lado derecho de la ecuación diferencial correspondiente
campo_vectorial=plot_slope_field(2*y+1,(x,x0-R,x0+R), (y,y0-R,y0+R))

#ahora vamos a resolver la ecuación diferencia
#creamos una función abstacta "f" dependiente de x
f=function("f")(x)
#definimos la ecuación diferencia 
#observe que ahora usamos la función f en lugar de la variable y
#diff(f,x) es f'(x)
diffeq =  diff(f,x) == 2*f+1 
#resolvemos el problema de valores iniciales con el método desolve
f(x)=desolve(diffeq, [f,x], [x0,y0])
#mostramos en pantalla la solucion
show("solucion f(x)=",f(x))
#creamos la isoclina correpondiente
curva_sol=plot(f, (x,x0-R,x0+R), color="blue")

#unimos todas nuestras graficas en una sola
grafica=campo_vectorial+curva_sol+punto
#y finalmente la mostramos en pantalla
grafica.show(ymin = y0-R, ymax = y0+R)

$$\newcommand{\Bold}[1]{\mathbf{#1}}\verb|solucion|\phantom{\verb!x!}\verb|f(x)=| -\frac{11}{2} \, e^{\left(2 \, x + 20\right)} - \frac{1}{2}$$