For L a given L-function (such as the Ramanujan-tau Dirichlet series) I would like to compute L(s) at 2.5*10^5 values of s equidistributed in a square region of the complex plane in a reasonably short time (meaning, say, less than an hour.) Is there a Maple function that will do this, either in the original software or available on the web?