from math import * h = 16383 tau = 300 tmax = 4096 y0 = h / (1 - exp(-tmax / tau)) for t in range(tmax): y = y0 * (1 - exp(-t / tau)) print(round(y)) for t in range(tmax): y = h - y0 * (1 - exp(-t / tau)) print(round(y))