- /
-
Gromit's Rocket
on 25 Oct 2021
- 6
- 14
- 0
- 0
- 273
L=@linspace;
x = linspace(-1.2,1,400); y=x;
[X,Y]=meshgrid(x, y);
C=X+i*Y;
m=5*ones(size(C));
z = C; %Inital Cond. z gets C
%n=1; %Delta brot, intro n
for n= 0:5
f=asin(C.^(tan(((conj((z)+C.^2))+C.^2).^n)));
z=atanh(sqrt(((C.^(sqrt(asinh(tan(log10((f.^f).^n)))))).^C)+C.^2));%
j=abs(z)>2;
z(j)=NaN; %As xpl
m(j)=n;
end
imagesc(m')
colormap hot %imshow(m,[])
axis equal off