restart:
with(plots):with(plottools):
with(LinearAlgebra):
<Text-field style="Heading 1" layout="Heading 1">2d</Text-field>
A:=Matrix([[2,5],[4,1]]);
A3:=Matrix([[2,5,0],[4,1,0], [0,0,1]]);
#A:=Matrix([[4,2],[2,1]]);
#A3:=Matrix([[4,2,0],[2,1,0], [0,0,1]]);
N:=30:point2:=[]:line2:=[]:
for k from 0 to N-1 do
x0:=Vector([sin(2*Pi*k/N),cos(2*Pi*k/N)]);
x1:=A.x0;
p0:=convert(x0,list);
p1:=convert(x1,list);
point2:=[op(point2),disk(p0,0.1,color=blue, style=patchnogrid)];
point2:=[op(point2),disk(p1,0.1,color=red, style=patchnogrid)];
line2:=[op(line2),line(p0,p1)];
end do:
d:=6: display(point2,line2,view=[-d..d,-d..d]);
<Text-field style="Heading 1" layout="Heading 1">3d </Text-field>
N:=30:point2:=[]:line2:=[]:
for k from 0 to N-1 do
x0:=Vector([sin(2*Pi*k/N),cos(2*Pi*k/N),0]);
x1:=A3.x0+Vector([0,0,3]);
p0:=convert(x0,list);
p1:=convert(x1,list);
point2:=[op(point2),sphere(p0,0.15,color=blue, style=patchnogrid)];
point2:=[op(point2),sphere(p1,0.15,color=red, style=patchnogrid)];
line2:=[op(line2),line(p0,p1)];
end do:
d:=6: display(point2,line2,view=[-d..d, -d..d, 0..3], scaling=constrained);
d:=6: display(point2,line2,view=[-d..d, -d..d, 0..3], scaling=constrained, axes=box);