Homework 7b numpy

Here you can find the solution of the homework 7b of the 'Herramientas computacionales' course at uniades (Section 6) Course, Homework.

In [1]:
import numpy as np
import matplotlib.pylab as plt
In [2]:
def Leib(n):
    return ((-1)**n)/(2*n+1)
In [3]:
x=np.linspace(0,10,11)
In [4]:
y=Leib(x)
In [5]:
y2=np.cumsum(y)
In [6]:
plt.plot(x,y2*4)
plt.show()
In [7]:
print(y2[-1]*4,np.pi,y2[-1]*4-np.pi)

3.232315809405594 3.141592653589793 0.09072315581580082
In [8]:
x=np.linspace(0,100,101)
In [9]:
y=Leib(x)
In [10]:
y2=np.cumsum(y)
In [11]:
import matplotlib.pylab as plt
plt.plot(x,y2*4)
plt.show()
In [12]:
print(y2[-1]*4,np.pi,y2[-1]*4-np.pi)

3.1514934010709914 3.141592653589793 0.009900747481198291
In [13]:
x=np.linspace(0,1000,1001)
In [14]:
y=Leib(x)
In [15]:
y2=np.cumsum(y)
In [16]:
import matplotlib.pylab as plt
plt.plot(x,y2*4)
plt.show()
In [17]:
print(y2[-1]*4,np.pi,y2[-1]*4-np.pi)

3.1425916543395442 3.141592653589793 0.0009990007497511222