Mandelbrot set: Difference between revisions

m
Line 5,289:
 
Z, dZ = zero(C), zero(C)
TS, ST, D = zeros(size(C)), zeros(size(C)), zeros(size(C))
 
for k = 1:n
M = abs.(Z) .< r
TS[M], ST[M] = T[M] .+ 1, S[M] + exp.(-abs.(Z[M])), T[M] .+ 1
Z[M], dZ[M] = Z[M] .^ 2 + C[M], 2 * Z[M] .* dZ[M] .+ 1
end
 
heatmap(TS .^ 0.1, c=:jet)
savefig("Mandelbrot_set_1.png")
 
heatmap(ST .^ 0.1, c=:jet)
savefig("Mandelbrot_set_2.png")
 
Line 5,329:
 
Z, dZ = zero(C), zero(C)
TS, ST, D = zeros(size(C)), zeros(size(C)), zeros(size(C))
 
for k = 1:n
M = abs.(Z) .< r
TS[M], ST[M] = T[M] .+ 1, S[M] + exp.(-abs.(Z[M])), T[M] .+ 1
Z[M], dZ[M] = Z[M] .^ 2 + C[M], 2 * Z[M] .* dZ[M] .+ 1
end
 
heatmap(-TS' .^ 0.1, c=:nipy_spectral)
savefig("Mercator_map_1.png")
 
heatmap(-ST' .^ 0.1, c=:nipy_spectral)
savefig("Mercator_map_2.png")
 
305

edits