Здравствуйте уважаемые эксперты, подскажите пожалуйста ошибку в программе, программа должна выводить график, а она выдает ошитбку при компиляции что невозможно обнаружить модуль graph

Приложение:
uses crt,graph;
function F1(x:real):real;
begin
F1:=cos(x)/sin(x);{левая часть}
end;
function F2(x:real):real;
begin
F2:=1/x-x/2;{правая часть}
end;
function F(x:real):real;
begin
F:=F1(x)-F2(x);{исследуемая функция}
end;
{процедура уточнения корня на интервале мотодом половинного деления - бисекции}
procedure Bisec(a,b,e:real;var x:real);
var c:real;
begin
if F(a)*F(b)<0 then
begin
repeat
c:=(a+b)/2;
if F(a)*F(c)<=0 then b:=c
else a:=c;
until b-a<e;
x:=(a+b)/2;
end
else x:=a-1;
end;
var a,b,x,e:real;{границы интервала, значение корня, точность}
k:array[1..20] of real;{массив корней}
n,i:integer;{количество корней, счетчик циклов}
gd,gm,x0,y0,xi,y1,y2:integer;{переменные для работы в граф. режиме}
mx,my:real;{массшатабы для перевода реальных координат в экранные}
s:string[10];{строка для перевода чисел в текст для вывода на экран в граф. режиме}
begin
clrscr;
writeln('F(x)=ctg(x)-1/x+x/2');
e:=0.00001;n:=0;
a:=0.5;b:=pi-e;{период функции ctg=pi}
while b<=20 do
begin
Bisec(a,b,e,x);{вычисляем корень на отрезке}
n:=n+1; {считаем}
k[n]:=x;{пишем в массив}
a:=a+pi;{переходим в новый отрезок}
b:=b+pi;
end;
writeln('Корни уравнения не интервале 0.5 20:');
for i:=1 to n do
writeln('x',i,'=',k[i]:9:6);
write('Press Enter...');
readln;
gd:=0;
initgraph(gd,gm,'');{переходим в графический режим}
x0:=30;y0:=getmaxY div 2;{начало координат}
mx:=(getmaxX-15-x0)/20;
my:=(getmaxY-20)/(2*abs(F2(20)));{массштабы = размеры экрана по осям/максимальные значения}
{ось Х}
line(x0,y0,getmaxX,y0);
outtextXY(getmaxX-10,y0-15,'X');
for i:=1 to 21 do
begin
line(x0+round((i-1)*mx),y0-3,x0+round((i-1)*mx),y0+3);
str(i-1,s);
outtextXY(x0+round((i-1)*mx),y0+10,s);
end;
{ось У}
line(x0,5,x0,getmaxY-5);
outtextXY(x0-20,5,'Y');
for i:=1 to round(abs(F2(20))) do
begin
line(x0+3,y0-round(i*my),x0-3,y0-round(i*my));
str(i,s);
outtextXY(x0-10,y0-round(i*my),s);
line(x0+3,y0+round(i*my),x0-3,y0+round(i*my));
str(i,s);
outtextXY(x0-20,y0+round(i*my),'-'+s);
end;
{подписи значений корней}
settextstyle(0,1,1);
setcolor(11);
for i:=1 to n do
begin
insert('.',s,length(s)-5);
str(k[i]:8:6,s);
outtextXY(x0+round(k[i]*mx),y0-80,s);
setlinestyle(0,1,3);
circle(x0+round(k[i]*mx),y0,1);
end;
{построение графика}
x:=0.5;
while x<=20 do
begin
xi:=x0+round(x*mx);
if sin(x)<>0 then
begin
y1:=y0-round(F1(x)*my);
putpixel(xi,y1,12);
end;
y2:=y0-round(F2(x)*my);
putpixel(xi,y2,10);
x:=x+0.01;
end;
readln
end.