- /
-
Möbius strip
on 17 Oct 2021
- 3
- 30
- 0
- 0
- 259
n=repmat(linspace(0,2*pi,100),2,1);
t=repmat([-1; 1]/2,1,100);
s=1+t./2.*cos(n/2);
x=s.*cos(n);
y=s.*sin(n);
z=t/2.*sin(n/2);
for i=2:100
P=@(v,j)[v(1,j-1) v(1,j) v(2,j) v(2,j-1)];
C=ones(1,4)*i/100;
patch(P(x,i),P(y,i),P(z,i),C);
hold on;
end
view(-90,20);
camroll(10);
axis equal off;