Some new GP features
ramanujantau
The function ramanujantau can now compute the coefficients
of other newforms of level 1, for weight 12, 16, 18, 20, 22, or 26.
? mf=mfinit([1,16,1],1);
? S=mfeigenbasis(mf)[1];
? mfcoef(S,101)
%62 = -817641571654098
? ramanujantau(101,16)
%63 = -817641571654098
? mfcoef(S,10001)
%64 = -870376615320443566571062660188
? ramanujantau(10001,16)
%65 = -870376615320443566571062660188
Here, ramanujantau is much faster than mfcoef.