Программа 17. Задача 26

>	restart;read "C:\\Podpr1.m"; with(plots):

Угловая скорость

>	omega[1]:=4: epsilon[1]:=0:

Координаты

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

>	x[2]:=4: y[2]:=3:

>	x[3]:=4: y[3]:=-1:

>	x[4]:=8: y[4]:=-1:

Уравнения 3х угловых скоростей

>	eq1:=sum(omega[i]*(x[i]-x[i+1]),i=1..3):

>	eq2:=sum(omega[i]*(y[i]-y[i+1]),i=1..3):

>	R:=solve({eq1,eq2},{omega[2],omega[3]}):assign(R):

>	omega2=omega[2],omega3=omega[3];

Уравнения 3х угловых ускорений

>	eq3:=sum(epsilon[i]*(x[i]-x[i+1]),i=1..3)-sum(omega[i]^2*(y[i]-y[i+1]),i=1..3):

>	eq4:=sum(epsilon[i]*(y[i]-y[i+1]),i=1..3)+sum(omega[i]^2*(x[i]-x[i+1]),i=1..3):

>	solve({eq3,eq4},{epsilon[2],epsilon[3]});

>	display(ЛиниЯ(1,2),ЛиниЯ(2,3),ЛиниЯ(3,4),Опора(1,0.1),Опора(4,0.1),axes=NONE,scaling=constrained, title=`Четырехзвенник`);

>