# -*- coding: utf-8 -*-
"""
Created on Fri Mar 31 09:23:11 2023

@author: 2017
"""

#Implementerer Simpsons metode

import numpy as np

def f(x):
    return x**4
#np.exp(x)
#np.sin(x)
#4 * np.sqrt(1-x**2)
#x**9
a=0
b=1
n=1000   #antall doble intervaller

m=2*n  #antall enkle intervaller
d = (b-a)/m   #avstand mellom punktene vi evaluerer funksjonen i.

sum  =  0  

for i in range(n):
   sum = sum + 4* f(a + d +2*d*i)  + 2 * f(a + 2*d*i)
   
sum = sum + f(b) - f(a) 

sum = sum *(1/3)*d

print("Estimatet for integralet er: " , sum )