Найти координаты точек касания к окружности

uses Math;//модуль для вычисления арксинусов
procedure TForm1.Button1Click(Sender: TObject);
var x,y:double; //коодинаты точки реальные
    xe,ye:integer;//экранные
    xo,yo:double;//координаты центра окружности реальные
    xoe,yoe:integer;//экранные
    r:double; //радиус окружности
    x1,y1,x2,y2:double;//координаты точек касания реальные
    x1e,y1e,x2e,y2e:integer;//экранные
    xk1,yk1,xk2,yk2:integer;//экранные координаты концов касательных
    a,b,d,d1:double;//переменые для вычисления точек касания
    c:integer;//центр начала кординат
    k:integer;//положение точек касания относительно оси Х
    i:integer;//счетчик цикла
    m:double;//масштаб
begin
if not TryStrToFloat(Edit1.Text,x)or(abs(x)>10)then
 begin
  ShowMessage('Неверно введена координата X');
  Edit1.Clear;
  Edit1.SetFocus;
  exit;
 end;
if not TryStrToFloat(Edit2.Text,y)or(abs(y)>10)then
 begin
  ShowMessage('Неверно введена координата Y');
  Edit2.Clear;
  Edit2.SetFocus;
  exit;
 end;
if not TryStrToFloat(Edit3.Text,r)or(r<1)or(r>5)then
 begin
  ShowMessage('Неверно введен радиус кружности');
  Edit3.Clear;
  Edit3.SetFocus;
  exit;
 end;
if not TryStrToFloat(Edit4.Text,xo)or(abs(xo)>5)then
 begin
  ShowMessage('Неверно введена координата X центра окружности');
  Edit4.Clear;
  Edit4.SetFocus;
  exit;
 end;
if not TryStrToFloat(Edit5.Text,yo)or(abs(yo)>5) then
 begin
  ShowMessage('Неверно введена координата Y окружности');
  Edit5.Clear;
  Edit5.SetFocus;
  exit;
 end;
d:=sqrt(sqr(x-xo)+sqr(y-yo));//расстояние от точки до центра окружности
if d<=r+1 then
 begin
  ShowMessage('Нельзя построить касательные');
  exit;
 end;
//экранные координаты начала координат  и масштаб
Image1.Width:=Image1.Height;
c:=Image1.Height div 2;
m:=(c-10)/10;
// Вычисление точек касания
d1:=sqrt(d*d-r*r); //расстояние от точки до точки касания
a:=arcsin((xo-x)/d);
b:=arcsin(r/d);
if y>yo then k:=-1 else k:=1;//точка выше или ниже оси
//координаты 1 точки касания
//реальные
x1:=x+d1*sin(a-b);
y1:=y+k*d1*cos(a-b);
//экранные
x1e:=c+round(x1*m);
y1e:=c-round(y1*m);
//координаты 2 точки касания
//реальные
x2:=x+d1*sin(a+b);
y2:=y+k*d1*cos(a+b);
//экранные
x2e:=c+round(x2*m);
y2e:=c-round(y2*m);
//координаты точки экранные
xe:=c+round(x*m);
ye:=c-round(y*m);
//координаты центра окружности экранные
xoe:=c+round(xo*m);
yoe:=c-round(yo*m);
//рисуем
 with Image1.Canvas do
 begin
 //коодинатную сетку
  moveto(c-round(10*m),c);
  lineto(c+round(10*m),c);
  moveto(c,c-round(10*m));
  lineto(c,c+round(10*m));
  pen.Style:=psDot;
  for i:=1 to 10 do
   begin
    moveto(c+round(i*m),c-round(10*m));
    lineto(c+round(i*m),c+round(10*m));
    moveto(c-round(i*m),c-round(10*m));
    lineto(c-round(i*m),c+round(10*m));
    moveto(c-round(10*m),c-round(i*m));
    lineto(c+round(10*m),c-round(i*m));
    moveto(c-round(10*m),c+round(i*m));
    lineto(c+round(10*m),c+round(i*m));
    moveto(c-round(10*m),c+round(10*m));
    textout(c+round(i*m)-10,c+5,inttostr(i));
    textout(c-round(i*m)+5,c+5,inttostr(-i));
    textout(c-20,c-round(i*m)+5,inttostr(i));
    textout(c-25,c+round(i*m),inttostr(-i));
   end;
  pen.Width:=2;
  pen.Color:=clRed;
  brush.Style:=bsClear;
  ellipse(xoe-round(r*m),yoe-round(r*m),xoe+round(r*m),yoe+round(r*m));
  pen.Width:=1;
  pen.Color:=clBlack;
  brush.Color:=clBlack;
  ellipse(xe-2,ye-2,xe+2,ye+2);
  pen.Color:=clBlue;
  moveto(xe,ye);
  lineto(x1e,y1e);
  moveto(xe,ye);
  lineto(x2e,y2e);
  end;
 end;

r-вводится
.......................
d:=sqrt(sqr(x-xo)+sqr(y-yo));//расстояние от точки до центра окружности
....................................................
// Вычисление точек касания
d1:=sqrt(d*d-r*r); //расстояние от точки до точки касания
a:=arcsin((xo-x)/d);

r:=sqrt(5);
d:=sqrt(x*x+y*y);//расстояние от точки до центра окружности
................................................
// Вычисление точек касания
d1:=sqrt(d*d-r*r); //расстояние от точки до точки касания
a:=arcsin(-x/d);

 a:=arcsin((xo-x)/d);