def p3d(f, *, nx=-2, px=2,
              ny=-2, py=2,
              D=1, is_log=False, top=3, hi_res=600):

    from matplotlib.cm import get_cmap
    cm = get_cmap("sheerluck")

    def _scale(x):
        if is_log:
            if x <= 0: return infinity
            return log(x)
        else:
            return x / D

    def _arg(x, y):
        val = f(x + I * y)    # f(z)
        a = arg(val) + pi     # f(z).arg
        τ = 2 * pi
        return n(a / τ)       # [-π..π] -> [0..1]

    def _abs(x, y):
        val = f(x + I * y)    # f(z)
        a = val.abs().n()     # f(z).abs
        a = _scale(a)
        if a < top:
            return a
        return infinity

    def _real(x, y):
        val = f(x + I * y)    # f(z)
        a = val.real().n()    # f(z).real
        a = _scale(a)
        if -1 < a < top:
            return a - 0.03
        return infinity

    P = plot3d( _abs,
                (x, nx, px),
                (y, ny, py),
                plot_points=hi_res,
                color=(_arg, cm))

    Q = plot3d( _real,
                (x, nx, px),
                (y, ny, py),
                plot_points=hi_res,
                color="blue",
                opacity=0.5)

    return P + Q


def thats_all_folks():
    P = complex_plot( 1/x^2,
                      (-9.5, 9.5), (-6.5, 6.5),
                      ticks=[list(range(-10, 10)),
                             list(range(-10, 10))],
                      aspect_ratio=1,
                      plot_points=700,
                      contoured=True,
                      figsize=11, cmap="Reds")

    Q = circle((0,0), 2.85, fill=True,
               rgbcolor=(0.0078, 0.1529, 0.2902))

    return P + Q