Contour Function

using AbstractPlotting

r = range(-10, stop = 10, length = 512)
z = ((x, y)-> sin(x) + cos(y)).(r, r')
contour(r, r, z, levels = 5, colormap = :viridis, linewidth = 3)