Here, ω k are the experimental frequency values, g 0 represents the Ohmic resistance (τ 0 = 0), g l for l ≥ 1 is a discretization of the continuous distribution R p τ γ(τ ) with N l=1 g l = R p , and the corresponding τ l are user-selectable, logarithmically spaced time constant values τ l that cover at least the frequency range of the experiment. Equation 3 can be expressed in matrix notation using vector ω containing all experimental frequency values ω k and vector τ containing the τ l values. τ 0 can also be included in τ, albeit g 0 is calculated by only using Tikhonov regularization in standard form (vide infra) . Then, a discretized kernel