def lorenz(x: float, x0: float, fwhm: float, a: float, o: float) -> float:
    """
    Lorenzian function
    Args:
        x: frequency(s) to evaluate function
        x0: center frequency of resonance
        fwhm: full-width-half-max of resonance
        a: amplitude of function
        o: offset of function from 0

    Returns:
        o + a /((x-x0)**2 + (fwhm)**2)
    """
    return o + a / ((x-x0)**2 + (fwhm)**2)