Dass du mit der Restkapazität rechnest, ist schon ok. Aber irgendwas ist danach schiefgelaufen. HIer mal als Illustration dein Ansatz. Wenn du zoomst, siehst du, dass 2000 ungehindert überschritten wird.
~plot~200*1,045^x;[[500]]~plot~
5% der Restkapazität müsste sein 0.05 ( 2000 - f(t)) und die werden addiert zu f(t).
Das wäre mit folgender Rekursion zu haben.
f(t+1) = f(t) + 0.05(2000 - f(t))
Nun ist das aber vorläufig noch keine Funktionsgleichung.