  We suggest an approximate method of finding a conformal mapping of an annulus
onto an arbitrary bounded doubly connected polygonal domain. The method is
based on the parametric Loewner--Komatu method. We consider smooth one
parameter families $\mathcal{F}(z,t)$ of conformal mappings of concentric
annuli onto doubly connected polygonal domains $\mathcal{D}(t)$ which are
obtained from a fixed doubly connected polygonal domain $\mathcal{D}$ by making
a finite number of rectilinear slits of variable lengths; meanwhile we do not
require the family of domains $\mathcal{D}(t)$ to be monotonous. The integral
representation for the conformal mappings $\mathcal{F}(z,t)$ has unknown
(accessory) parameters. We find a PDE for this family and then deduce a system
of PDEs describing dynamics of accessory parameters and conformal modulus of
$\mathcal{D}(t)$ when changing the parameter $t$. We note that the right-hand
sides of equations in the system of ODEs involve functions standing for speeds
of movement of the end points of the slits. This allows us to control
completely dynamics of the slits and to seek their agreed change, if we have
more than one slit. We also give some results of numerical calculations which
show the efficiency of the suggested method.
