- /
-
Precision Fractal
on 13 Oct 2021
- 4
- 67
- 0
- 0
- 255
w = 1000;
h = 1000;
N = w*h;
[wx, hy] = meshgrid(1:w, 1:h);
wx = reshape(wx, N, 1);
hy = reshape(hy, N, 1);
rF = hypot(wx - (w + 1)/2, hy - (h + 1)/2);
I = reshape(1 - 2*10^-9*rF.^2 - 2*10^-14*rF.^4, h, w);
figure, imagesc(double(uint16(65535*I)) - 65535*I);
axis image off