  We present a general numerical approach to solve the equations for bubble
wall profiles in models with more than one scalar field and CP violating
phases. We discuss the algorithm as well as several problems associated with it
and show some profiles for demonstration found with our method.
