- /
-
Flower or Fractal
on 5 Oct 2021
- 1
- 9
- 0
- 0
- 174
c = jet(4);
x(1)=0.1;y(1)=0;
a = .02;
b = .97;
for k=1:5e5
x(k+1)=y(k)+a*x(k)-20/(x(k)^3+1)+1+0.2*exp(-y(k)^2);
y(k+1)=-b*x(k);
end
scatter(x+y,x-y,5,c(2,:),"filled")
ax=gca;
ax.Visible='off';