  We present a method for numerical computation of conformal mappings from
simply or doubly connected domains onto so-called canonical domains, which in
our case are rectangles or annuli. The method is based on conjugate harmonic
functions and properties of quadrilaterals. Several numerical examples are
given.
