from numpy import zeros, float32, float64

n = 20
x = zeros(n, dtype=float32)
y = zeros(n, dtype=float32)
x[0], y[0] = 1.0, 1.0
x[1], y[1] = 1.0/3.0, 1.0/3.0

for i in range(2,n):
   x[i] = (13.0/3.0)*x[i-1] - (4.0/3.0)*x[i-2]
   y[i] = y[i-1]/3.0
   print "%6d %20.10e %20.10e" % (i, x[i], y[i])

n = 30
x = zeros(n, dtype=float64)
y = zeros(n, dtype=float64)
x[0], y[0] = 1.0, 1.0
x[1], y[1] = 1.0/3.0, 1.0/3.0

for i in range(2,n):
   x[i] = (13.0/3.0)*x[i-1] - (4.0/3.0)*x[i-2]
   y[i] = y[i-1]/3.0
   print "%6d %20.10e %20.10e" % (i, x[i], y[i])