using AbstractPlotting x = LinRange(-1, 1, 20) y = LinRange(-1, 1, 20) z = x .* y' contour(x, y, z, levels = 0, linewidth = 0, fillrange = true)