do 10 i = 1, 200 x = -i z = amillr(x) write (*, '(1x,e16.8, 2x, e16.8)') x, z 10 continue end