  We give a new algorithm of slow continued fraction expansion related to any
real cubic number field as a 2-dimensional version of the Farey map. Using our
algorithm, we can find the generators of dual substitutions (so-called tiling
substitutions) for any stepped surface for any cubic direction.
