Electrochemical impedance spectroscopy is attracting more attention due to an increasing production of power sources. One highly popular tool for diagnosing diverse power sources is distribution function of relaxation times (DRT), which has led to numerous approaches for extracting DRT from impedance data. The majority of these are based on the numerical approximation of integral. However, herein we have applied an analytical approximation of the EIS integral. For the first time, we have employed Levenberg-Marquardt algorithm (LMA) to extract the applicable DRT from impedance data by using the Jacobian matrix that was obtained without any discretization errors. Although LMA was previously used to fit EIS data by DRT characteristics, the DRT profile was not applicable due to discretization errors. In this work, LMA was applied as it has an automatic update of the regularization λ parameter. Tests conducted in this work have shown that LMA is capable of extracting DRT from ZARC and FRAC synthetic data.