Вращающийся куб

uses GraphABC;
 
const
  n = 8;
  R = 1;
  W = 500;
  H = 500;
  Dist = 10;
  
type
  Point3D = class
    X, Y, Z: real;
    constructor(px, py, pz: real);
    begin
      X := px;
      Y := py;
      Z := pz;
    end;
  end;
 
var
  Points3D: array [0..n - 1] of Point3D;
  Points2D: array [0..n - 1] of Point;
 
{sx = xSize/2+x*dist/(z+dist),
 sy = ySize/2-y*dist/(z+dist).}
 
procedure Project(i: integer);
var
  sx, sy: Integer;
begin
  try
    Points2D[i].X := Trunc(W / 2 + Points3D[i].X * dist / (Points3D[i].Z + dist));
    Points2D[i].Y := Trunc(H / 2 - Points3D[i].Y * dist / (Points3D[i].Z + dist));
  except end;
end;
 
begin
  SetWindowIsFixedSize(true);
  SetWindowSize(W, H);
  LockDrawing();
  SetSmoothingOff();
  
  for var i := 0 to n - 1 do Points3D[i] := new Point3D(0, 0, 0);
  for var i := 0 to n - 1 do Points2D[i] := new Point(0, 0);
  
  while true do
  begin
    for var a := 0 to 360 do
    begin
      ClearWindow();
      for var angle := 0 to 3 do
        for var i := 0 to 3 do
        begin
          Points3D[i].X := R * Cos(DegToRad(angle * 90 + a));
          Points3D[i].Y := R * Sin(DegToRad(angle * 90 + a));
          Points3D[i].Z := R;
          Project(i);
        end;
      
      for var angle := 0 to 3 do
        for var i := 4 to n-1 do
        begin
          Points3D[i].X := R * Cos(DegToRad(angle * 90 + a));
          Points3D[i].Y := R * Sin(DegToRad(angle * 90 + a));
          Points3D[i].Z := -R;
          Project(i);
        end;
      
      for var i := 0 to n - 1 do
        for var j := 0 to n - 1 do
          if i <> j then
            Line(Points2D[i].X, Points2D[i].Y, Points2D[j].X, Points2D[j].Y);
      Redraw();
      Sleep(100);
    end;
  end;
end.

  R = 9;

uses GraphABC; 
 
const
  L = 100;
  distance = 400;
 
type
  pix3D = record
    X,Y,Z : Real;
    
    constructor; begin (X, Y, Z) := (0, 0, 0); end;
    constructor (X, Y, Z : Real); begin (Self.X, Self.Y, Self.Z) := (X, Y, Z); end;
    
    function rotationX(angle : Real) := new pix3D(X, Y*cos(angle)+Z*sin(angle), -Y*sin(angle)+z*cos(angle));
    function rotationY(angle : Real) := new pix3D(X*cos(angle)-z*sin(angle), Y, X*sin(angle)+Z*cos(angle));
    function rotationZ(angle : Real) := new pix3D(X*cos(angle)+Y*sin(angle), -X*sin(angle)+Y*cos(angle), Z);
    
    function winX := Window.Center.X + Round(X * distance / (Z + distance));
    function winY := Window.Center.Y + Round(Y * distance / (Z + distance));
  end;
 
var pix : array of pix3D = (
  new pix3D(+L,+L,+L), new pix3D(+L,+L,-L), new pix3D(+L,-L,+L), new pix3D(+L,-L,-L),
  new pix3D(-L,+L,+L), new pix3D(-L,+L,-L), new pix3D(-L,-L,+L), new pix3D(-L,-L,-L)
);
 
procedure Li(a, b : pix3D) := Line(a.winX, a.winY, b.winX, b.winY);
 
procedure DrawCube(angleX, angleY, angleZ : Real);
begin
  var r : array of pix3D;
  SetLength(r, 8);
  for var i := 0 to 7 do r[i] := pix[i].rotationX(angleX).rotationY(angleY).rotationZ(angleZ);
  Li(r[0],r[1]); Li(r[0],r[2]); Li(r[1],r[3]); Li(r[2],r[3]); Li(r[0],r[4]); Li(r[1],r[5]);
  Li(r[4],r[5]); Li(r[2],r[6]); Li(r[4],r[6]); Li(r[3],r[7]); Li(r[5],r[7]); Li(r[6],r[7]);
end;
 
begin
  Randomize;
  var (rotX, rotY, rotZ) := (0.0, 0.0, 0.0);
  LockDrawing;
  repeat
    Window.Clear;
    rotX += Random/40; rotY += Random/40; rotZ += Random/40;
    DrawCube(rotX, rotY, rotZ);
    Redraw;
    Sleep(10);
  until False;
end.

var rotx:=0.0;
var roty:=0.0;
var rotz:=0.0;

uses Graph3D;
 
begin
  var куб := Cube(0,0,1,2,Colors.Green);
  while True do
  begin
    куб.Rotate(v3D(0,0,1),1);
    Sleep(10);
  end;
end.