We present a method for numerical calculation of two dimensional distributions of the attempt relaxation times and activation energies from the temperature dependence of the experimental dielectric permittivity dispersion. We introduce empirical attempts to account for broad and/or asymmetric dispersions with the idea of using a weighted collection of Debye relaxation times. Then we present a modification of the aforementioned idea including attempt relaxation time and activation energy using the Arrhenius law, which significantly complicates the computation of the aforementioned distribution. Incorporating the activation energy and the attempt relaxation time into the equation transforms the discretized matrix equations into tensor equations. We rework the tensor equations into simpler matrix equations, thus permitting us to solve the presented discretized integral equation by using existing Least Distance Problem solving methods. Also, we present a regularization method and a way to choose the regularization parameter based on a best fit criterion. In the end we discuss the method showing some simulated results and experimental results. We then point out some problems involved in the calculations and propose methods to reduce their significance.