from __future__ import print_function import numpy as np N = 100 x = np.zeros(N+3) x[0] = 1. x[1] = 2./3. for n in range(0,N+1): x[n+2] = -(4/3.0)*x[n+1] + (4./3)*x[n] print('x[%d]: %g' % (n, x[n]));