- /
-
Torus Onion
on 18 Oct 2021
- 7
- 41
- 2
- 0
- 279
p = pi;
b = [3:5;2.8:4.8;4:6;1.05:.25:1.55];
for a = b
[u,v]=meshgrid(0:.01:2*p/a(4),-p/a(4):.01:p/a(4));
c = cos(u);
d = cos(v);
e = sin(u);
n = (a(1)+a(2)*d);
surf(c.*n,e.*n,-5*log(1-0.3157*v)+a(3)*sin(v)+2*d,'FaceColor','w',EdgeAlpha=.05)
hold on
view([33 11])
end
axis off
lightangle(gca,-45,25)