def geometric_progression(a1, n, q): print(a1) for i in range(1, n + 1): power = i - 1 t = a1 * q ** power print(t)