Кирсанов М.Н. 
Решебник.Теоретическая механика.
 М.:ФИЗМАТЛИТ, 2002.
Иллюстрации и программы к книге

13.8. Задача 6

>    restart:

>    read "ris.m":

>    nam:=["O","A","O1","A1","B"]:
a:=25:d:=1:h:=70:b:=35:

                    Координаты

>      x[1]:=0:y[1]:=0:
  x[3]:=0:y[3]:=h:
  y[5]:=b:y[6]:=b:y[7]:=b:
    
       Количество кадров K

>    K:=12: with(plots):with(plottools):

           Создаем все кадры

>    for i from 0 to K do

>     t:=i/K:

>      phi:=sin(2*pi*i/K)*pi/6+pi/6:
  x[2]:=a*cos(phi):y[2]:=a*sin(phi):
  x[4]:=x[2]:      y[4]:=y[2]+h:
  x[5]:=x[2]:      x[6]:=x[5]+h:
  x[7]:=x[2]+2:
  P[i]:=display(Cir(2,1),Cir(4,1),
Lin(1,2),Lin(3,4),Lin(7,6),Line(2,4,8),box(5,2,6),
seq(TEXT([x[j]-4,y[j]],nam[j]),j=1..5)
):
od:

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

>    PP:=display(seq(P[i],i=0..K),insequence=true,
                            thickness=2,
                            scaling=constrained,
                            axes=none):

Warning, the name changecoords has been redefined

Warning, the name arrow has been redefined

>    display(PP,Опора(1,1),Опора(3,1),Поверхность(30,b-1,11,2),
Поверхность(30,b+1,11,-2));

[Maple Plot]

>