Овал Кассини

uses GraphABC,crt;
Var i,x0,y0,x1,y1,k,b,c: integer;
    x,y: real;
Begin
  setWindowSize(700,500);
   x0:=320; y0:=200; k:=30;
   line(x0,20,x0,380);   line(20,y0,620,y0);    {оси координат}
   line(x0,20,x0-10,30);  line(x0,20,x0+10,30);
   line(630,y0,620,y0-10);  line(630,y0,620,y0+10);
   for i:=-10 to 10 do
       begin
       line(x0+(i*k),y0-10,x0+(i*k),y0+10);
       textout(x0+(i*k),y0+10,inttostr(i));
       end;
   for i:=-5 to 5 do
       begin
       line(x0-10,y0+(i*k),x0+10,y0+(i*k));
       if i<>0 then textout(x0+10,y0+(i*k),inttostr(-i));
       end;
   setPenColor(clRed);
   write('Введите значение b: ');
   read(b);
   write('Введите значение c: ');
   read(c);
   x:=-10;  { начальное значение аргумента }
   While x<=10 do
   begin
     y:=sqr(sqr(exp(ln(b)*4)+4*c*c*x*x))-x*x-c*c); { функция }
     x1:=trunc(x0+x*k);  y1:=trunc(y0-y*k); { экранные координаты }
     circle(x1,y1,2);
     circle(x1,-y1,2);
     x:=x+0.001;
   end;
End.

x:=-10; 
 While x<=10 do
   begin
    z:=sqrt(b*b*b*b+4*c*c*x*x)-x*x-c*c;//выражение под первым корнем
    if z>=0 then
     begin
       y:=sqrt(z);
       setpixel(x0+round(x*k),y0-round(y*k),clRed);
       setpixel(x0+round(x*k),y0+round(y*k),clRed);
     end;
    x:=x+0.001;
   end;