

In [1]:

    
%pylab
%matplotlib inline









    



Using matplotlib backend: MacOSX
Populating the interactive namespace from numpy and matplotlib

Trapes



In [2]:

    
n=4 #velg antall delintervaler
a=0
b=1
Deltax=(b-a)/n
integralet=exp(-a**2) #her skriver vi f(a)
for x in arange(a, b, n):
    integralet=integralet+2*(exp(-x**2))#her skriver vi integralet + 2*f(x)

integralet=integralet + exp(-b**2)#her skriver vi integralet + f(b)
integralet=Deltax*integralet/2 #Gange med Delta x og delt på 2
print(integralet)









    



0.420984930146

Riemann sum

Riemann sum for $f(x)=x^2$ i $[1,2]$:



In [3]:

    
n=4 #velg antall delintervaler
a=1
b=2
Deltax=(b-a)/n
integralet=a**2#her skriver vi f(a)
for x in arange(a, b, n):
    integralet=integralet+x**2#her skriver vi integralet + f(x)

integralet=Deltax*integralet #Gange med Delta x 
print(integralet)

0.5



In [ ]: