def f(x):              #Her definerer vi en fukjson.   
    return  x**9       # To stjerner er tegnet for potensopphøying i Python   

"""Vi implementerer Simpsons medtode.  Vi estimerer integralet av f(x) fra x=a til X=b 
med m doble intervaller.  m=2m er da antall intervaller vi benytter. Bredden på hvert av de 2m 
interallene er gitt ved d.
"""
a=0                    
b=1
m=1000000
d= (b-a)/(2*m)

I=f(a) - f(b)    #Vi trekker i fra f(b) siden to kopier av f(b) blir lagt til når vi gjører for-løkken. 

for j in range(m):                                         # j gjennomløper 0,1,2,  ...  , n-1
 I = I+  4*f( a+d + 2*d*j   ) + 2*f( a+2*d + 2*d*j   )    
 #Vi oppdaterer summen med bidrag fra punkten nummer 2j+1 og 2j +2
 
S = I*(b-a)/(6*m)  #Her regner vi ut tilnærmingen til integralet.

print(S)   #Resultatet skrives ut.