>>469 Yep! sum_{n=0}^k n^a= =(Bernoullipolynomial_(a+1)(k+1)-Bernoullinumber_(a+1))/(a+1)= =harmonicnumber(k,-a)= =-Hurwitzzeta(-a,k+1)+Riemannzeta(-a) for a in Z+