default(parisize,"10M"); \\ 8M overflows
tauvec(N) = Vec(q*eta(q + O(q^N))^24);
N = 10^3; vector(N, n, ramanujantau(n)) == tauvec(N)
ramanujantau(10^6+3)
ramanujantau(10^7+3,16)
ramanujantau(10^7+3,24)
ramanujantau(-1)
ramanujantau(factor(-1))
ramanujantau(factor(0))
ramanujantau(factor(6))
ramanujantau(factor(6),24)
ramanujantau(factor(-6),24)

\\ ERRORS
ramanujantau(x)
