restart: with(linalg): with(plots):line:=[R,0,u];T:=matrix([[1,0,0,-R],[0,1,0,0],[0,0,1,0],[0,0,0,1]]);Ry:=matrix([[cos(v/2),0,sin(v/2),0],[0,1,0,0],[-sin(v/2),0,cos(v/2),0],[0,0,0,1]]);Rz:=matrix([[cos(v),-sin(v),0,0],[sin(v),cos(v),0,0],[0,0,1,0],[0,0,0,1]]);M:=evalm(Rz&*inverse(T)&*Ry&*T);mobius:=[seq(convert(evalm(M&*[seq(line[i],i=1..3),1]),list)[i],i=1..3)];R:=2: plotline:=spacecurve(line,u=-1..1,color=green,thickness=2): plotmobius:=plot3d(mobius,u=-1..1,v=0..2*Pi,color=grey):display(plotline,plotmobius,scaling=constrained,axes=normal);