710.sagews

# desplegar grafica de campo direccional (ecuaciones diferenciales)
x,y=var('x,y')
f(x,y) = exp(x*y)
f(x,y) = exp(-(x^2+y^2))
grad_f = [diff(f,var) for var in (x,y)]
plot_vector_field(grad_f,(-2,2),(-2,2))