contours

using AbstractPlotting

x = LinRange(-1, 1, 20)
y = LinRange(-1, 1, 20)
z = x .* y'

vbox(
    contour(x, y, z, levels = 50, linewidth =3),
    contour(x, y, z, levels = 0, linewidth = 0, fillrange = true)
)