In the paper the correspondence between a formal multiple power series and a special type of branched continued fractions, the so-called ‘multidimensional regular C-fractions with independent variables’ is analysed providing with an algorithm based upon the classical algorithm and that enables us to compute from the coefficients of the given formal multiple power series, the coefficients of the corresponding multidimensional regular C-fraction with independent variables. A few numerical experiments show, on the one hand, the efficiency of the proposed algorithm and, on the other, the power and feasibility of the method in order to numerically approximate certain multivariable functions from their formal multiple power series.