Transiting planet

Daniel Pereira
on 27 Oct 2021
- 7
- 12
- 0
- 0
- 251
r=meshgrid(0:.1:10);
t=r'*2*pi/10;
P=exp(-r.^4/5);
x=r.*cos(t);
y=r.*sin(t);
contourf(y,x,abs(P),1e3,'edgec','n');
x=rand(10,2)*20-10;
hold on;
plot(x(:,1),x(:,2),'.w');
plot(-1,-1,'ok','markersize',5,'markerfacecolor','k');
colormap hot;
axis equal off;
set(gcf,'color','k');