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

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

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

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

    Определение координат центра тяжести плоской фигуры

  restart:
Число точек на дуге и число угловых точек
N:=20:K:=4:
pi:=evalf(Pi):  n:=K+N:

     Матрица координат
   F:=i*pi/N: X:=1,0,4,4: Y:=0,2,3,0:
     Точки на дуге   
   Xcircl:=seq(3+cos(F),i=1..N):
   Ycircl:=seq(sin(F),i=1..N):
     Матрица координат   
   T:=matrix(2,n+1,[X,Xcircl,1,Y,Ycircl,0]):
   with(plots):polygonplot(T,scaling=constrained);

     Площадь
   Z:=T[2,i]*T[1,i+1]-T[1,i]*T[2,i+1]:
   A:=add(Z,i=1..n)/2;

     Координаты центра тяжести
      Xc:=add(Z*(T[1,i]+T[1,i+1]),i=1..n)/6/A;
  Yc:=add(Z*(T[2,i]+T[2,i+1]),i=1..n)/6/A;