Plotting poly logarithmic function
I am trying to plot the below function f(s,t):=(2π)sΓ(s)polylog(−s,e−2πt) in the range (s,t)∈(k,∞)×(1,∞) where k is a fixed positive real number. Thus, for my purposes, it is a real valued function on (subset of ) R2. (Although, the function polylog(s,z) is known to be complex analytic with respect to both s and z for arbitrary complex numbers s and |z|<1.)
Following is the code which is working fine in sage Cell server for smaller values of s.
s,z,t=var('s,z,t')
z(t) = exp(-2*pi.n()*t)
plot3d(((2*pi)^s)/gamma(s)*polylog(-s,z(t)),(s,10,150.1),(t,1,100))
However, the same code fails when I run larger values of s, even 200. I need to study the limiting behaviour of this function in terms of both parameters. Can someone tell me how to work around this "math range error"? PS I am new to Sage.
s, z, t = var('s,z,t')
z(t) = exp(-2*pi.n()*t)
plot3d(((2*pi)^s)/gamma(s)*polylog(-s,z(t)),(s,10,1500.1),(t,1,100))
opt/sagemath-9.0/local/lib/python3.7/site-packages/sage/repl/rich_output/display_manager.py:592: RichReprWarning: Exception in _rich_repr_ while displaying object: math range error
RichReprWarning,
The function is getting very close to zero when s and t gets bigger:
Maybe the exponents get so big that the plot method can not handle it? (that would be a bug)