Найти область определения функции D(f), ее нули, экстремумы и множества - Pascal ABC

uses graphABC;
function F(x:real):real;
begin
F:=-x*x*x*x+4*x*x;
end;
var xn,xk,x,mx,my,dx:real;
    x0,y0,i:integer;
begin
x0:=windowwidth div 2;
y0:=windowheight div 2;
xn:=-3;xk:=3;{интервал по Х}
mx:=(x0-30)/xk;{масштаб по Х}
my:=(y0-20)/-F(xn);{по У}
for i:=1 to 6 do
 begin
  line(x0+round(i*mx/2),y0-3,x0+round(i*mx/2),y0+3);
  line(x0-round(i*mx/2),y0-3,x0-round(i*mx/2),y0+3);
  line(x0-3,y0-round(i*my*10),x0+3,y0-round(i*my*10));
  line(x0-3,y0+round(i*my*10),x0+3,y0+round(i*my*10));
  textout(x0+round(i*mx/2),y0+10,floattostr(i/2));
  textout(x0-round(i*mx/2),y0+10,floattostr(-i/2));
  textout(x0-20,y0-round(i*my*10),inttostr(i*10));
  textout(x0-25,y0+round(i*my*10),inttostr(-i*10));
 end;
line(0,y0,windowwidth,y0);{оси}
line(windowwidth,y0,windowwidth-10,y0+5);
line(windowwidth,y0,windowwidth-10,y0-5);
line(x0,0,X0,windowheight);
line(x0,0,x0-5,10);
line(x0,0,x0+5,10);
textout(x0+5,y0+10,'0');
{подписи концов осей}
textout(windowwidth-10,y0-20,'X');
textout(x0+10,0, 'Y');
{график}
x:=xn;
setpencolor(clBlue);
setpenwidth(2);
dx:=0.1;
while x<=xk do
 begin
  if x=xn then moveto(x0+round(x*mx),y0-round(F(x)*my))
  else lineto(x0+round(x*mx),y0-round(F(x)*my));
  x:=x+dx; {наращиваем х}
 end;
setfontsize(10);
setfontcolor(clBlue);
textout(20,0,'y=-x^4+4x^2');
textout(20,20,'D(f)=R');
textout(20,40,'W(f)=[-oo;4)');
textout(20,60,'min x=0 y=0');
textout(20,80,'max x1=-sqrt(2) y=4 x2=sqrt(2) y=4');
textout(20,100,'корни x1=-2 x2=0 x3=2');
end.