Here is my solution, there is one problem that i do not understand and hope
another can explain. When col=20 W is filled with NaN + NaN*im but for col=10
it work like matlab, I think it is because of overflow.
using PyPlot
col = 10;
m = 400;
cx = -.6;
cy = 0;
l = 1.5;
x = linspace(cx-l,cx + l,m);
y = linspace(cy-l,cy + l,m);
X = repmat(x,1,m)
Y = repmat(y,1,m)'
Z = zeros(m,m);
C = X + im * Y;
W = zeros(Z)
for k = 1:col;
Z = Z.*Z+C
W = exp(-abs(Z));
println(Z[100,100])
end
pcolor(W,cmap=ColorMap("copper"))
By changing the calculation to happen with BigFloats it works.
Z = zeros(BigFloat, m,m);
pcolor(float64(W),cmap=ColorMap("copper"))
