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()