x,y,t=var("x,y,t")
x0=-10
y0=-6
punto=point((x0,y0), color="red", size="26")
R=10
campo_vectorial=plot_slope_field(4*y,(x,x0-R,x0+R), (y,y0-R,y0+R))

f=function("f")(x)
diffeq =  diff(f,x) == 2*y+1 # Aquí debe sustituir "2*y+1" el lado derecho de la ecuación diferencial correspondiente
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_vectorial+curva_sol+punto).show(ymin = y0-R, ymax = y0+R)

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()