Программа 2

>

restart:

>	with(plots):with(plottools):

>	K:=48: #Количество кадров

>	R1:=30: R2:=20: OA:=R1+R2:#Размеры

>	Водило:= PLOT3D(cuboid([-3,-3,0],[OA+4,3,2])):

>	Колесо1:=PLOT3D(cylinder([0,0,5],R1,18),cylinder([0,0,-2],1,16)):

>	Колесо2:=PLOT3D(cylinder([0,0,7],R2,6),cylinder([0,0,-2],1,16)):

>	for i from 0 to K do

>	Амплитуда:=0.3:

>	phi:=sin(6.28*i/K)*Амплитуда:

>	phi1:=phi*(R1+R2)/R2:

>	P1:=rotate(Колесо2,0,0,-phi1):

>	P2:=rotate(Водило,0,0,-phi):

>	x:= OA*cos(phi): y:= OA*sin(phi):

>	P[i]:=display(Колесо1,translate(P1,x,y,0),P2);

>

od:

     Изображение механизма в движении

>	display(seq(P[i],i=0..K),orientation=[55,101],insequence=true, scaling=constrained,axes=normal);

>