Hey all,
 I have created a nested while loop to investigate the Riemann zeta function: 

for q to 10 do
 m := 10^(-q);
  for s from 2 to 10 do
   n := 1;
   b := 1;
    while b >= m do
     b := evalf(abs(Zeta(s)-euclid(s, n)));
     n := n+1
    end do;
  print(m, s, n-1)
 end do
end do 

it...