def zpn(z, c, n=2):
    """
    nth order polynomial generating a beautiful set

    Parameters
    ----------
    z: Complex
        scalar, vector or matrix to evaluate the function on elementwise
    c: Complex
        scalar, vector or matrix, the starting value(s)
    """
    return z**n + c