GraphABC + процедуры и функции (PascalABC.NET)

[Vova_P]

Здравствуйте! Хотел бы попросить помощи в доработке программы. Основа есть но выдает ошибку времени выполнения. Заранее спасибо.

Вот сама программа.Необходимо построить два графика на одном экране учитывая область определения и разрывы функций.

program sub_8;

uses
  GraphABC;

var
  a, b: real;
  n: integer;

function F(x: real): real;
begin
  F := sin(x) / cos(x) + 1 / sqrt(x * x - 9)
end;

function G(x: real): real;
begin
  G := exp(sqrt(x)) + Ln(x)
end;

procedure MaxMin(a, b, h: real; var yMin, yMax: real);
var
  x, y, a1, b1, b2, b3: real;
begin
  if a < -3 then
  begin
    yMin := F(a1);yMax := yMin;x := a1;
    while x < b1 + h / 2 do
    begin
      y := F(x);
      if y < yMin then
        yMin := y;
      if y > yMax then
        yMax := y;
      x := x + h;
    end;
  end;
  if b2 > 3 then
  begin
    while x < b2 + h / 2 do
    begin
      y := F(x);
      if y < yMin then
        yMin := y;
      if y > yMax then
        yMax := y;
      x := x + h;
    end;
  end;
  if b3 > 0 then
  begin
    while x < b3 + h / 2 do
      y := G(x);
    if y < yMin then
      yMin := y;
    if y > yMax then
      yMax := y;
    x := x + h;
  end;
  
end;


procedure GrFunc(a, b: real; n: integer);
var
  yMin, yMax: real;
  x, y, h, Mx, My: real;
  xe, ye, xep, yep, ox, oy,
  xemin, yemin, xemax, yemax, i, xemin1, yemin1, xep2, xe2, yep2, ye2: integer;
  function prx(x: real): integer;
  begin
    prx := round(xemin + (x - a) * Mx)
  end;
  
  function pry(y: real): integer;
  begin
    pry := round(yemin + (yMax - y) * My)
  end;



begin
  h := (b - a) / n; MaxMin(a, b, h, yMin, yMax);
  xemin := 20; xemax := WindowWidth - 20;
  yemin := 20; yemax := WindowHeight - 20;
  Mx := (xemax - xemin) / (b - a);
  My := (yemax - yemin) / (ymax - ymin);
  x := a; y := f(a);
  xep := prx(x); yep := pry(y);
  yep2 := pry(G(x));xep2 := prx(x);
  pen.Color := clGreen;
  for i := 1 to n do
  begin
    x := x + h; y := f(x);
    xe := prx(x);ye := pry(y);
    line(xep, yep, xe, ye);
    xep := xe; yep := ye;
    line(xep2, yep2, xe2, ye2);
    xep2 := xe2;yep2 := ye2;
  end;
  
      {выводятся оси}
  if yMin * yMax < 0 then
  begin
    oy := pry(0);pen.Color := clBlack;
    Line(xemin, oy, xemax, oy);
    TextOut(xemax - 15, oy + 10, 'Ox');
  end;
  if a * b < 0 then
  begin
    ox := prx(0);
    Line(ox, yemin, ox, yemax);
    TextOut(ox + 10, yemin + 15, 'Oy');
  end;
  {================================}
  xemin1 := xemin; yemin1 := yemin;
  pen.Color := ClGray; pen.Style := psDot;
  while xemin <= (xemax + round(ox / 15)) do
  begin
    line(xemin, yemin1, xemin, yemax);
    xemin := xemin + round(ox / 15);
  end;
  
  while yemin <= (yemax + round(oy / 15)) do
  begin
    line(xemin1, yemin, xemax, yemin);
    yemin := yemin + round(oy / 15);
  end;
  {=================================}
  var q, w, e, r: string;
  str(a, q);
  textout(xemin1, oy - 10, q );
  str(b, w);
  textout(xemax, oy - 10, w);
  yMax := round(yMax);
  str(yMax, e);
  textout(ox - 15, yemin1 + 5, e);
  yMin := round(yMin);
  str(yMin, r);
  textout(ox - 15, yemax - 10, r);
end;

begin
  n := 1000;
  write('Введите концы отрезка а..b');
  readln(a, b);
  GrFunc(a, b, n)
end.

[Vova_P]

Прошу не пугаться размеров,там большое ОДЗ просто