tmin=0
tmax=100
pas=0.001
smin=0
smax=0.95
pas_s=0.05
#.......................................
nb=(tmax-tmin)/pas
i=CC.0
#.......................................
def Zeta(x,y,couleur) :
    z=zeta(x+i*y)
    pz=(z.real(),z.imag())
    return(point(pz,rgbcolor=couleur,pointsize=1))
#.......................................
#for s in srange(smin,smax,pas_s):
#  g=Graphics()
#  for t in srange(tmin,tmax,pas):
#    g+=Zeta(s,t,hue(tmin+t*(1-tmin)/tmax))
#  nom='Graphe_Zeta'+str(s)+'_'+str(tmin)+'_'+ str(tmax)+'_'+ str(floor(1/pas))+'.pdf'
#  g.save(nom)
s=0.95
g=Graphics()
for t in srange(tmin,tmax,pas):
    g+=Zeta(s,t,hue(tmin+t*(1-tmin)/tmax))
nom='Graphe_Zeta'+str(s)+'_'+str(tmin)+'_'+ str(tmax)+'_'+ str(floor(1/pas))+'.pdf'
g.save(nom)