Заполнение окружности окружностями

procedure TForm1.Button1Click(Sender: TObject);
const
  N = 10000;
var
  mas,mas1 : array[1..2,1..N] of Integer;
  i, j, m: Integer;
  b : boolean;
begin
  R1:=StrToInt(Edit1.Text);
  R2:=StrToInt(Edit2.Text);
  k:=1;
 
  with PaintBox1.Canvas do
  Ellipse(0,0,2*R1,2*R1);
 
  for i:=1 to 20*Round(Sqr(R1)/sqr(R2)) do
  begin
    x0:=Random(2*R1-R2);
    y0:=Random(2*R1-R2);
 
    if Sqrt(sqr(x0 - R1)+sqr(y0 - R1))+R2<=R1 then
    begin
      mas[1,k]:=x0;
      mas[2,k]:=y0;
      k:=k+1;
    end;
  end;
 
  m:=0;
  with PaintBox1.Canvas do
  for i:=1 to k do
  begin
    b:=false;
    //Проверяем очередную окружность на непересечение с уже нарисованными
    for j:=1 to m do
    if Sqrt(sqr(mas1[1,j] - mas[1,i])+sqr(mas1[2,j] - mas[2,i]))<2*R2 then
    begin
      b:=True; //пересекается
      break;
    end;
 
    if Not b then
    //если не пересекается ни с одной ранее нарисованной
    //тоже выдадим
    begin
      inc(m);
      mas1[1,m]:=mas[1,i]; //mas1 - массив в котором накапливаются нарисованные окружности
      mas1[2,m]:=mas[2,i];
      Ellipse(mas[1,i]-R2,mas[2,i]-R2,mas[1,i]+R2,mas[2,i]+R2);
    end;
  end;
end;

procedure TForm1.Button2Click(Sender: TObject);
var
  x,y,i,k,r: Integer;
  a,fi : Extended;
begin
  //Нужно добавить в описание Uses ссылку на модуль Math !!!!!!!!
  R1:=StrToInt(Edit1.Text);
  R2:=StrToInt(Edit2.Text);
  If R2>R1 Then
  Begin
    i:=R1; R1:=R2; R2:=i;
  End;
 
  x0:=R1; y0:=R1; //Центр большой окружности
  PaintBox1.Canvas.Ellipse(0,0,2*R1,2*R1);
 
  r:=R1-R2;
  Repeat
    //Считаем угол поворота...
    If r<R2 Then a:=2*Pi Else
    a:=ArcCos(1-2*Sqr(R2/r));
    //Считаем количество окружностей
    k:=Trunc(2*Pi/a);
    //Рисуем окружности...
    fi:=0;
    With PaintBox1.Canvas Do
    For i:=1 To k Do
    Begin
      x:=Round(r*Sin(fi))+R1;
      y:=Round(r*Cos(fi))+R1;
      Ellipse(x-R2,y-R2,x+R2,y+R2);
 
      fi:=fi+a;
    End;
    r:=r-2*R2;
  Until r<0;
end;