nb10=N(pi,digits=500)
x=floor(nb10*12^272)
decimales=x.digits(12)

couleur666=[(0,0,0),(1,0,5),(2,0,4),(3,0,3),(4,0,2),(5,0,1),(5,3,1),(4,3,2),(3,3,3),(2,3,4),(1,3,5),(0,3,0)]
couleur=[(couleur666[k][0]*0.2,couleur666[k][1]*0.2,couleur666[k][2]*0.2) for k in range(12)]

liste_points=[[-10+k1,6-k2] for k2 in range(13) for k1 in range(21)]
tableau=sum(point(liste_points[k],rgbcolor=couleur[decimales[272-k]],pointsize=100) for k in range(273))
tableau.show(aspect_ratio=1,axes=False)