def elem_Q(w, Q, n):
    """
    Simulation Function: -Q-
    
    Inputs
    ----------
    w = Angular frequency [1/s]
    Q = Constant phase element [s^n/ohm]
    n = Constant phase elelment exponent [-]
    """
    return 1 / (Q * (w * 1.0j) ** n)
