Рисование параболы

uses
  graphABC;
 
var
  k, xq, yq, x0, y0, px, py: integer;
  
  a, x, y, b, x1,  x2, mx, my, dx: real;
 
procedure coord;//Координатна сітка
var
  n: integer;
begin
  xq := k;
  yq := k;
  MoveTo(WindowWidth div 2, windowHeight);//Вісь Y
  lineTO(WindowWidth div 2, 0);
  lineto((WindowWidth div 2) + 5, 5);
  line(WindowWidth div 2, 0, (WindowWidth div 2) - 5, 5);
  TextOut(WindowWidth div 2 + 10, 10, 'y');
  
  MoveTo(0, WindowHeight div 2);//Вісь Х
  lineTO(windowWidth, WindowHeight div 2); 
  lineto(windowWidth - 5, windowHeight div 2 + 5);
  line(windowWidth, WindowHeight div 2, windowWidth - 5, (windowHeight div 2) - 5);
  TextOut(windowWidth - 10, WindowHeight div 2 + 5, 'x');
  TextOut(WindowWidth div 2 + 5, WindowHeight div 2 - 16, '0');
  
  n := WindowHeight div yq;//відстань між штрихами
  while n < windowHeight do //штрихи на Y<0
  begin
    line(WindowWidth div 2 - 5, WindowHeight div 2 + n, windowwidth div 2 + 5, WindowHeight div 2 + n);
    n := n + WindowHeight div yq
  end;
  while n > 5 do //штрихи на Y>0
  begin
    line(WindowWidth div 2 - 5, WindowHeight div 2 - n, windowwidth div 2 + 5, WindowHeight div 2 - n);
    n := n - WindowHeight div yq
  end;
  
  n := Windowwidth div xq;//відстань між штрихами
  while n < windowwidth do // штрихи на Х<0
  begin
    line(Windowwidth div 2 + n, WindowHeight div 2 - 5, WindowWidth div 2 + n, windowHeight div 2 + 5 );
    n := n + Windowwidth div xq;
  end;
  
  while n > 5 do // штрихи на Х>0
  begin
    line(Windowwidth div 2 - n, WindowHeight div 2 - 5, WindowWidth div 2 - n, windowHeight div 2 + 5 );
    n := n - Windowwidth div xq;
  end;
end;
 
procedure Input;
begin
  write('(y²=2аx+b):');
  write('Введіть а:');
  readln(a);
  write(a, ' Введіть b:');
  readln(b);
  K := 10;
end;
 
procedure graph;
begin
  x0 := Windowwidth div 2;
  Y0 := WindowHeight div 2;
  x1 := -k / 2;
  x2 := k / 2;
  dx := 0.001;
  mx := Windowwidth / (x2 - x1);//довжина в пікселях по осі Х/реальна довжина
  my := WindowHeight / k;//довжина осі по У/(максимум прийнятих значень
  x := x1;
  if a < 0 then begin
    while x <= x2 do
    begin
      px := x0 + round(x * mx); //зміна по осі Х
      y := sqrt(abs(a * 2 * x + b)); //графік в реальних координатах
      py := y0 - round(y * my);//в экранних
      if x < -b / (2 * a) then  putpixel(px, py, clred);
      x := x + dx; //нарощуєм х
    end;
    x := x1;
    while x <= x2 do
    begin
      px := x0 + round(x * mx); //зміна по осі Х
      y := -sqrt(abs(a * 2 * x + b)); //графік в реальних координатах
      py := y0 - round(y * my);//в экранних
      if x < -b / (2 * a) then putpixel(px, py, clred);
      x := x + dx; //нарощуєм х
    end;
  end;
  
  if a > 0 then begin
    while x <= x2 do
    begin
      px := x0 + round(x * mx); //зміна по осі Х
      y := sqrt(abs(a * 2 * x + b)); //графік в реальних координатах
      py := y0 - round(y * my);//в экранних
      if x > -b / (2 * a) then  putpixel(px, py, clred);
      x := x + dx; //нарощуєм х
    end;
    x := x1;
    while x <= x2 do
    begin
      px := x0 + round(x * mx); //зміна по осі Х
      y := -sqrt(abs(a * 2 * x + b)); //графік в реальних координатах
      py := y0 - round(y * my);//в экранних
      if x > -b / (2 * a) then putpixel(px, py, clred);
      x := x + dx; //нарощуєм х
    end;
  end;
  TextOut(3, 3, 'Стрілки для масштабування, F12 для виходу')
end;
 
procedure Resize;//обробник зміни розміру вікна
begin
  ClearWindow;
  graph;
  coord;
  Redraw;
end;
 
procedure KeyDown(Key: integer);
begin
  case Key of
    VK_F12: halt;
    VK_Up:    k := k + 1;
    VK_Down:  if k > 1 then k := k - 1;
  end;
  ClearWindow;
  graph;
  coord;
  Redraw;
end;
 
begin
  input;
  ClearWindow;
  graph;
  coord;
  // Прив'язка обрабників до подій
  OnKeyDown := KeyDown;
  onResize := resize;
end.