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

uses crt,graph;
var gd,gm,x0,y0,i,xi,yi:integer;
    xn,xk,mx,my,x,y:real;
    s:string;
begin
gd:=0;
initgraph(gd,gm,'');
x0:=getmaxX div 2;
y0:=getmaxY-40;
outtextXY(40,20,'Grafik functii y=3+2/x+3/x^2');
outtextXY(40,30,'na intervale [-6,6]');
outtextXY(40,40,'oblast opredelenya x<>0');
{Ё*вҐаў*« Ја*дЁЄ*}
xn:=-6;
xk:=6;
{¬*бив*Ўл Ї® ®бп¬}
mx:=(x0*2-40)/(xk-xn);
my:=(y0-20)/20;
{®бм •}
line(0,y0,getmaxX,y0);
outtextXY(x0+5,y0+10,'0');
outtextXY(x0+round(xk*mx)+5,y0-15,'X');
for i:=1 to round(xk) do
 begin
  str(i,s);
  line(x0+round(i*mx),y0+3,x0+round(i*mx),y0-3);
  outtextXY(x0+round(i*mx)+5,y0+10,s);
  line(x0-round(i*mx),y0+3,x0-round(i*mx),y0-3);
  outtextXY(x0-round(i*mx)-5,y0+10,'-'+s);
 end;
{Ћбм “}
line(x0,y0,x0,10);
outtextXY(x0-10,10,'Y');
for i:=1 to 20 do
 begin
  str(i,s);
  line(x0+3,y0-round(i*my),x0-3,y0-round(i*my));
  outtextXY(x0-15,y0-round(i*my),s);
 end;
{Ја*дЁЄ}
x:=xn;
while x<=xk do
 begin
  xi:=x0+round(x*mx);
  if abs(x)>0.001 then
   begin
    y:=3+2/x+3/(x*x);
    yi:=y0-round(y*my);
    putpixel(xi,yi,12);
   end;
  x:=x+0.001;
 end;
readln
end.