Построить график ( 1/x )

uses graphABC;
function f1(x:real):real;
begin
f1:=1/x
end;
function f2(x:real):real;
begin
f2:=sqrt(9-x*x)
end;
var c,i:integer;
    d,t1,t2,x1,y1,x2,y2,x3,y3,x4,y4,x,y,m:real;
begin
//определим размер экрана
setwindowsize(500,500);
//расположим его в центре экрана
centerwindow;
//определим центр графика
c:=windowwidth div 2;
m:=(c-20)/3;//его масштаб
circle(c,c,round(3*m));//рисуем окружность
line(0,c,2*c-10,c);//и оси
setpenwidth(1);
line(c,10,c,2*c);
//делаем разбивку на осях
for i:=1 to 3 do
 begin
  line(c+round(i*m),c-3,c+round(i*m),c+3);
  line(c-round(i*m),c-3,c-round(i*m),c+3);
  textout(c+round(i*m),c+10,inttostr(i));
  textout(c-round(i*m)+5,c+10,inttostr(-i));
  line(c-3,c-round(i*m),c+3,c-round(i*m));
  line(c-3,c+round(i*m),c+3,c+round(i*m));
  textout(c-20,c-round(i*m),inttostr(i));
  textout(c-25,c+round(i*m),inttostr(-i));
 end;
textout(c-10,c+10,'0');
line(c,0,c-5,10);
line(c,0,c+5,10);
line(c-5,10,c+5,10);
floodfill(c,5,clBlack);
textout(c+8,0,'Y');
line(2*c,c,2*c-10,c-5);
line(2*c,c,2*c-10,c+5);
line(2*c-10,c-5,2*c-10,c+5);
floodfill(2*c-5,c,clBlack);
textout(2*c-8,c-20,'X');
//найдем точки пересечения графика y=1/x и окружности x^2+y^2=9
//для чего решим биквадратное уравнение x^4-9x^2+1=0
d:=(81-4);
t1:=(9+sqrt(d))/2;
t2:=(9-sqrt(d))/2;
x1:=-sqrt(t1);y1:=1/x1;
x2:=-sqrt(t2);y2:=1/x2;
x3:=sqrt(t2);y3:=1/x3;
x4:=sqrt(t1);y4:=1/x4;
setpenwidth(3);
//построим левую ветку графика y=1/x
x:=x1;
moveto(c+round(x*m),c-round(f1(x)*m));
while x<=x2 do
 begin
  lineto(c+round(x*m),c-round(f1(x)*m));
  x:=x+0.01;
 end;
//построим правую ветку графика y=1/x
x:=x3;
moveto(c+round(x*m),c-round(f1(x)*m));
while x<=x4 do
 begin
  lineto(c+round(x*m),c-round(f1(x)*m));
  x:=x+0.01;
 end;
//заштрихуем области
setbrushstyle(bsFDiagonal);
floodfill(c+round(2*m),c-round(2*m),clBlack);
floodfill(c-round(0.5*m),c+round(0.5*m),clBlack);
floodfill(c+round(0.5*m),c+round(0.5*m),clBlack);
line(c-round(3*m),c,c+round(3*m),c);
//обведем внешний контур области
x:=-3;
moveto(c+round(x*m),c-round(f2(x)*m));
while x<=x1 do
 begin
  lineto(c+round(x*m),c+round(f2(x)*m));
  x:=x+0.001;
 end;
x:=x2;
moveto(c+round(x*m),c+round(f2(x)*m));
while x<3 do
 begin
  lineto(c+round(x*m),c+round(f2(x)*m));
   x:=x+0.001;
 end;
x:=x3;
moveto(c+round(x*m),c-round(f2(x)*m));
while x<=x4 do
 begin
  lineto(c+round(x*m),c-round(f2(x)*m));
  x:=x+0.001;
 end;
end.