r, theta = symbols("r, theta")
plot_contour(sin(2 * r) * cos(theta), (theta, 0, 2*pi), (r, 0, 7),
    {"levels": 100}, polar_axis=True, aspect="equal")