743.sagews

var('x y')
f(x,y) = x^2-y^2
plot3d(f(x,y), (x, -3 ,3), (y, -3, 3))

plot_vector_field(f.gradient(), (x, -3, 3), (y, -3, 3), aspect_ratio=1) + contour_plot(f(x,y), (x, -3, 3), (y, -3, 3), fill=False, contours=9)

g(x, y) = exp(-(x^2+y^2)/4)
plot3d(g(x,y), (x, -3 ,3), (y, -3, 3))

plot_vector_field(g.gradient(), (x, -3, 3), (y, -3, 3), aspect_ratio=1) + contour_plot(g(x,y), (x, -3, 3), (y, -3, 3), fill=False, contours=9)