Исследовать область определения функции и построить график - Turbo Pascal

uses graph;
function F(x:real):real;
begin
F:=x*x*x+2*x*x+x;
end;
 
var a,b,mx,my,x:real;
    x0,y0,i:integer;
    s:string;
begin
x0:=0;
initgraph(x0,y0,'');{инициализация графического режима}
a:=-5;
b:=5;{интервал по Х}
x0:=getmaxX div 2;{начало координат по Х}
mx:=(x0-30)/b;{масштаб по Х}
y0:=getmaxY div 2;{начало координат по Y}
my:=(y0-40)/F(b);
line(0,y0,getmaxX,y0);{оси}
line(x0,0,x0,getmaxY);
for i:=1 to 5 do{максимальное количество засечек в одну сторону}
 begin
  line(x0-3,y0-round(i*my*50),x0+3,y0-round(i*my*50));{засечки на оси У}
  line(x0-3,y0+round(i*my*50),x0+3,y0+round(i*my*50));
  {подпись оси У}
  str(i*50,s);
  outtextXY(x0-35,y0-round(i*my*50),s);{соответственно засечкам}
  outtextXY(x0-40,y0+round(i*my*50),'-'+s);
  line(x0+round(i*mx),y0-3,x0+round(i*mx),y0+3); {засечки на оси Х}
  line(x0-round(i*mx),y0-3,x0-round(i*mx),y0+3);
  {подпись оси Х}
  str(i,s);
  outtextXY(x0+round(i*mx)+5,y0+10,s);
  outtextXY(x0-round(i*mx)+5,y0+10,'-'+s);
 end;
{центр}
outtextXY(x0-15,y0+10,'0');
{подписи концов осей}
outtextXY(getmaxX-10,y0-10,'X');
outtextXY(x0+5,10, 'Y');
{график}
x:=a;
setcolor(12);
while x<=b do
 begin
  if x=a then moveto(x0+round(x*mx),y0-round(F(x)*my))
  else lineto(x0+round(x*mx),y0-round(F(x)*my));
  x:=x+0.001;
 end;
{название}
outtextXY(100,10,'y=x^3+2x^2+x');
outtextXY(100,30,'int.[-5;5]');
readln
end.