Creo que para este problema la forma más directa es hacer un bucle que implemente exactamente lo que describes. Para ello no necesitas numpy, puedes hacerlo con listas "normales" de Python. No se me ocurre forma de lograr lo mismo con una sola línea de numpy que lo haga vectorialmente.
Pero eso sí, a tu descripción se la puede optimizar un poco, ya que donde dices:
- paso 1: sustraer del valor anterior (x) otro valor (y)
- paso 2: multiplicar el resultado por 0.01
- paso 3: sumar los dos anteriores
Realmente puedes decir:
- paso único: computar (x-y)*1.01
Ya que el paso 1 te da (x-y), el paso 2 te da (x-y)*0.01 y el paso 3 te da (x-y) + (x-y)*0.01, que se puede agrupar como (x-y)(1+0.01) y por tanto es lo mismo que (x-y)*1.01
Con esto en mente, el código sería:
a = [1, 2, 3, 3, 4, 4, 4, 4, 4]
result = [100]
previo = result[0]
for n in a:
previo = (previo - n) * 1.01
result.append(previo)
y al terminar tienes en result
:
[100.0,
99.99,
98.9699,
96.929599,
94.86889499,
91.7775839399,
88.655359779299,
85.501913377092,
82.31693251086291,
79.10010183597154]