using Plots

function Z(n, c, z=c, p=2)
    k = 0
    while (k < n)
        if (abs2(z) > 4) return k end
        k += 1
        z = z^p + c
    end
    return 0
end

const matrix = [i * im + j for i in range(-1.5, stop=1.5, length=1000), j in range(-1.5, stop=1.5, length=1000)]

(a, b) = rand(2)
c = 2(a - .5 + (b - .5) * im)

println(c)
heatmap((x -> Z(60, c, x, i)).(matrix))