- /
-
Bessel pentagon
on 20 Oct 2021
- 2
- 21
- 2
- 0
- 210
[x,y]=meshgrid(-10:.01:10);
b=@(p,q) -besselj(0,(((x-p*2).^2+(y-q*2).^2).^0.5));
mesh(x,y,b(2,0)+b(0.618,1.9021)+b(-1.618,1.1756)+b(-1.618,-1.1756)+b(0.618,-1.9021))
colormap(gca,'colorcube')
view(0,90)
axis equal off