Рассматривается задача о построении оптимальной упаковки из фиксированного числа n > 1 кругов в общем случае различного радиуса в плоское компактное множество М. Считается, что для каждого элемента упаковки задано положительное число такое, что радиус круга равен его произведению на общий для всей упаковки параметр г. Критерием оптимальности выбран максимум r, что приводит в том числе и к увеличению плотности упаковки - отношения ее площади к площади фигуры М. Основу метода решения задачи составляет итерационное изменение координат центров элементов упаковки Sn, дающее возможность увеличивать радиусы кругов. Разработанные вычислительные процедуры реализуют имитацию отталкивания центра каждого элемента упаковки от близко лежащих других центров и от границы множества М. Исследованы дифференциальные свойства функции двух переменных (х, у), значение которой равно максимальному радиусу круга упаковки, располагающегося с центром в точке (х, у). При это координаты центров остальных элементов упаковки считаются фиксированными. При программной реализации используется конструкция чебышевского центра компактного множества. Создан программный комплекс, с его помощью рассмотрен ряд примеров для множеств М различной геометрии. Выполнена визуализация полученных результатов. We consider the problem of constructing an optimal packing of a fixed number n > 1 of circles, generally, of different radii in a planar compact set M. It is assumed that a positive number is given for each element of the packing such that the radius of the circle equals the product of this number by a parameter r, which is common to the whole package. The optimality criterion is the maximum of r, which leads, in particular, to an increase in the packing density, which is the ratio of the area of the packing to the area of M. In the proposed solution method, we iteratively change the coordinates of the centers of the packing elements Sn, which makes it possible to increase the radii of the circles. The developed computational procedures imitate the repulsion of the center of each element of the packing from nearby centers of other elements and from the boundary of M. We study the differential properties of a function of two variables (x, y) whose value is the maximum radius of the circle of the packing centered at the point (x, y), where the centers of the remaining elements are assumed to be fixed. The software implementation employs the notion of Chebyshev center of a compact set. A software complex is created and a number of examples are considered for sets M of different geometry with the use of this complex. The results are visualized.