- /
-
Rosette 2253
on 21 Oct 2021
- 8
- 49
- 2
- 0
- 278
[X,Y]=meshgrid(-3:5e-3:3);
a=0.1+0.03i;
fz=sqrt(abs((X+i*Y).^4+(X-i*Y).^5+a*((X+i*Y).^5.*(X-i*Y)+(X+i*Y).*(X-i*Y).^7)+a*((X+i*Y).^4./(X-i*Y).^6+(X-i*Y).^4./(X+i*Y).^6)));
%fz(fz>9)=9;
pcolor(fz);
cm=colormap("lines");
cm(end,:)=cm(end-1,:);
colormap(cm);
shading flat
axis([200 900 100 900],'square','off')