In [1]:

    

from sympy.plotting import plot


    

In [2]:

    

i = Integral(log((sin(x)**2+1)*sqrt(x**2+1)),(x,0,y))


    

In [3]:

    

i


    

In [4]:

    

i.evalf(subs={y:1})


    

In [5]:

    

plot(i,(y, 1, 5))


    

In [6]:

    

s = summation(1/x**y,(x,1,oo))


    

In [7]:

    

s


    

In [9]:

    

plot(s, (y, 2, 10))


    

In [10]:

    

p = plot(summation(1/x,(x,1,y)), (y, 2, 10))


    

In [11]:

    

p[0].only_integers = True

p[0].steps = True

p.show()