Кривая Минковского

program krivaya_minkovskogo;
 
uses
  GraphABC;
 
var
  sx, sy, hu, hv, fi: real;
 
procedure Init;
begin
  hu := 3;
  hv := window.height;
  sx := window.width;
  sy := window.height;
  fi := Pi / 3;
end;
 
procedure Outline(x1, y1, x2, y2: real; c: color);
var
  u1, v1, u2, v2: word;
begin
  u1 := round(x1 * sx + hu);
  v1 := round(-y1 * sy + hv);
  u2 := round(-y2 * sy + hv);
  pen.Color := c;
  line(u1, v1, u2, v2);
  
end;
 
procedure Draw(x1, y1, x2, y2: real; n: word);
var
  x3, x4, x5, y3, y4, y5, dx, dy: real;
begin
  if (n > 0) then
  begin
    dx := (x2 - x1) / 3;dy := (y2 - y1) / 3;
    x3 := x1 + dx;y3 := y1 + dy;
    x5 := x3 + dx;y5 := y3 + dy;
    x4 := (x5 - x3) * cos(fi) + (y5 - y3) * sin(fi) + x3;
    Draw(x1, y1, x3, y3, n - 1);
    Draw(x3, y3, x4, y4, n - 1);
    Draw(x4, y4, x5, y5, n - 1);
    Draw(x5, y5, x2, y2, n - 1);
  end
  else
    Outline(x1, y1, x2, y2, clred);
end;
 
 
begin
  Init;
  Draw(0, 0, 1, 0, 5);
  Draw(0, 0, 1 / 2, sqrt(3) / 2, 5);
  Draw(1 / 2, sqrt(3 / 2), 1, 0, 5);
end.