Построение ограниченых областей плоскости

uses graph,wincrt;
function F1(x:real):real;
begin
F1:=sqrt(2*x*x+7*x+13);
end;
function F2(x:real):real;
begin
F2:=-9*x*x+15;
end;
function F3(x:real):real;
begin
F3:=5*x+10;
end;
var max,min,mx,my,h:real;
    a,b,x,y,x1,y1,x2,y2,x3,y3,x4,y4:real;
    x0,y0,i,n:integer;
    s:string;
begin
initgraph(x0,y0,'');
a:=-1.3;
b:=1.3;{интервал по Х}
min:=F2(a);
max:=min;
x:=a;
while x<=b do
 begin
  if F1(x)<min then min:=F1(x);
  if F2(x)<min then min:=F2(x);
  if F1(x)>max then max:=F1(x);
  if F2(x)>max then max:=F2(x);
  x:=x+0.01
 end;
if abs(min)>max then max:=abs(min);
x0:=getmaxX div 2;
mx:=(x0-30)/b;{масштаб по Х}
y0:=getmaxY-70;
my:=(y0-30)/max;
line(0,y0,getmaxX,y0);{оси}
line(x0,0,X0,getmaxY);
{settextstyle(0,0,2);}
for i:=1 to trunc(max)+1  do{максимальное количество засечек в одну сторону}
 begin
  line(x0-3,y0-round(i*my),x0+3,y0-round(i*my));{засечки на оси У}
  line(x0-3,y0+round(i*my),x0+3,y0+round(i*my));
  {подпись оси У}
  settextstyle(0,0,1);
  str(i,s);
  outtextXY(x0-20,y0-round(i*my),s);{соответственно засечкам}
  outtextXY(x0-20,y0+round(i*my),'-'+s);
  line(x0+round(i*mx/10),y0-3,x0+round(i*mx/10),y0+3); {засечки на оси Х}
  line(x0-round(i*mx/10),y0-3,x0-round(i*mx/10),y0+3);
  {подпись оси Х}
  str(i/10:0:1,s);
  settextstyle(0,1,1);
  outtextXY(x0+round(i*mx/10),y0+10,s);
  outtextXY(x0-round(i*mx/10),y0+10,'-'+s);
 end;
settextstyle(0,0,1);
{центр}
outtextXY(x0+5,y0+10,'0');
{подписи концов осей}
outtextXY(getmaxX-10,y0-10,'X');
outtextXY(x0+5,10, 'Y');
{график}
x:=a;
while x<=b do
 begin
  putpixel(x0+round(x*mx),y0-round(F1(x)*my),12);
  putpixel(x0+round(x*mx),y0-round(F2(x)*my),9);
  putpixel(x0+round(x*mx),y0-round(F3(x)*my),10);
  if x<0 then
   begin
    if abs(F1(x)-F2(x))<0.1 then
     begin
      x1:=x;
      y1:=F1(x)
     end
    else if abs(F2(x)-F3(x))<0.1 then
     begin
      x2:=x;
      y2:=F2(x)
     end
   end
  else if x>0 then
   begin
    if abs(F2(x)-F3(x))<0.1 then
     begin
      x3:=x;
      y3:=F2(x)
     end
    else if abs(F1(x)-F2(x))<0.1 then
     begin
      x4:=x;
      y4:=F1(x)
     end
   end;
  x:=x+0.001;
 end;
{название}
setcolor(12);
outtextXY(130,10,'y=sqrt(2x^2+7x+13)');
setcolor(9);
outtextXY(130,30,'y=-9x^2+15');
setcolor(10);
outtextXY(130,50,'y=5x+10');
h:=10/my;
n:=round((y3-y1)/h);
for i:=n downto 1 do
 begin
  x:=x1;
  while x<=x4 do
   begin
    x:=x+0.001;
    y:=y1+h*(n-i);
    if (x>x1)and(x<x4)and(y>f1(x))and(y<f2(x))and(y<f3(x))
    then putpixel(x0+round(x*mx),y0-round(y*my),13);
   end;
 end;
readkey
end.