- /
-
Coronal mass ejection
on 15 Oct 2021
- 1
- 24
- 0
- 0
- 259
x=[0.1;0.1];
c=hot;
rng(1234)
p=2*randn(1,3);
y=x;
for i=1:3e5
z(1)=p(1)+p(2)*sin(y(1)*y(2));
z(2)=p(3)*(cos(y(1)-y(2)))+real(log(y(2)*y(1)^3));
if i>1000
x(:,i-1000)=z;
end
y=z;
end
scatter(x(1,:),x(2,:),1,x(1,:).*x(2,:),'filled');set(gca,'visible','off');colormap(c(256:-1:1,:))