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

uses graph;
const xn=-4;
      xk=4;
function F(x:real):real;
begin
F:=3*x*x;
end;
var mx,my,dx,y1:real;
    x:real;
    x0,y0,i:integer;
    s:string;
begin
x0:=0;
initgraph(x0,y0,'');
x0:=getmaxX div 2;
mx:=(x0-30)/xk;{масштаб по Х}
y0:=getmaxY-60;
my:=(y0-20)/F(xk);{по У}
line(0,y0,getmaxX,y0);{оси}
line(x0,0,X0,getmaxY);
for i:=1 to round(F(xk)/10)+1 do{максимальное количество засечек в одну сторону}
 begin
  line(x0-3,y0-round(i*my*10),x0+3,y0-round(i*my*10));{засечки на оси У}
  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*10,s);
  {подпись оси У}
  outtextXY(x0-25,y0-round(i*my*10),s);{соответственно засечкам}
  {подпись оси Х}
  str(i,s);
  outtextXY(x0+round(i*mx),y0+10,s);
  outtextXY(x0-round(i*mx),y0+10,'-'+s);
  end;
{центр}
outtextXY(x0+5,y0+10,'0');
{подписи концов осей}
outtextXY(getmaxX-10,y0-10,'X');
outtextXY(x0+5,10, 'Y');
{график}
x:=xn;
dx:=0.001;
while x<=xk do
 begin
  putpixel(x0+round(x*mx),y0-round(F(x)*my),10);
  x:=x+dx;
 end;
{название}
outtextXY(50,20,'Grafik functii y=3x^2');
outtextXY(70,40,'na intervale [-4;4]');
outtextXY(90,60,'ODZ (-oo;+oo)');
readln
end.

initgraph(x0,y0,'');