- /
-
Rotating Newton's Fractal (4th power)
on 29 Nov 2023
- 8
- 73
- 0
- 1
- 164
drawframe(1);
Write your drawframe function below
function drawframe(f)
a=exp(1i*pi*2/47*(f-1));
X=meshgrid(linspace(-2,2,301));
Z=X+1i*X';
for i=1:20
Z = Z-a*(Z.^4-1)./(4*Z.^3);
end
imagesc(X(1,:),X(1,:),angle(Z));
colormap(hsv*0.83);
axis equal off;
end