- /
-
Snail Shell
on 19 Oct 2022
- 8
- 59
- 3
- 0
- 225
u = linspace(0,2,512);
v = u.';
s = (1-v/6)/3 ;
t = .2-s.*(1+cospi(u));
x = t.*cospi(2*v);
y = t.*sinpi(2*v);
z = v/4-s.*sinpi(u);
surf(x,y,z,y)
shading interp
axis off equal
colormap(parula)
lightangle(80,-40)
lightangle(-90,60)
view([-150 10])