from sympy.plotting import plot
i = Integral(log((sin(x)**2+1)*sqrt(x**2+1)),(x,0,y))
i
i.evalf(subs={y:1})
plot(i,(y, 1, 5))
s = summation(1/x**y,(x,1,oo))
s
plot(s, (y, 2, 10))
p = plot(summation(1/x,(x,1,y)), (y, 2, 10))
p[0].only_integers = True
p[0].steps = True
p.show()