function moebius(r,g,a,p,f,alpha)
%MOEBIUS(R,G,A,P,F,ALPHA) représente un ruban de Moebius.
% R est le rayon du ruban.
% G est l'inclinaison du ruban au départ.
% A est la largeur du ruban et P est le pas.
% F est le nombre de rotation(s) de la bande sur
% elle-même. F=1,2,3,...
% ALPHA est l'angle de projection.

% jbr, 3-3-97.

s=-a:p:a;
ca=cos(alpha);
sa=sin(alpha);
plot(0,0)
hold on
ra=r+a;
axis([-ra ra -ra ra])
axis('square')
axis('off')
cla
ov=[];
ow=[];
for t=0:p:2*pi+p
	% Coord. des points
	phi=g+t*f/2; %angle de rot. de la dte sur elle-même
	x=(r+s.*cos(phi)).*cos(t);
	y=(r+s.*cos(phi)).*sin(t);
	z=s.*sin(phi);
	v=y;
	w=-x.*sa+z.*ca;
	plot(v,w)
	if ~isempty(ov) 
		for i=1:length(v)
			plot([ov(i) v(i)],[ow(i) w(i)])
		end
	end
	ov=v;
	ow=w;
end
figure(gcf)
hold off
