- /
-
Rosette 1313
on 17 Oct 2021
- 11
- 98
- 4
- 0
- 280
[X,Y]=meshgrid(-3:5e-3:3);
a=0.1+0.02i;
fz=sqrt(abs((X+i*Y).^5+(X-i*Y).^5+a*((X+i*Y).^6.*(X-i*Y)+(X+i*Y).*(X-i*Y).^6)+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 interp
axis([200 900 300 999],'square','off')