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

13.2. Задача 8

>

restart:

>	read "ris.m";

               Подписи шарниров на рисунке

>	nam:=[O,A,B,C,D]:            Pазмеры ( cм) OA:=33: DB:=63:

>	AB:=30: BC:=23:

                Координаты опоры O

>	x[1]:=0:y[1]:=0:

               Координаты опоры D

>	x[5]:=AB*cos(pi/6)+DB:

>	y[5]:= -OA+AB*sin(pi/6):   x[4]:=(AB+BC)*cos(pi/6):

 Количество кадров K

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

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

>	for i from 0 to K do

       Закон движения  АO

>	Амплитуда:=pi/60:

>	phi:=3*pi/2+sin(2*pi*i/K)*Амплитуда:            Координаты шарнира А                         N1

>	x[2]:= OA*cos(phi): y[2]:= OA*sin(phi): koord(5,2,3,DB,AB): y[4]:=y[3]-sqrt(BC^2-(x[3]-x[4])^2): P[i]:=display(seq(Cir(i,1),i=2..4), box(4,2,5), seq(Lin(i,i+1),i=1..3), Lin(3,5),seq(TEXT([x[j]+3,y[j]+4],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,seq(Опора(i,1),i=[1,5]),Стенка(x[4]-3,-38,15,2));

>