Введение в Maple и рассмотрение задач теоретической механики. Эксаревская М.Е - 24 стр.

    кривошип := polygon([point1, point2, point3, point4],
                        linestyle=1, thickness=2);
    pointA := [OA*cos(phi), OA*sin(phi)];
    lineOA := line([0, 0], pointA, color=red, linestyle=3);
    psi := arcsin(OA*sin(phi)/AB);
    pointB := [pointA[1] + AB*cos(psi), 0];
    pointC := [pointB[1] - (AB-AC)*cos(psi),
               pointB[2] + (AB-AC)*sin(psi)];
    lineAB := line(pointA, pointB, color=red, linestyle=3);
    point1 := [pointA[1] + h*sin(psi), pointA[2] + h*cos(psi)];
    point2 := [pointB[1] + h*sin(psi), pointB[2] + h*cos(psi)];
    point3 := [pointB[1] - h*sin(psi), pointB[2] - h*cos(psi)];
    point4 := [pointA[1] - h*sin(psi), pointA[2] - h*cos(psi)];
    шатун := polygon([point1, point2, point3, point4],
                     linestyle=1, thickness=2);
    шарнирO := disk([0, 0], 3.5*h, color=green);
    шарнирA := disk(pointA, 3.5*h, color=green);
    шарнирB := disk(pointB, 3.5*h, color=green);
    h := 5*h;
    point1 := [pointB[1] + 2*h, pointB[2] + h];
    point2 := [pointB[1] + 2*h, pointB[2] - h];
    point3 := [pointB[1] - 2*h, pointB[2] - h];
    point4 := [pointB[1] - 2*h, pointB[2] + h];
    ползун := polygon([point1, point2, point3, point4], linestyle=1,
                      color=magenta, thickness=1);
    h := 3*h;
    textA := textplot([pointA[1] + h, pointA[2] + h, 'A']);
    textB := textplot([pointB[1] + h, pointB[2] + h, 'B']);
    textC := textplot([pointC[1] + h, pointC[2] + h, 'C']);
    display(lineOA, шарнирO, lineAB, шарнирA, ползун, шарнирB,
            textA, textB, textC,
            polygon([[0, 0], [-h, -h], [h, -h]],
                    linestyle=1, thickness=2),
            disk(pointC, h/5, color=green),
            scaling=constrained);
  end proc:

> draw := proc(phi)
    local
      OA, AB, AC, h, psi,
      point1, point2, point3, point4, pointA, pointB, pointC,
      lineOA, lineAB, кривошип, шатун, ползун,
      шарнирO, шарнирA, шарнирB, textA, textB, textC;
    OA := 35; AB := 70; AC := 45; h := .35;
    point1 := [h*cos(phi+Pi/2), h*sin(phi+Pi/2)];
    point2 := [h*cos(phi+Pi/2) + OA*cos(phi),
               h*sin(phi+Pi/2) + OA*sin(phi)];
    point3 := [point2[1] + 2*h*cos(phi+3*Pi/2),
               point2[2] + 2*h*sin(phi+3*Pi/2)];
    point4 := [h*cos(phi+3*Pi/2), h*sin(phi+3*Pi/2)];
    кривошип := polygon([point1, point2, point3, point4],
                        linestyle=1, thickness=2);
    pointA := [OA*cos(phi), OA*sin(phi)];
    lineOA := line([0, 0], pointA, color=red, linestyle=3);
    psi := arcsin(OA*sin(phi)/AB);
    pointB := [pointA[1] + AB*cos(psi), 0];
                                 24