I am working on cole cole model which basically exhibits how the permittivity varies with respect to frequency and is given by;

Where, ε_∞ is the higher permittivity, ε_s is the static permittivity ε_s>ε_∞,
ε_(o) is 8.854e-12,

ω=2πf,

α is the relaxation time 0≪α≪1 σ is the conductivity (can be a CONSTANT) I tried one of the examples given here https://lmfit.github.io/lmfit-py/examples/example_complex_resonator_model.html but was not able to follow along. Currently I have frequency vs real permittivity data. I want to fit the model to my measured data by using lmfit and guess the values from the fit. Basically i want to fit the real part of permittivity to the model.