from __future__ import print_function import numpy as np N = 50 x = np.zeros(N) x[0] = 190 x[1] = 260 for i in range(2,N): x[i] = 0.25*x[i-1] + 1.25*x[i-2] print('x[%d]: %g' % (i, x[i]));