x,y=var("x,y")
x0=-10
y0=-6
punto=point((x0,y0), color="red", size="22")
R=10
campo_vec=plot_slope_field(y,(x,x0-R,x0+R), (y,y0-R,y0+R),headaxislength=3, headlength=3)

f=function("f")(x)
diffeq =  diff(f,x)==f
show(diffeq)
f(x)=desolve(diffeq, [f,x], [x0,y0])
show(f(x))
curva_sol=plot(f, (x,x0-R,x0+R), color="blue")
(campo_vec+curva_sol+punto).show(ymin = y0-R, ymax = y0+R)

$$\newcommand{\Bold}[1]{\mathbf{#1}}D[0]\left(f\right)\left(x\right) = f\left(x\right)$$

$$\newcommand{\Bold}[1]{\mathbf{#1}}-6 \, e^{\left(x + 10\right)}$$

plot_vector_field( (sin(y), y-x**2), (x,-3,3), (y,-1,8) )

u,v=var("u,v")
u0=0
v0=0
R=10
campo_vec=plot_vector_field( (v,-4*u), (u,u0-R,u0+R), (v,v0-R,v0+R) )
campo_vec.show()

u,v=var("u,v")
u0=0
v0=0
R=10
campo_vec=plot_vector_field( (v,4*u), (u,u0-R,u0+R), (v,v0-R,v0+R) )
campo_vec.show()