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

13.8. Задача 5

>    restart:

>    read "ris.m":

>    a:=30:d:=1:h:=25:b:=22:

                    Координаты

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

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

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

>    for i from 0 to K do

>     t:=i/K:

>      phi:=sin(2*pi*i/K-1)*pi/7+pi/2:
  x[2]:=a*cos(phi):
  y[2]:=a*sin(phi):
  y[3]:=y[2]+d:y[4]:=y[3]:y[5]:=y[3]:y[6]:=y[3]+h:
  P[i]:=display(Cir(2,1),Cir(2,0.8),
  Lin(1,2),Line(3,4,8),Line(5,6,8)):
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),Стенка(-1,33,8,1),Стенка(1,33,8,-1));

[Maple Plot]

>   

>