I am using the below function to compute the Hurwitz zeta function from Riemann zeta function. But I am not getting the correct results when compared with the value of Wolfram alpha Hurwitz zeta function. I am suspecting some issues with floating point errors but not sure. Can someone help me out with the explanation?
function h=hzeta(s,q)
z = zeta(s)
h = z - sum((1:(q-1)).^(-s));
where $q=1-500$ and $s=1-7.$
Thanks,