Write a Python function to return stella octangula numbers

Solution

def stella_octangula_number(n):
    if n >= 0:
        return n*(2**n - 1)