r, theta = symbols("r theta")
p = plot_parametric_region(r * cos(theta), r * sin(theta),
    (r, 1, 2), (theta, 0, 2*pi/3),
    {"color": "k", "linewidth": 0.75},
    n1=5, n2=15, aspect="equal")  # doctest: +SKIP