Numerical-stochastic procedure developed for an optimization of underground cable systems is presented. The procedure is based on: (1) the finite element mesh generator realized by means of the advancing front method, (2) the non-linear coupled electric-thermal model, (3) the multiquadric approximations of ampacities including the successive zooming steps, (4) the stochastic optimization method - differential evolution.