plot sqrt((RiemannSiegelZ[(2*pi*x/ln(2))])^2+((DivisorSum[x, # &]-x)/x)^2), x from 1 to 100