def cir_RsQ(w, Rs, Q, n):
    """
    Simulation Function: -Rs-Q-

    Inputs
    ----------
    w = Angular frequency [1/s]
    Rs = Series resistance [Ohm]
    Q = Constant phase element [s^n/ohm]
    n = Constant phase elelment exponent [-]
    """
    return Rs + 1/(Q*(w*1j)**n)