This paper presents a boundary integral equation method with the adjoint generalized Neumann kernel for computing conformal mapping of unbounded multiply connected regions and its inverse onto several classes of canonical regions. For each canonical region, two integral equations are solved before one can approximate the boundary values of the mapping function. Cauchy's-type integrals are used for computing the mapping function and its inverse for interior points. This method also works for regions with piecewise smooth boundaries. Three examples are given to illustrate the effectiveness of the proposed method.