# 0811.py
execfile('../libreria/alfa.py')

def h (j,k): a=(j+2*k+3*j*k)%41; return a/40.0

w=imm(41,41,h)
w.scrivibyte('0811.pgm')

z=imm()
z.leggibyte('0811.pgm')
print z.colonne, z.righe
# 41 41

for j in xrange(2):
  for k in xrange(2): print h(j,k),z(j,k)
# 0.0 0.0
# 0.05 0.0470588235294
# 0.025 0.0235294117647
# 0.15 0.149019607843