Moebius

t=linspace(-1,1,20)';
x=linspace(0,%pi,40);
factor=2+ t*cos(x);
X=factor*diag(cos(2*x));
Y=factor*diag(sin(2*x));
Z=t*sin(x);
plot3d2(X,Y,Z);