Кирсанов М.Н. 

Решебник.Теоретическая механика.

 М.:ФИЗМАТЛИТ, 2002.

Иллюстрации и программы к книге

13.8. Задача 1

>

restart:

>	read "ris.m":

>	R:=3: a:=6:

                    Координаты

>	y[1]:=R:             Количество кадров 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+3*pi/2:   x[1]:=-R*phi:   x[2]:=x[1]+a*cos(phi):   y[2]:=y[1]+a*sin(phi):   F:=arrow([x[2],y[2]+4],evalm([0,-20]*0.2),0.1,0.6,0.2,color=red):   P[i]:=display(F,Cir(1,R),Cir(2,0.2),Cir(1,0.2),   cir4(1,R,phi),   Line(1,2,3)): 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,Поверхность(-18,0,8,0.5));

>