  We present a method to derive the analytical expression of the roughness of a
fractal surface whose dynamics is ruled by cellular automata. Starting from the
automata, we write down the the time derivative of the height's average and
variance. By assuming the equiprobability of the surface configurations and
taking the limit of large substrates we find the roughness as a function of
time. As expected, the function behaves as $t^\beta$ when $t\ll t_\times$ and
saturate at $w_s$ when $t\gg t_\times$. We apply the methodology to describe
the etching model \citep{Bernardo}, however, the value of $\beta$ we obtained
are not the one predicted by the KPZ equation and observed in numerical
experiments. That divergence may be due to the equiprobability assumption. We
redefine the roughness with an exponent that compensate the nonuniform
probability generated by the celular automata, resulting in an expression that
perfectly matches the experimental results.
