% Dibuix d'una banda de Moebius
% (codi amb Matlab i apte també per a Octave)
% Geometria Diferencial, FME, UPC
% 2021/02/12

% Paràmetres de la banda de Moebius
r = 11; % radi de la banda
t = linspace(0, 2*pi, 101);
zor = linspace(-1, 1, 21); % ample de la banda

% Taula de punts de la banda
for k = 1:length(t)
  A = [
    cos(t(k)) -sin(t(k)) 0
    sin(t(k)) cos(t(k)) 0
    0 0 1
  ];
  B = [
    1 0 0
    0 cos(t(k)/2) -sin(t(k)/2)
    0 sin(t(k)/2) cos(t(k)/2)
  ];

  for l = 1:length(zor)
    phi = A*(B*[0; r; zor(l)] + [0; r; 0]);
    x(k, l) = phi(1);
    y(k, l) = phi(2);
    z(k, l) = phi(3);
  end
end

% Dibuixem la superfície
%xlabel('x');
%ylabel('y');
%zlabel('z');
surf(x, y, z);
print("01.png", "-dpng")
%axis equal